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Abstract 

Cataclysmic Variables (CV) are binary systems where a main-sequence star transfers mass onto a 
white dwarf (WD). According to standard evolutionary theory, angular momentum loss drives CVs 
to initially evolve from longer to shorter orbital periods imtil a minimum period is reached 80 
minutes). At roughly this stage, the donors becomes degenerate, expand in size, and the systems 
move towards longer orbital periods. Theory predicts that 70% of all CVs should have passed their 
minimima period and have sub-stellar donors, but imtU recently, no such systems were known. I 
present one CV showing evidence of harbouring a sub-dwarf donor, SDSS J1507-I-52. Due to the 
system's imusually short orbital period of ^ 65 minutes, and very high space velocity, two origins for 
SDSS J1507+52 have been proposed; either the system was formed from a young WD/brown-dwarf 
binary, or the system is a halo CV. In order to distinguish between these two theories, I present UV 
spectroscopy and find a metallicity consistent with halo origin. 

Systems close to the minimum period are expected to be faint and have low accretion rates. 
Some of these CVs show absorption in their spectra, implying that the underlying WD is exposed. 
This yields a rare opportunity to study the WD in a CV. I introduce two new systems showing WD 
signatures in their light curves and spectra, SDSS J1457-I-51 and BW Sculptoris. Despite the fact that 
CVs close to the mioimimi period should be faint, we find systems that are much too bright for 
their orbital periods. Such a system is T Pyxidis - a recurrent nova with an imusually high accretion 
rate and a photometrically determined period < 2 hours. The system is ~ 2 times brighter than 
any other CV at its period. However, to confirm the status of this unusual star, a more reliable 
period determiaation is needed. Here, I present a spectroscopic study of T Pyxidis confirmiag its 
evolutionary status as a short-period CV. 

In this thesis, I discuss what implications these systems may have on the current understanding 

of CV evolution, and the importance of studying individual systems in general. 
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Chapter 1 

Cataclysmic Variables and their 
Evolution 

In this Chapter, I will give an overview of cataclysmic variables (CVs) and how 
we think they evolve. Key physical concepts such as Roche-lobe overflow, loss of 
angular momentum and nova eruptions will be introduced, as well as characteristic 
observational features in CVs, in particular, their variability. I will then give an 
overview of the current status of theoretical and observational understanding of 
these objects. 

1.1 CVs AS Accreting Systems 

Cataclysmic variables are interacting binary stars in which a donor star loses mass 
to a primary white dwarf (WD). The donor is generally near the main sequence 
(except in rare cases where the donor may already have undergone a significant 
amount of nuclear evolution). In Figure IlTI I present the standard picture of a CV, 
showing how matter from the donor is transferred to the WD by means of an ac- 
cretion disc. In the disc, viscous processes efficiently transport angular momentum 
outward and matter inward towards the WD. As a result, gravitational potential 
energy is converted into radiation, and the disc is responsible for a large portion 
of the detected optical light. The gas emerges from the donor and forms a stream 
that impacts the disc at a location called the bright spot. Close to the WD, kinetic 
energy is converted into radiation as the disc material is slowed down before it is 
accreted by the WD. The boundary layer where this occurs can emit up to half of 
the total luminosity of the system. In a CV, the two stars are very close together, 
and the separation is comparable to the diameter of our Sun. This means that CVs 
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have short orbital periods (on the order of hours). As a cor\sequence of this close 
interaction, the donor is heavily deformed due to tidal and rotational distortion, 
and it is forced to move synchronously with the orbit. 

The interaction between the two stars can be understood by considering the 
Roche-lobe approximation, in which the full gravitational potential of the binary 
is approximated as that of two point masses moving in circular orbits. This is ex- 
pected to be a good approximation, since the donor is forced to rotate synchron- 
ously with the orbit, resulting in the removal of any initial eccentricity. The orbital 
angular momentum of the binary system is given by 



where a is the binary separation, and Mj and M2 are the masses of the primary 
star and donor star, respectively. Figure 11.21 shows a schematic view of the Roche 
potential for a particular mass ratio, q - Mj/M^ - 2. The inner Lagrangian point 
(Li) is the point where the individual Roche lobes of the two stars are in contact. 
Mass transfer will occur if one of the stars overfills its Roche lobe. This can happen 
either if the Roche lobe shrinks in size, or if the star expands its radius. Cataclys- 
mic variables are semi-detached systems, which means that the donor star slightly 
overfills its Roche lobe. The internal gas pressure of the donor then causes it to lose 
mass to the primary through the inner Lagrangian point. 

Roche-lobe overflow in CVs can push the donor star out of thermal equilibrium 
(TE). A star is in TE when the energy generated in the core is equal to the amount 
of energy that is transported to the surface and radiated away. The donor can only 
stay in thermal equilibrium if the time scale for the mass loss, = M2/M2, is 
larger than the thermal time scale, also called Kelvin-Helmholtz time scale, Tth = 
GMj/LjRj- This is defined as the time scale on which the donor is able to adjust 
its structure so as to stay in thermal equilibrium. If the mass-loss rate is very low, 
Tj^^ » Tth, the donor will manage to stay in TE, and its radius will be the same 
as that of a main-sequence star. If the mass-loss rate is very high, « Tth, the 
process of mass loss will be adiabatic, and the donor will expand its radius (as- 
suming that the donor is a low-mass star with a large convective envelope). As it 
turns out, for a CV donor, ~ Tt^, which means that the donor is not quite able 




(1.1) 
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Figure 1.1: Author's impression of a cataclysmic variable. 



to adjust its radius fast enough in response to mass-loss. This means that donors 
in CVs are slightly out of thermal equilibrium, so even though their radii do shrink 
in response to mass loss, this shrinkage does not quite manage to keep up with the 
ongoing mass loss. As a result, CV donors become larger than equal-mass main- 
sequence stars. Observations show that donor stars can be oversized by up to 30% 



( Patterson et al 



2nn5al: 



Kniggell2nn6h 



The radius of the donor (R2) can be compared to the average radius of the 
Roche lobe volume, which only depends on the binary separation and the mass 
ratio q = M2/M1. In the range 0.1 < q < 0.8, this vo lume-averaged radius can be 
approximated to an accuracy of 2% (|Paczynskilll97ll ) by: 



^«0.462p- 



1/3 



(1.2) 



Also, Newton's generalisation of Kepler's third law shows that the orbital period 
(P) depends on the binary separation and the masses of the two components. 



P = 



G(Mi +M2)' 



(1.3) 
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Figure 1.2: Roche-lobe potential for a binary system with a mass ratio q = M2/M1 = 2. If 
one of the stars fills its Roche lobe, material will escape through the inner Lagrangian point 
(here Li), and be accreted by the companion. 



(Source: htt}) ://hemel. wmrnemen.com/Informatie/Sterr en/hoof dstukS.htmlmtr^ 



where G is the gravitational constant. Combining Equations ll.2l and 11.31 yields the 
period-density relation for Roche-lobe filling stars 



P2 = 



3M2 



95.4 1 
~G~P' 



(1.4) 



If we assume that the donor is a low-mass star near the main sequence, M2/M0 = 
/(R2/R0), where f ~\, then the mass and radius of the donor only depend on the 
orbital period of the system. This results in the mass-radius/ period relation 



(1.5) 



In a cataclysmic variable, material from the donor emerges from the Lagrangian 
point at approximately the speed of sound, ~ 10 kms \ and forms an accretion 
stream towards the primary white dwarf. Due to the orbital motion of the system, 
the Lagrangian point moves perpendicular to the line connecting the primary and 
donor, at a speed of about 100 km s"^ 

The stream material has too much angular momentum to accrete directly onto 
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the white dwarf. Instead, the material settles down into the orbit with the low- 
est possible energy for its angular momentum, which is a circular orbit. This is 
the so-called circularisation radius (Rq)- The Keplerian velocity at Rq is - 



VGMTTRc- However, the stream material emerging from the LI point has a ve- 
locity of InRi^i/P. Assuming that the angular momentum of the stream mater- 
ial at the LI radius is conserved at Rq, the circularisation radius is expressed as 
Rq = 4:n^Rl^/GMiP^. Combined with EquationO this yields 

(l+q)Rt, 

Rq = 5 . (1.6) 

Material rotating at Rq loses angular momentum as it spirals down to the WD. 
Closer to the primary, matter orbits faster than further away from it, so any form of 
viscosity acting on the gas rotating near Rq will tend to heat it and spread it into a 
disc. For the total angular momentum to be conserved, a small amount of material 
carrying angular momentum with it must be transferred bac k to the outer parts of 



the disc, causing the disc to extend further out from the WD (|Pringlelll981l) . In the 
outer edge of the disc, angular momentum is fed back to the binary orbit via tidal 
interactions. However, standard molecular viscosity is too inefficient to be relevant 
to CVs; the associated time scale for the transfer of material through the disc would 
be much longer than those observed in CVs. 



Shakura & SunyaevI (|l973h in 



To further our understanding of this problem, 
troduced a so-called a-disc model in which turbulence in the gas is assumed to 
provide the relevant source of viscosity. Their model assumes a geometrically thin 
disc where the disc radius (R^) is much larger than the disc height (H). The kin- 
ematic viscosity (v) at a given radius is taken to be v = av^H, where Vg is the speed 
of sound, and the parameter a is expected to be less than 1, where corresponds to 
no accretion. 

The specific physical origin of this turbulent viscosity is still under debate, but 
one of the most likely sources, at lest for highly ionized discs (suc h as those in nova- 



ike C Vs), is the magneto-rotational instability (re-)discovered by 



Balbus & Hawley 



( 19911) (earlier discussed by Evgeny Velikhov and Subrahmanyan Chandrasekhar 
in 1960). According to this, a weakly magnetised disc, orbiting a central compact 
object will be highly unstable and develop into a turbulent flow, resulting in a re- 
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distribution of the angular momentum. 

The total luminosity of the disc can be expressed as 

(1-7) 

where and Ri is the mass and radius of the WD. However, this is only half of the 
total accretion energy (La^c - GMM/R). The other half of the accretion luminosity 
is released in the so-called boundary layer, which forms the region where the ma- 
terial is transferred from the inner edge of the accretion disc onto the primary WD. 
This energy release occurs because the white dwarf spins at an equatorial speed of 
a few hundred kms"^ (see for instance Chapter |5]l but, just above the WD surface, 
the Keplerian velocity of the material in the disc is about ten times higher, and this 
kinetic energy needs to be dissipated. The dissipation is expected to occur within 
the narrow boundary where material is slowed down and accreted onto the white 



dwarf (e.g.|Pringle|ll98lD. 

If the WD has a strong magnetic field, it will partially or completely prevent the 
formation of an accretion disc, and material will follow the field lines and accrete 
directly onto the magnetic poles of the WD. In this thesis, I will almost exclusively 
discuss CVs which are assumed to have dynamically insignificant magnetic fields. 



1.2 CVs AS Variable Stars 

One fundamental observational characteristic of CVs is their variability. CVs vary 
on many different time scales, from long-term variations of years, down to short- 
term variability of seconds. When studying these events, we are limited by the fact 
that we watch the current CV population from a snapshot in time compared to their 
long-term time scale for evolution. Also, since we are unable to directly resolve the 
light from the different components within the CV, it is quite a challenging task to 
distinguish what mechanisms /locations are responsible for causing any particular 
observed variability. 
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1.2.1 Long-Term Variability 

Classical Nova Eruptions 

As a result of the process of accretion, 
material from the donor star slowly 
builds up on the surface of the WD, 
creating a hydrogen-rich layer that is 
growing in size. The accreted mater- 
ial is heavily compressed due to the 
high gravity of the WD, and electrons 
become degenerate under the intense 
pressure. As the pressure increases, 
so does the temperature. However, 
in a degenerate gas, the pressure does 
not depend on temperature, and excess 
heat can not be dissipated by thermal 
expansion of the gas. The temperat- 
ure will continue to increase until the 
hydrogen layer reaches a critical pres- 
sure (hence a critical envelope mass), 
after which the gas ignites and deton- 
ates in a nuclear-fusion reaction. This 
type of runaway thermo-nuclear reac- 
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Figure 1.3: Schematic figure showing the 
eruption behaviour of Classical novae and 
dwarf novae (from Ritter 1992). 



tion on the surface of the white dwarf is called a classical nova eruption and is 
responsible for the most violent events observed in CVs. In the explosion, a large 
portion of the accreted hydrogen layer is ejected in an expanding shell surround- 
ing the nova. Left in the centre is the WD core, normally composed of carbon and 
oxygen (more massive WDs show evidence of neon/ oxygen cores - so-called neon 
novae). After the eruption, the WD will start the slow process of accreting more 
hydrogen-rich material until enough material is accreted for a new classical nova 
eruption to take place once more. One o f the earliest, but still useful reviews on 



classical nova eruptions was presented by 



Gallagher & Starrfieldl (|l97; 



In theory, the expected mass required to ignite the hydrogen layer on the WD 
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surface should be roughly equal to the total mass of the expelled n ova shell. A WD 



of mass 1 Mq will need to accumulate « 10 ^ Mq, before ignition ( Fujimoto 



1982). 



Assuming an accretion rate of M 1 0^^ Mq yr~^, the time scale between two such 
eruptions is about 1 0^ - 1 0^ years. 

CVs that only have one observed eruption are called Classical Novae (CNe). 
Observationally, CNe with more massive WDs are brighter during the eruptions, 
but decline faster, compared to CNe with lower WD masses. They are therefore 
divided into classes depending on how fast they fade, with the speed class usually 
being defined by the quantity t2, which is defined as the time it takes for the nova 
to decline 2 magnitudes below maximum brightness. In the fastest novae, t2 < W 



days, and they decline with a rate up to 0.2 mag d ^. As. 



and drops with a rate of 0.008 - 0.013 mag d ^ dWarner 



ow no va has f2 > 150 days 



1995ar ). Spectroscopically, 



due to a redistribution of the flux to shorter wavelengths, novae become bluer as 
they decline. After the initial decline in t he visual w avelength range, the flux in 



the ultraviolet region rises (see Figure 1 in lBodell2010l ). The top panel in Figure [L3l 



shows the general behaviour and speed classification of classical nova eruptions. 

There is another subclass of CVs that also undergo classical nova outbursts, 
the Recurrent Novae (RNe). They are observationally defined as novae for which 
more than one outburst has been observed, implying a nova recurrence time scale 
of approximately 10 years - 100 years. According to standard nova models (e.g. 



Yaron et al. 



20050 . this can only occur for systems which have a high accretion rate 
(of order 10"^ Mq y'^) and a massive WD (Mi > 1 Mq). In RNe, the WD is expected 
to gain more mass between eruptions than it ejects during them. This could make 
the already high WD mass exceed the Chandrasekhar limit (the limit for which the 
degeneracy pressure cannot support the WD anymore, Mch ~ 1.4 M©). Therefore, 
RNe are considered as candidate supernova Type la progenitors. 

One intriguing possibility in RNe is that the expanding shells of material from 
subsequent eruptions may interact with each other. For instance, bright knots are 
observed in the shell surrounding T Pyxidis which are t hought to be forme d as 



material from different shells collides and ignite the gas pchaefer et al. 



20101 ). In 



Chapter m I present a spectroscopic analysis of the recurrent nova T Pyxidis. 
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Figure 1.4: The stability of a disc annulus can he understood from the dependence of vis- 
cosity (v) on the surface density (E), which is represented by the so-called S-curve (Figure 
from Connon-Smith 2007). A cool, un-ionized disc would be situated on the stable lower 
branch (A - B), while a hot, ionized disc would be located at the upper stable branch (C - 
D). Systems located between the two stable states (B - D) have thermally unstable discs and 
exhibits dwarf nova eruptions. CVs belonging to the class dwarf novae will move between 
the two states, in the direction indicated by the arrows. 



Dwarf Nova Eruptions 

Dwarf nova eruptions are less violent than the classical nova eruptions described 
in the previous section. However, they occur more frequently and are caused by 
instabilities in the accretion disc. As more material builds up in the disc, the disc 
heats up and eventually a critical density is reached, causing the stored material to 
be rapidly transported onto the WD. As a result, a large amount of gravitational 
potential energy is released. 

Dwarf nova eruptions can be understood in the corttext o f the thermal-v iscous 



disc instability m odel (DIM) (first proposed by lOsakilll974l. see reviews by 



1996 : 



Lasoti 



1 



Osaki 



20011 ). Generally, material from the secondary reaches the disc at a rate 



Msec 3nd will flow through the disc until it is accreted by the primary WD. The rate 
at which material flows through the disc (M) is a function of the viscosity (v). If 
M < Mggj, material will start to build up in the disc. However, if M > Mgec/ gas is 
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removed from the disc. The stability of a disc annulus depends on the viscosity (v), 



and the surface density (E), and can be represented 
a curve is shown in Figure 11.41 (from 



3V a s o called S-curve. Such 



Cormon-Smith 



20070 . In a disc annulus, T 



and M increase with v, which means that the y-axis in Figure [L4] can be expressed 
in any of those units. Each point on the curve indicates the thermal equilibrium 
for a given disc annulus, i.e the locus of points in the E versus M, where viscous 
heating equals radiative cooling. Since this thermal equilibrium is established on 
a very short time scale, disc annuli can effectively only exist with parameters that 
place them on the S-curve. In other words, disc annuli can only move along the 
S-curve in response to changes in M^ec or temporary fluctuations in the physical 
parameters. Any disc annulus for which the externally imposed Msec intersects the 
lower or upper branch of the S-curve is stable. 

Systems for which M^ec intersects the unstable middle branch display frequent 
disc instabilities and are called dwarf novae (DNe). Consider a system starting 
at point A for which Mggc intersects the unstable branch. In such a system, a disc 
annulus is cold and un-ionized. Here M < M^^c, which means that the disc as a 
whole accumulates mass. As more mass is being stored, the surface density of the 
disc annulus increases along with the mass-flow rate, pushing the disc annulus up 
along the path from A to B. Eventually, the critical surface density Ej^ax (point B) 
is reached. Here, M is still less than M^ec ^^id therefore, E still increases. At this 
point, the disc annulus heats up and the gas becomes partially ionised as it moves 
towards the next stable state, which is on the upper branch (point C). Local heating 
in a disc annulus can trigger heating in adjacent annuli, and eventually cause the 
whole disc to heat up globally, resulting in accumulated disc material to be released 
onto the WD. In point C, the gas in the disc annulus becomes fully ionised. Here, 
M > Msec which results in a cooling as the disc annulus moves down along the 
stable branch from C down to D. 

Any disc annulus for which Mgg,, intersects the unstable branch between points 
B and D is unstable, and even if an annulus were to start with M = Mgg^, on this 
branch, any small fluctuation in M^ec or ^ would be sufficient to push the sys- 
tem away from the intersection point, and towards either point B or D. The critical 
surface densities Ej^j^ and "^uiaxf represent the value above which no cold equilib- 
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rium state exists (I^max)/ ^rid the value below which no hot equilibrium state exists 
(^min)- These limits are specific for every s ystem and d epend on, for instance, the 



Lasota 



200l|). 



viscosity parameter and the WD mass (e.g. 

In systems with lower Mggc/ the outbursts are triggered from a location close 
to the WD and are then spread through the disc (so called inside-out outbursts), 
while systems with higher M^ec have so call ed outside-in erup tions that start from 
the outer edge of the disc and move inward jOsaki et al.ll200ll) . 

DNe have typical inter-outburst time scales ranging from weeks to years. The 
eruptions can last a couple of days up to a couple of weeks, and correspond to 
increases in the optical luminosity by 3 - 5 magnitudes. In quiescence, the DNe 
have low accretion rates. DNe dominate the short-period CV population below the 
period gap (P £ 2 hours; see Section [T.3.3ll . They are further divided into subclasses 
depending on their outburst behaviour. A system that belongs to the U Gem class 
of stars displays normal dwarf nova eruptions. Z Cam stars have accretion rates 
high enough to sometim es get caugh t at a more or less constant brightness level on 
the upper stable branch (|Osakilll996|) . SU UMa stars occasionally show superout- 
bursts in-between the normal dwarf nova eruptions. These superoutbursts last 
longer and are thought to be triggered due to a combined tidal-thermal insta bil- 



ity, as the disc reaches the 3:1 resonance radius (jWhitehurst 



1988; Osaki 



1989). A 



particularly interesting and important feature of these superoutbursts, the so-called 
superhumps, is discussed in more detail in Section ll.2.21 Figure [L3] illustrates the 
different outburst behaviours of U Gem, Z Cam and SU UMa stars. 

Whether or not a disc is stable is determined by Mgg^. Systems with a constant 
high Msec will constantly be on the upper stable branch in the S-curve, thus they 
do not exhibit dwarf nova eruptions. Such systems are called nova-like variables 
and include all systems where no eruption has been observed. Nova-like variables 
have hot WDs, high accretion rates and preferentially populate the orbital-period 
region between 3-4 hours. 
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1.2.2 Variability on Orbital Timescales 

The orbital modulation dominates the light curves of eclipsing CVs and provides 
the most accurate orbital period measurements. In non-eclipsing systems, the bright 
spot still causes modulations on the orbital period. The reflection effect caused by 
irradiation of the donor star by the hot WD and disc, can also produce a signal on 
the orbital time scale, although such a signal would roughly be in anti-phase with 
any signal from the bright spot. Also, tidal deformation of the disc can produce 
an ellipsoidal disc shape, causing superhumps, which are detected on orbital times 
scales in systems with very low mass ratio. 

A key marker for observationally identifying orbital modulations is that such 
signals usually are coherent, i.e. stable in frequency and phase. Also, orbital signals 
are always present in quiescence, but might be less obvious during outbursts as the 
orbital signal then might be depleted temperately due to strong disc emission. 



Superhumps 



Superhumps are high-amplitude signals with periods a few percent {~ 3%) longer 



than the orbita' 



VW Hyi (|Vogt 



period and were first detected during a dwarf nova outburst of 



19741) . As described above in Section [1.2.11 dwarf nova outbursts 



are caused by thermal instabilities in the disc. In addition, systems with very low 
mass ratio and low accretion rates show instabilities that are thought to arise from 
an eccentric insta bility at the 3:1 resonance in the tidally unstable accretion disc 



(Whitehurst 



1983). 



More specifically, dwarf novae that display superhumps are called SU UMa 
stars. These systems show both normal dwarf nova outbursts lasting a few days, 
and superoutbursts lasting more than ten days. Generally, superhumps are only 
present during the superoutbursts and are, in most cases, not seen during the nor- 
mal outbursts (see bottom panel in Figure [Oi l . 



To explain the existence of superoutbursts JOsakil (|l989r ) presented the thermal- 
tidal instability model, which combines the thermal disc instability (DIM) and 
the tidally driven eccentric instability. In this model, the disc is circular and tidally 
stable at the early phase of the cycle for the superoutbursts. However, during every 
normal outburst, not all the accreted material is transferred onto the WD, resulting 
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BJD - 24S5000 



Figure 1.5: Kepler light curve of V 344 Lyme showing several normal outbursts and one 



superoutburst (from 



Still et al. 



2010) . 



in an accumulation of mass in the accretion disc. Along with an increased disc 
mass, the total angular momentum of the disc increases, which results in a gradu- 
ally increasing disc radius. Eventually, during a normal outburst, the radius of the 
disc reaches a critical value at the 3:1 resonance with the donor, which causes the 
disc to become tidally unstable. The normal outburst has turned into a superout- 
burst. When the disc even tuall y cools down, the disc becomes circular and tidally 



stable again (see ISmakI (|2009ah for a modification to the theory of superhumps). 
During the superoutbursts, the ellipsoidal disc precesses in relation to the orbital 
motion, and the beat frequency is detected as a superhump signal, with a frequency 
typically a few per cent less than the orbital frequency. The Kepler light curve of 
V344 Lyrae demonstrates a real example o f how a norma l outburst turns into a 
sup eroutburst (see Figure 11.51 obtained from I 



Patterson et al. 



Still et al, 



2010 ). 



(|2005al) have systematically observed more than 200 dwarf novae, 
and find that all systems with mass ratios q < 0.25 exhibit superhumps, while sys- 
tems with q > 0.36, do not. Between these two limiting mass ratios, there is a smooth 
monotonic decline in the fraction of systems showing superhumps. Measurement 



of the superhump period a 



be estimated (e.g. 



Patterson 



ows m ass ratios as well as the component masses to 



1998b[) , since the precession rate, Ppj-g^ = P^J^ - P^^, is 



expected to be proportional to the mass of the donor star, and hence also to the 
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mass ratio. Also, the fractional period excess of sup erh ump period e scales with 



orbital period, e = (Pgh - Porb)/^'orb (e.g. 



Murray 



2000 and 



Pearson 



2007 ). From the 



mass-period relation, e is therefore also a function of q. From the e - g relation, the 



donor r nass can be est i mated by assuming a WD mass; for example. 



Patterson et al. 



( 2005ah : 



Kniggd (|2006h : lKjiiggeet all (|2011h adopt M^d = O.ZSMq, based on obser- 



vational estimates in eclipsing systems. 

In this section, I have only discussed the most common so-called positive su- 
perhumps. In addition, there are also negative superhumps, which produce power 
excess at slightly higher frequency than the orbital frequency. This type of super- 
hump is thought to be associated with a dis c tilt that is undergoing retrog rade pre- 



cession (jHarvey et al. 



1995 : 



Smak 



2009b 



and 



Montgomery & Martin 



20101) . For fur- 



ther information about differen t tvpes of superhumps , see Appendix Al - Hump 



zoology in cataclysmic variables of 



Patterson et al. 



(2002b). 



Although superhumps are most commonly observed during superoutburst, sim- 
ilar signals have also been detected for a few systems in quiescence (for instance 



V344 Lyrae, 
Chapter |3] 



Still et al.1 



20101) . One such system - BW Sculptoris - is discussed in 



1.2.3 Short-Term Variability 



Signals On the White Dwarf Spin Period 



In magnetic systems, variability on time scales of tens of seconds to tens of minutes 
is commonly seen and is usually understood as arising from the accreting magnetic 
poles on the rotating W D. This phenomenon was first detected in the system DQ 



Herculis (jWalker 



1956h and is sometimes referred to as DQ Her modulations. 



Quasi-Periodic Oscillations and Dwarf Nova Oscillations 

Oscillations on similar time scales to the DQ Her modulations have been found in 
some nova-like variables, as well as in dwarf novae during outburst. However, 
neither class of systems are thought to harbour strongly magnetic white dwarfs. 



Th ese osci 
by IWarner 



ations became known as dw arf nova oscillations (DNOs) (reviewed 



2004 



Warner & Woudt 



20081) . DNOs are further distinguished from 



DQ Her modulations by exhibiting considerably lower coherence. The so-called 
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quality factor (Q) is often used as a measure of the coherence of the signal, where 
Q = \ DNOs are characterised by quality factors ranging from 10^ < Q < 10^, 
which ca n be compared to the f actors typical of DQ Her modulations, which are 
Q~10^^ ( Warner &Woudt 200^ . DNOs typically have periods of 5 s - 40 s, where 
the period is somewhat different depending on CV subclass. In some systems, a 
second DNO is detected at a slightly longer period than the first one. As sugges- 
ted by their relatively low quality factors, DNO periods are not constant. In fact, 
they appear to follow a period-luminosity relation, in the sense t hat the minimum 



DNO period corresponds to when the system is at a maximum M (jWarner & Woudt 



20081) . Thus, DNO periods are shorter during the rise of a dwarf nova eruption and 



longer during the declining phase. 

Quasi-periodic oscillations (QPOs) appear at longer periods than t he DNOs 

and are characterised by a yet much lower coherence (typically 5 < Q < 20: IWarner & Woudt 



20081) . QPOs were first discovered as, and n amed for, a broad exc ess of power seen 



during a dwarf nova eruption of RU Pegasi (jPatterson et al.lll977l) . Since then, they 
have been commonly found with periods of 10 - 20 minutes in the high-state light 
curves of many nova-like variables, as well as in dwarf novae during outburst. 

DNOs and QPOs can exist simultaneously. In particular, in systems showing 
double DNOs, the beat period between the two DNOs is roughly equal to the QPO 
period. The ratio ^'qpo/^'dno is constant in systems with high M, indicating that the 
two signals are related an d that they might be cau sed by the same (or related) un- 



derlying physical process (jWarner & Woudt 



20081) . However, QPOs can also exist 



independently of the DNOs and have been detected in dwarf novae during quies- 
cence. 

The origin of DNOs and QPOs is poorly understood, but could be connected 



to the ma gnetic fi e 



magnetic. 



d of t he W D, even for systems tha t are considered to be non- 



Warnei] (|l995bh and 



Warner & Woudtl (|2002l ) constructed what they call 



the Low Inertia Magnetic Accretor model (LIMA), in which the rapid release of 
material onto the WD during outbursts, causes hot, rapidly spinning equatorial 
belts to form on the WD surface. Any magnetic field anchored in these belts would 
create accretion shocks, similar to those seen in intermediate polars. DNOs would 
occur as the equatorial belts spin up and down in response to the varying accretion 
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Figure 1.6: General model showing the deformation of a star caused by g-modes. Here the 
number of modes is m = 6 in the azimuthal direction, and I = 3 in the longitude direction. 
(Source: \http://www.astro.wisc.edu/~townsend/resource/news/nrp-surface-cutaway.pn^ 



rate in the outbur st cycle. The conn ection between DNOs and equatorial belts was 
first discussed by IPaczyhskil (|l978h . As a consequence of the radiation from these 
belts, QPOs are caused by reprocessing in a traveling wave at the inner edge of the 
accretion disc, resulting in a vertical thickening of the disc, that propagate through 
the disc in prograde direction in respect to the motion of the disc material. QPOs 



would then occur as this trave 



primary WD jWarner & Woudt 



ing w ave obscures and reflects the light from the 



2002h . 



Non-Radial Pulsations 

Non-radial pulsations (NRPs) are gravity wave pulsations within the star itself and 
are commonly found in isolated white dwarfs of DA type, so called ZZ Ceti stars. 
These stars have hydrogen-rich atmospheres, and pulsations occur as the WD cools 
and passed through a phase of pulsational instability, detected mainly as g-modes 



( Gianninas et al. 



20061 ). The oscillations themselves are described by three quantum 
numbers, the radial order n, degree /, and the azimuthal number m. In practice, n 
specifies the number of nodes between the centre and surface of the star, m describe 
modes in the azimuth direction, and / in longitude direction. Figure 11.61 shows a 
general model of g-modes in a star where m = 6 and / = 3. 

NRPs are detected with a range of periods from a few seconds to a couple of 
minutes in isolated WDs. During the last decade, similar signals, generally inter- 
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Figure 1.7: Power spectrum showing the observed NRPs in GW Lib as presented by 
Warner & van Zyl (1998). GW Lib was the first CV to be found to harbour a non-radially 
pulsating WD. 



preted as non-radial WD pulsations, have also been detected in faint cataclysmic 
vari ables. The first CV propo sed to harbour a pulsating white dwarf was GW Lib- 



rae. 



Warner & van Zyll (|l998h . found rapid, periodic, and non-commensurate sig- 



nals in its light curve, suggesting non-radial pulsations of the underlying white 
dwarf (see Figure [L7|) . In most CVs, the accretion energy tends to dominate the 
luminosity, and the white dwarf itself, shining with My ~ 10 - 12, is seldom seen. 
However, for some of the most intrinsically faint CVs, spectroscopy and time-series 
photometry does reveal signatures of the underlying white dwarf, such as broad 
absorption features in the spectrum, sharp eclipses and, sometimes, non-radial 
pulsations in the light curve. These signals have now been detected in about a 
dozen CVs, all quiescent systems of low lum inosity. In this thesi s, th ese systems are 



Szkody et al. 



12010) and 



Mukadam et al. 



called GW Lib stars, after the first discovery. 
( 20091) present recent reviews of this group of stars. 

WDs in CVs are different from isolated ones since they undergo accretion, which 
gives them atmospheres of approximately solar composition. The white dwarfs in 
CVs are hotter and are also found to be spinning faster compared to isolated ones 



( Szkody et al. 



20091) . In isolated WDs with pure hydrogen atmospheres, pulsations 



are only observed in stars with temperatures located within a so-called instability 
strip in the log^ - Tgff plane, spa nning the temperature range Tgff = 10900 K - 12200 
K (see Figure [L8] obtained from iGiarminas et al.ll2006l) . Howe ver, there is no clea r 



instability strip for the GW Lib stars (see Figure [L9] taken from lSzkody et al 



20101), 
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Figure 1.8: The instability strip for isolated WDs of type DA, where pulsating WDs are 
mar ked as black and red op en circles. Solid dots represent non-pulsating DA stars (figure 



from 



Gianninas et al. 



2006) 



and pulsations are found in systems with WD effective temperatures up to at least 
15000 K. 



NRPs typically have periods between 80 s - 1300 s (IWarner & Woudtl200m . with 
a very low coherence in the signal. Since their periods are of the same order as both 
DNOs, QPOs and DQ Her modulations, one has be careful when identifying the 
origin of signals found in a light curve. An important indicator is the coherence of 
the signal, which is specific for each type of modulation. 



Flickering 

Low-amplitude stochastic variability on the order of 0.1 magnitudes on time scales 
of seconds up to minutes is seen in almost all CV light curves. This variability 
is called flickering and is associated with the process of accretion. Though the 
physical origin of flickering is not well understood, it may arise preferentially in 
the bright-spot region, where the accretion stream i mpacts the d isc and/ or in the 
boimdary layer in the innermost region of the disc (|Bruchll2000l ). Some faint sys- 
tems with low accretion rates and faint discs show very little flickering (e.g. SDSS 
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Figure 1.9: The Tf,fj- logg distribution ofpul satinz (marked as d ots) and non-pulsating 



(marked as triangles) WDs in CVs (figure fi'om lSzkody et al. 



2om . 



J1457+51 and BW Sculptoris, see ChapterlD. 

1.3 CVs AS Close Binaries: Secular Evolution 
1.3.1 Pre-CV Evolution 

The evolution leading to the formation of a cataclysmic variable starts with a binary 
system containing two main-sequence stars. The separation between the two stars 



is about fl > 50 Rq (e.g. 



Warner 



1995ah . but for the system to become a cataclysmic 



variable, a substantial amount of angular momentum needs to be removed to de- 
crease the size of the orbit. The more massive of the two stars will be the first one 
to evolve into a cool red giant and fill its Roche lobe. Due to Roche-lobe overflow, 
the giant then starts to lose mass to its companion. The centre of mass is therefore 
moved closer to the companion, resulting in a decreased size of the giant's Roche 
lobe. However, it cannot decrease its radius to stabilise the mass transfer, since the 
equilibrium radius of a giant d oes not deperi d on its total mass, but only on the 



mass of its degenerate core (e.g. 



Warner 



1995a|). This will lead to a run-away mass- 



transfer phase. During this phase, gas is transferred faster than it can be accepted 
by the companion, and eventually the accreted gas will overflow the Roche lobe 
of the companion and form a common envelope (CE) enclosing both stars. The 
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binary can stay in this phase for about 10^ - 10^ years. As the stars spiral around 
their centre of mass inside the envelope, they lose orbital energy due to frictional 
interactions with the envelope material. As a result, the orbital separation shrinks 
drastically. The CE phase is expected to end when the deposited energy exceeds 



the envelope's own binding energy (e.g. IWebbinldll984l '). At this point, the envel- 
ope is ejected and forms a planetary nebula. Left in the now close binary system 
are a primary WD and a main sequence secondary. This second detached phase 
can last 10^ - 10^ years until the loss of angular momentum due to gravitational ra- 
diation or magnetic braking (see below) has shrunk the orbit until the point where 
the secondary fills its Roche lobe. This marks the beginning of a new phase of mass 
transfer and the start of the system's life as a stable, semi-detached CV. 

1.3.2 Mechanisms for Angular-Momentum loss 

In general, the driving mechanism for the evolution of cataclysmic variables is loss 
of angular momentum through the processes of gravitation wave emission (GR) 
and/ or magnetic braking (MB). 

Gravitational Radiation 

For systems with periods below ~ 2 ho urs, the main mechanism operating is usu- 



ally taken to be gravitation al radiation ( IKraft et al.lll962l:lPaczynski & Sienkiewicz 



1981 



Rappaport et al.lll982n . The rate at which angular momentum is lost depends 



on the binary separation a and the masses of the two components Mj and M2, 



/ _ 32G3 MiM2(Mi+M2) 

7 ~ ^5 • 



(1.8) 



The corresponding ma ss-loss rate isM«10 -10 



-10 



line with observations ( Patterson 



2009) 



1998b 



Littlefair et al. 



Vli-:^ V ^, which is broadly in 



2008 



Townsley & Gansicke 



Magnetic Braking 

The observationally infer red accretion rates for systems with periods above 3 hours 

(M « IQ-''- 10-^ Moy-i 



Patterson 



1984; 



Townsley & Gansicke 



2009 



Knigge et al. 
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20111) are much higher than can be accounted for by GR. Therefore, another, more 



efficient mechanism for loss of angular mom entum is needed to expla i n the high 



accretion rates in systems with P > 3 hours. 



Verbunt & Zwaan (1981 



(19811); 



Paczynski 



Rappaport et alj (|1983l) . proposed mass transfer driven by a mechanism called 
magnetic braking. In this scenario the donor has a magnetic field, as well as a 
weak, ionised wind that travels along the field lines out to the Alfven radius. Up 
to this radius, the magnetic pressure in the wind material exceeds the gas pres- 
sure. The wind is thus forced to co-rotate with the star out to this radius, produ- 
cing a magnetic braking torque on the secondary star. Due to this, the rotational 
speed of the secondary is reduced. Also, in a close binary system, tidal effects force 
the donor to rotate synchronously with the binary orbit. As a result, the magnetic 
braking torque ultimately removes orbital angular momentum from the system. 
This leads to a decreased orbital period. As long as the process of magnetic brak- 
ing continues, the orbit continues to shrink, and the secondary loses mass steadily. 
Magnetic braking in CVs is now fairly widely accepted as the main mechanism for 
angular-momentum loss in systems with periods above 3 hours. 



1.3.3 The Period Gap 

There is a dearth of systems with orbital periods between « 2 and 3 hours. Ac- 
cording to the standard evolutionary model for CVs, this is explained by assuming 
that CVs evolve through this period range as detached systems. The mass-period 
relation for main-sequence donor stars (Equation I1.5L shows that for an orbital 
period of 3 hours, the mass of the donor is ~ 0.3 Mq. At this stage, the secondary 
becomes fully convective. Since the magnetic field in a low-mass main-sequence 
star is thought to be anchored at the interface between the convective envelope 
and radiative core, the standard model for CV evolution posits that the field will 
disappear (completely or partially) as the star becomes fully convective. With a 
weakened or non-existent magnetic field, the process of magnetic braking is dis- 
rupted or considerably reduced, which means that the accretion rate (M) drops 
substantially. As a result, the donor adjust its radius closer to the one it would 
have for thermal equilibrium. The donor will therefore shrink and detach from the 
Roche lobe, leading to an interruption of mass transfer. The CV will stay in this 
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phase for approximately 10^ years (e.g. lKnigge et al.ll201ll) . This is long enough for 
the secondary to get back to its equilibrium radius. However, angular-momentum 
loss from the process of gravitational radiation is still active, shrinking the binary 
orbit, and so the donor is able to catch up with the Roche lobe radius again at an 
orbital period of about ~ 2 hours. Mass transfer is once again resumed, now driven 
by gravitational r adiation. Emp irically, this period gap spans from the range 2.15 < 
Park < 3.18 hours (IKniggJbooeL . 

The disrupted mag netic-braking theory just outli n ed has been discussed in more 



detai 



and /or used by 



(1993!) and 



Rappaport et al. 



Howell et al. 



200 ih . 



(1982, 



19831), 



Spruit & Ritteil (|l98: 



3), 



Kolb 



1.3.4 The Minimum Period 

As the donor loses even more mass, the thermal time scale increases more quickly 
than the mass-loss time scale, so that eventually Xj^^ < T^h, and the donor is driven 
further and further out of thermal equilibrium again. Adopting the power-law 
mass-radius relation Rj oc (where a ~ 0.8 corresponds to low-mass main se- 
quence stars in thermal equilibrium), and combining this with the period-density 
relation oc M2R2^, gives oc Mj"^^. Differentiating logarithmically yields 



d(ln P[t]) l3a-l\d(\nM2[t]) 



dt 



dt 



(1.9) 



and we obtain a relation between the orbital period-derivative and the total mass- 
loss rate from the donor (mass loss both from wind and mass transfer). 



P _ 3a - 1 M2 
P ~ 2 Mj' 



(1.10) 



This shows that the system reaches a period minimum, P = 0, when it has been 



pushed so far out of thermal equilibrium that t he mass radius index (eva 



along the CV evolution track) becomes a = 1/3 (Knigge 



2006 



Knigge et al 



The minimum per iod is theoretically predicted to occur at P « 65 min (jKolb 



Howell et al 



uated 



2011). 



1993 : 



200 ih . The latest point at which a system can reach its minimum 



period is roughly when the donor mass falls below the hydrogen-burning limit 
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and the donor becomes a sub-stellar object. This is because a sub-stellar object re 
sponds differently to continued mass-loss and expands or stays cor istant in size 



(a £ 0). The hydrogen-burning limit for isolated objects is « 0.072 M© (|Baraffe et al, 
1994 



Kolbl (|1993l ^ and 



Howell et al. 



19971 ^ theoretically predicted that ^ 99% of all 



CVs should have periods below the period gap and that ^ 70% should have passed 
their minimum period (these systems are often referred to as period bouncers). 
The reason for these high percentages is the drastically longer evolutionary time 
scale as systems evolve closer to the minimum period. As a result, a larger por- 
tion of systems should be found close to Pmm- This pile-up of short -period sys - 



tem s is often called the p eriod spike and is predicted, for instance, by 



and 



Kolbl (|l993l ) 



Kolb & Baraffd (|l999h (a further discussion is presented below in Section [T.3.5ll . 
Observationally, systems close to Prnin typically have M = (3 - 10) • 10"^^ MqY'^ 



( Patterson 



1998bl: 



Littlefair et al.„2006h 



1.3.5 Theory versus Observations 

Observations of cataclysmic variables are broadly consistent with the standard evol- 
utionary theory described in the previous sections. Yet there are some observa- 
tional facts than cannot be reproduced by standard theory. The most striking dif- 
ference is t he minimum period. Theoretically, it is predicted to be at « 65 minutes 
dKolbll993|). but obse rvations indicate a minimum period of Pmin = 82.4 + 0.7 minutes 



dGansicke et al, 



20091) . If one allows for tidal distortion and an inflated donor radius 



associated with magnetic activity, the theoretically predicted pe riod becomes some - 



what longer, but is still inconsistent with the observed v alue ( tCnigge et al 



Also, according to the theoretical predictions of 



Kolb 



19931) and 



Howell et al. 



mm 



19971), 



only about 1% of the observed CVs should have periods above the period gap. This 



is in confl ict with observations that show a far larger 



Pretorius et al. 



2007a 



Pretorius & Kniggell2008bh . 



popu 



ation of long-period sys- 



tems (e.g. 

The two main tracers of the evolutionary state of a CV are the orbital period 
and the donor mass. Luckily, it is often a rather simple task to find the orbital 
period of a system (see Chapter |2]|. To a first approximation, the orbital period is 
enough to tell a younger (long-period) system apart from an older (short-period) 
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system. In Figure 11.101 I present the current period distribution o f non- 



Ritter & Kolbl ({201(h . version 



magnetic 



CVs. Periods were obtained from the catalogue by 
7.14. The white area shows all non-magnetic systems while the dark grey area 
only shows the distribution of dwarf novae. An accumulation of systems close to 
the period minimum is seen. This is the period-spike predicted by theory. Earlier 



Gansicke et al. 



versio ns of the Ritter & Kolb catalogue did not show this peak, and 
( 20091) were the first to observationally confirm an accumulation of systems near the 
period minimum using the Sloan Digital Sky Survey (SDSS) CV sample. The reason 
why such a pile-up of systems close to the period minimum has not been observed 
before is simply because these short period systems are difficult to detect due to 
their low accretion rates and faint luminosities. Therefore short-period systems 
tend to be under-represented among observed CVs, although more recent and deep 
surveys (such as the SDSS) are beginning to uncover these systems. 

As already noted in Section 11.3.41 standard evolutionary theory predicts that 
70% of all CVs should have passed their minimum period and have sub-stellar 



donors ( Kolb 



1993 



Howell et al. 



19971) . However, few such systems have been suc- 



cessfully confirmed, and, until recently, almost no CVs containing donors with 
masses below the hydrogen-burning limit had been found. Nevertheless , a few 



CVs containing su b-stellar donors have now been convincingly identified by 



( 2006 . 



2007, 



Littlefair et al. 



20081) . The system discussed in Chapter |5l SDSS J1507+52, is one of 



these CVs, and turns out to be quite unusual in other respects as well. 

In order to distinguish between short-period systems that have not yet reached 
their minimum period and period bouncers, one has to consider both their orbital 
periods and donor masses. According to Equation 11.51 there is a relation between 
the orbital period and the mass of the donor. However, this is only valid for 
systems close to thermal equilibrium (main-sequence stars), whereas sub-stellar 
donors will follow quite a different period-mass relation. To account also for sub- 
stellar donors, the more general power-law mass-period relation can be employed 

oc Mj""^ (see Section ri.3.4L where main-sequence donors have a ~ 0.8, and 
sub-stellar donors have a ~ 0. A plot showing donor masses as a function of their 
orbital periods should then reveal the two branches of CV ev olution, and thi s mass - 
period relation for CV donors has been studied in detail by 



Littlefair et al. 



( 20081) : 
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Figure 1.10: The CV distribution from the Ritter & Kolb catalogue, version 7.14 (2010). 
The white region is the distribution of all non-magnetic systems, while the dark grey area 
only includes dwarf novae. The period gap is marked as a transparent area between 2.14 - 
3.18 hours. The two solid lines represent the theoretical minimum period at 65 minutes in 
grey (Kolb 1993), and the observational one at 82 minutes in black (Gansicke et al. 2009). 



Sirotkin & Kiinl (|2ninh : 



such relation is presented (Figure from 



Knigge et alJ OOllh: IPattersqnl OOllh . In Figure IlTTI one 



Knigge et al. 



201 ih , where the different col- 



ours and shapes of the data points indicate the type of system and how the masses 
were found. The solid line is the best fit to the data, while the dashed line shows 
the standard evolutionary model. The flat line between « 2.3 h - 3.2 h shows how 
the donor mass stays constant through the period gap. 

Figure IlTTI demonstrates the discrepancy between the theoretical and observa- 
tional minimum period. Many studies have been devoted to solve this problem, as 
well as to find a solution that would account for the fact that we have not found as 
many post-bounce systems as predicted by theory. Generally, an enhanced (i.e. 
stronger than GR) angular momentum loss rate below the gap and /or a larger 
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Figure 1.11: Figure from Knigge, Patterson & Baraffe (2010) showing the mass-period 
relation, both from observations and theory. The different colours and shapes correspond to 
different type of systems and how their masses were estimated. Filled squares correspond 
to short-period CVs while filled circles are long-period CVs. Black points represent systems 
where the masses were estimated from superhump-periods while the masses of the red points 
where found from eclipse analysis. 
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donor radius is proposed to solve the problem with the min imum period 



1998b 



Meyer & Meyer-Hofmeisteii2000 



Knigge et al.1120111 '). For example. 



Patterson 




. Littlefair et al. 



(12003) discussed the possibility that stellar activity/ star spots may significantly in- 
flate the donor radii in CVs, which would result in a minimum period closer to 
the observed one. Conversely iKnigge et alJ pOllh . discussed the implications of 
residual magnetic braking below the gap. They found that by allowing for a lar- 
ger donor radius (due to tidal and rotational deformation), increasing the effect of 
angular-momentum loss below the gap by a factor of about 2.5, and by lowering 
the effect of angular momentum loss above the gap (by a factor of 1.5), they were 
able to account for the main differences between theory and observations. Also, 
their model yields a minimum period consistent with observations. Since a lower 
level of angular-momentum loss above the gap would make systems stay above 
the gap for a longer time, their model are also able to explain why we find more 
systems above the period gap than predicted by theory. Although Knigge et al.'s 
(2011) model would require magnetic braking in fully convective donor stars, this 
is not necessarily in conflict with the disrupted magnetic braking model usually 
employed to explain the existence of the period gap. This is because a reduction 
in / by a factor of abou t 5 is sufficient to e nsure that mass transfer will cease at the 



upper edge of the gap (jKnigge et al 



20111) 
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1.4 CVs IN This Thesis 



Statis t ical analyses of observed CV s amples, such as those presented b y 



12008); 



Townsley &: Gansickel(|2009h : 



Gansicke et al. 



120091); 



3 



Patterso nl d2011h and 



Littlefair et al 



Knigge et al. 



( 201lh , are important and powerful tools when studying cataclysmic variables. 
However, to test proposed models quantitatively, it is valuable to also study in- 
dividual CVs. Systems that appear to be outliers could, in fact, be more important 
to CV evolution than initially thought. For instance, there could be systems that 
represent stages of rapid (and potentially destructive) evolution - as believed is the 
case for T Pyxidis presented below. 

In order to establish the evolutionary status of a system, the orbital period, mass 
ratio and component masses have to be accurately determined. Chapter|2]presents 
an overview of some methods that are commonly used to determine some the most 
important system parameters. 

In Chapters |3l|4]and|5l I present analyses of four systems close to the minimum 
orbital period that all are particularly interesting in the context of studying CV 
evolution. Two of them, SDSS J1457+51 and BW Sculptoris, are consistent with be- 
longing to the ordinary CV population, as described by the standard evolutionary 
models proposed for cataclysmic variables. The other two systems, T Pyxidis and 
SDSS J1507+52, are more exotic. 



1.4.1 Two Ordinary CVs 



SDSS J1457+51 and BW Sculptoris are two faint CVs with low accretion rates and 
thin discs. Both have short orbital periods « 80 minutes (Uthas et al. in press), 
which is right at the mean observed minimum period. Absorption from hydrogen 
and heli um is present in bo th their spectra, indicating that the underlying WD is 



exposed (jSzkody et al. 



20051) 



In Chapter|3j I present the discovery of non-coherent pulsations with main peri- 
ods of about 10 and 20 minutes in both systems, and discuss the possible origin 
of these signals. Chapter |3] also discusses an interesting signal detected in BW 
Sculptoris, which would be connected to the unusual phenomenon of quiescent 
superhumps. 
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1.4.2 Two Peculiar CVs 



The system discussed in Chapter IH T Pyxidis, is a luminous RN that accretes at a 
much higher rate than is expected for its ph otometrically estimated orbital period 
of just under 2 hours ( Patterson et al. 1998I) . T Pyxidis used to have a recurrence 
time scale between eruptions of about 20 years. However, the last eruption was in 
1966, which means that the system has passed its mean recurrence time by more 
than 20 years. 

According to standard evolutionary theory, a system with a period less than 2 
hours should be faint and have a low accretion rate. However, T Pyxidis is about 
two times more luminous than expected for a CV driven by GR at this period. 
In order to explain its high l uminosity, it has been suggested that it c ould be a 



1998 



Knigge et al 



200flh. Theoret - 



wind-driven supersoft source ([Patterson et al. . 
ically, RNe both have high accretion rates and high WD masses (|Yaron et al.ll2005l) . 
This implies that RNe, in general, are candidate progenitors for Supernova type 
la. Also, the current evolutionary timescale of T Pyx is only a few million years. 
If other CVs would go through similar phases of high accretion rates, this could 
speed up of their evolution, resulting in a reduction of the number of short-period 
systems. 

Compared to other short-period systems, the high luminosity and accretion rate 
in T Pyxidis would be highly unusual. However, a photometrically established 
period of a non-eclipsing system is never as reliable as a spectroscopically determ- 
ined one. Therefore, to confirm the current status of T Pyxidis as a short-period 
system, spectroscopic data is highly desirable. In Chapter IH I present a detailed 
analysis of such data of T Pyxidis, revealing its current evolutionary status. 



The eclipsing system SDSS J1507+52, which is the subject of Chapter |5l was first 
observed by the SPS S and quickly recognised to be odd due of its short orbital 
period of 67 minutes (jSzkody et al.l. l2005l ). Hydrogen absorption originat ing from 



the primary WD was detected in its optical spectrum (jSzkody et al. 



2005h . and to- 



gether with the short orbital period, this indicates a system with a low-mass sec- 
ondary, a low accretion rate and a faint disc. 



iLittlefair et al. 



(|2007l) performed eclipse analysis of SDSS J1507+52 and found a 
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donor mass of « 0.05 Mq, which would imply that the secondary is a sub-stellar ob- 
ject. However, if SDSS J1507+52 is placed on the plot presented in Figure IlTTI the 
system does not appear on either the evolutionary branch for perio d bouncers or 



Littlefair et al. 



the on e for main-sequence donors. In addition to the low donor mass, 
(|2007n found that its donor radius is much smaller than expected for sub-stellar 
donors at this mass. They also point out that these unusual donor properties may 
indicate a young system, in which the secondary has not yet had time to adjust its 
radius in response to mass loss. They therefore propose that mass transfer started 
only recently in SDSS J1507+52, i.e. that its progenitor was a detached WD-brown 
dwarf binary system. 



Patterson et al 



120081), at 



Another study of SDSS J1507+52 was carried out by 
approximately the same time. They found another peculiarity of the system - it has 
a very high space velocity, much like stars in the galactic halo. If SDSS J1507+52 is a 
halo CV, it should have a low m etallicity and therefore also a smaller donor radius 
and shorter orbital period. Thus 



Patterson et al. 



20081) suggest that the system is a 



member of the Galactic halo population. 

The effective te mperature for SDSS 11507 +52 was estima t ed to be « 11000 K 



- 11500 K by both 



Littlefair et al. 



(|2007h and 



Patterson et al 



(|2008|). At this low 



temperature, strong absorption features due to Fe II and III are expected in the UV 
spectrum. A simple test to distinguish between the two proposed theories would 
thus be to obtain UV spectra to find out if this absorption is present. In Chapter |5l 
I present the analysis of such a study of SDSS J1507+52, which settles the question 
of its origin as either a young WD-brown dwarf system, or a halo CV. 
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Methods and Techniques 

- FOR ANALYSIS OF CVS - 



The first step in confirming the current evolutionary status of a CV is to determine 
accurately its key system parameters, such as the orbital period and mass-ratio. In 
this chapter, I will discuss what methods are commonly used to determine some of 
these key systems parameters. I will also give a brief summary of the observational 
properties of CVs and of some common techniques used to analyse such observa- 
tional data. In particular, I will focus on those techniques and methods applied in 
Chapters IS landlS] 

2.1 Observational Properties of CVs 

Since telescopes are not able to separate spatially the light originating from the 

accretion disc from that coming from the WD or donor, a major part of studying 

CVs involves the development of methods to isolate the light originating from the 

different components of the system. Eclipsing systems, in particular, provide us 

with a unique opportunity for this. Figure 12.11 shows the optical light curve of 

the eclipsing system SDSS J1507+52 (see Chapter |5]l, together with a sketch of the 

phase-dependent CV geometry that is thought to give rise to these light curveill 

The light curve for this particular system shows a large hump during phases where 

the bright spot is visible. The very steep ingress and egress of the eclipse indicates 

a faint accretion disc. This all implies that the bright spot and WD dominates the 

flux distribution in SDSS J1507+52 (a detailed eclipse analysis of this system was 

^Both the light curve (obtained at the Nordic Optical Telescope during 2005) and phase models of SDSS 
J1507+52 are provided by Are Vidar Boye Hansen, Olesja Smirnova and Arturs Barzdis. 
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/ 



m m 




Figure 2.1: The optical light curve of SDSS J1507+52 obtained at the Nordic Optical 
Telescope, together with models showing the phase-dependent CV geometry (adapted from 
figures by Are Vidar Boye Hansen, Olesja Smirnova and Arturs Barzdis). 



presented by 



Littlefair et al 



2007 ) 



In addition to giving the most accurate orbital periods, analyses of eclipsing 
systems can provide us with a number of other system-parameter estimates. For 
instance, as shown in Figure IZTl the total duration of the overall eclipse is propor- 
tional to the size of the accretion disc. The size of the WD is found from the duration 
of the sharp ingress and egress of the eclipse. Also, the depth and duratio n of the 



eclips e is a function of the mass ratio and inclination (e.g Equation 2.64 in 



Warner 



1995ah . 



No eclipse analyses have been applicable to the studies I have performed in this 
thesis, and as such, no detailed deduction of this method is given. 



2.1.1 Multi-Wavelength Properties 

The various components in a CV have a different relative contribution to the radi- 
ation in different wavelength regions. Figure IZ2l shows an example o f the spectral 
energ y distribution (SED) of a CV with a high-state accretion disc (from lPringle & Wade 
19851) . The following are common to all CVs: 



- Accretion disc: The flux distribution from the accretion disc spans the whole 
range from IR to UV. It is usually the dominant feature in the optical wavelength 
region, making it difficult to observe any of the stellar components directly. The 
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Figure 2.2: The spectral energy distribution of a CV, comprising a high-state accretion disc 
(Pringle & Wade 1985). 

temperature in the disc varies from ~ 100000 K at the inner edge of the disc to ~ 3000 
K in the outer edge of the disc. The disc spectrum mostly contains emission lines 
from low ionisation states, such as H (from Paschen, Balmer and Lyman) and He I. 
A higher state of ionisation is often detected in systems with high accretion rates. 
Especially, He II is used as a tracer of high M s ystems. For instance , the system 



Pyxid is has an accretion rate of > 1 ^ Mq yr ^ ([Patterson et al. 



1998 



Selvelli et al, 



20081) ■ and shows very strong He II lines in its spectrum (see Figure [ 

Emission lines formed in the accretion disc exhibit double-peaked profiles due 
to strong rotational Doppler broadening. The splitting of an emission line will be 
more or less pronounced depending on the or bital inclinat ion and the spectral res- 



olution of the data (quantitatively explored by ISmakll 19811) . Furthermore, optically 
thick disc lines (as considered to be the usual case for disc emission of H and He II), 
exhibit a deeper central minimum for higher inclinations. For optically thin lines, 
the central minimum is shallower and ha s a softer shape. The mechanism of line 
formation is described in detail by Horne & Marshlll986, and is illustrated in Fig- 
ure 12.31 In the left figure (A), the radial velocity contour geometry of a Keplerian 
disc (seen face-on) is shown, while the figure to the right (B), shows its correspond- 
ing emission line profile viewed at an angle i. Keplerian rotation causes the lines of 
constant velocity to form a dipolar pattern aligned with the line of sight. The wings 
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2 
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Figure 2.3: Emission from different regions in a Keplerian disc are marked with i ts corres 



ponding velocity bin in the line profile (Figure adapted from lHorne & MarsH\l98a. showing 
a system with q = 0.15). 

of the line profile are shaped by gas rotating at high radial velocities close to the WD 
(see position 1 in Figure 1231) . The radiating surface area A in each velocity bin dV 
determines the strength of the line at each velocity V, and is proportional to V^^dV 
(since A oc RdR and R oc V^'^ => dR oc V^^dV), which explains the rapid decline in 
the line wings. The maximum line strength is reached when = +VKep(^D)sinf, 
which corresponds to the outer boundary R^ of the emission line region (position 
2). At this point, A starts decreasing, reaching a local minimum at V - where the 
flux arises from gas moving perpendicular to the line of sight (position 3), and a 
dip in flux, and apparent splitting of the line is seen. 

Since the emissivity in the disc is not constant (the inner parts of the disc is 
much hotter than the outer parts), the line strength does not only depend on the 
radiating surface area, but also on the disc emissivity. 



- White dwarf: A WD spectrum displays broad absorption lines (usually of hydro- 
gen), widened by the process of pressure broadening. Average WD temperatures 



j ust above the 



( Urban & Sion 



perio d gap are ~ 26000 K and just below the period gap ~ 18000 K 



200a) . These hig h temperatures imply that the WD is visible in the 



ultraviolet (UV) spectrum fe.g. lTownsley & Bildsten 2002). In short-period system 
with low M, observations in the UV can be used to isolate the WD light from the 
accretion disc and bright spot. 
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Figure 2.4: Optical spectrum of SDSS J1507+52 by Patterson, Thorstensen & Knigge 
(2008), showing deep double-peaked disc emission lines, superimposed on broad absorption 
lines originating from the primary WD. 

In the optical, spectral features from the WD are difficult to detect due to the 
strong emission from the accretion disc. However, there are systems in which ab- 
sorption features from the WD can be detected also in the optical. These stars typic- 
ally have low accretion rates and faint discs. Figure 123] shows the optical spectrum 
of the faint system SDSS J1507+52 which shows double-peaked disc emission lines 
superimposed on absorption from the WD (Chapter|5]). 



- Donor star: The donors in CVs are typically cool, near main-sequence stars of 
spectral types between G - M, and temperatures between » 3000 K - 5000 K (|Knigge 



20061: 



Knigge et al.ll201lh . Due to the process of accretion, their radii are larger com- 



pared to main-sequence stars of same mass. Also, their spectral types are found to 
be later than they would have been compared to what a Roche-lobe filling main- 



d have. This is because a larger donor 



sequence star at the same orbital period wou 
radius results in a lower density, p oc Mj/Rj ( Kniggell2006 

Infrared (IR) observations are the most promising way to uncover features ori- 



ginating from the cool donor (e.g. 



Hamilton et al 



2011 ). However, in some CVs, the 



donors might be detectable also in the optical, for instance, if the hot WD (and/ or 
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the inner part of the accretion disc) illuminates the front side of the donor. Such 
illumination could give rise to a reflection effect, caused by the increase in emission 
from the irradiation-heated front hemisphere. For instance, this has been detected 
in systems undergoing nova eruptions. 

- Bright spot: The impact point between the accretion stream and disc is often de- 
tected as an variation on the orbital period in the optical light curve. The bright 
spot is less prominent in the IR, and has not been detected in the UV at all. 



- Boundary layer: X-ray observations are sometimes used when studying the pro- 
cess of accretion, and especially for dwarf novae, when studying the extremely 
hot boundary layer close to the WD. The boundary layer is responsible for a large 
portion of the total emission from the system, and therefore, studying this region 
can be considered as highly important. In systems with high M, the boundary 
layer is optically thick and will emit in the soft X-ray range, while for systems with 
lower M, the boundary layer is optically thin and is observed typ ically as thermal 



bremsstrahlung radiation in the hard X-rays (jKuulkers et al. 



2003 ) 



2.2 Analysing Data 
2.2.1 Time-Series Analysis 

A common goal of astronomical time-series analysis is the detection and charac- 
terisation of periodic signals in noisy data. Such work is usually carried out by 
statistical analyses of the frequency content of the data. In this following section, I 
will provide a brief introduction to this type of analysis. 



The lowest frequency that can be detected in a light curve is determined by the 
length of the data set, while the highest detectable frequency in a light curve is de- 
termined by the sampling of the data. In frequency space, the Nyquist sampling 
theorem states that for evenly spaced data, the sampling frequency f s, should be 
least twice the highest frequency expected in the data, the so-called Nyquist fre- 
quency Fn (Fn - fs/2)- A frequency greater than will not appear at the correct 
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frequency. Figure |Z5l shows an example of this problem, where the sampling rate, 
marked as black dots, creates a false periodic signal (blue) of the original signal 
in the data (red). This problem is called aliasing, and refers to when periodic sig- 
nals in the source are not reconstructed at the correct frequencies in the time-series 
analysis. Any data set being sampled will suffer from aliasing. In particular, multi- 
night observations will, apart from the aliasing associated with the sampling rate 
of the observations during the night, also show frequencies introduced by the daily 



sampling. Figure 3.2(a) shows an example of aliasing present in the power spec- 
trum of SDSS J1457+51. Also, as a consequence of multi-night observations, the 
data are no longer evenly spaced in time, which means that the Nyquist frequency 
cannot be specified by a single value. 

To minimise the problem with aliasing, one needs to obtain data sets much 
longer than the period /periods expected in the data, and if the data are obtained 



over many nights, as densely and u niformly spaced in tirne as 



possible. The Cen 



ter of Backyard Astrophysics (CBA: 



Skillman & Patterson 



1993; 



Patterson 



1998ah . 



is a collaboration between small telescopes evenly spaced both in longitude and 
latitude. This is an efficient way of observing that results in minimal problems 
with aliasing, since it results in a light curve with fairly regular sampling, free from 
the problems introduced by daily observations from the same site. Another, more 
expensive way to reduce aliasing is to use space telescopes that are able to con- 
tinuously study a source without interruptions. Thi s is done for instance in the 



Cannizzo et al 



2010). How- 



Kepler mission (see the observations of V344 Lyrae by |i 
ever, often, the observing circumstances are such that aliasing cannot be completely 
avoided. In order to manage the problem with aliasing, one could, for instance, 
perform Monte Carlo simulations or bootstrap simulations on the data. When 
performing Monte Carlo simulations, for each simulation, the flux values are ran- 
domly re-distributed within their errors (see Section l2.3.8l l. In the case of bootstrap 
simulations, for each simulation, a random subset of the data is selected, altering 
the sampling pattern and hence also the aliasing pattern. In theory, for a large 
number of such simulations, the real signal would always be found at the same 
frequency, while those caused by aliasing would not. As a result, the real signal 
should statistically be best represented. In practice, one has to be aware that when 
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Figure 2.5: Model showing how the sampling rate (black dots) can produce an alias (blue) 
of the real signal of the data (red). 

Figure from }fittp://sv.wikipedia.org/wiki/Fil:AliasingSines.svg 



throwing away data points, one might also throw away points that signify the real 
period. 

The Fourier theorem states that any periodic signal in the light curve can be 
reproduced by summing a set of sine waves. As a result, the periodicities in a 
light curve can be found by fitting sine waves of different frequencies to the data, 
where the amplitude of each sine wave indicates how much signal is present at 
each frequency. The power of each frequency can be calculated as the square of the 
amplitude. In a power spectrum, the power is presented as a function of frequency. 

All data presented in this thesis are light curves unevenly distributed in time. 
The most commonly used met hod to find period signals i n such data is to construct 



a Lomb-Scargle periodogram (jLombl 



19761: 



Scargle 



1982 ). However, before such a 



periodogram can be constructed, the mean value needs to be subtracted from the 
data, since the sine waves being fitted have a mean of zero. In a Lomb-Scargle 
periodogram, the power at the frequency w of a light curve y(t„), is defined as 



(L^=i y{tn)cos{a){t„ - r))f (eJ)Li y(tn)sm(a){t„ - t))) 



2\ 



Lti cosHcv(t„ - T)) Zti sin^cvitj - T)) 

where a is total variance of the data, and t is 



(2.1) 



Ln=i sin2a)t„ 



tan(2a;T) = 



LJ^^I COSl(jOt„ 



(2.2) 
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In the case where a light curve contains periodic but non-sinusoidal signals, addi- 
tional higher frequency components are required to reconstruct the data. As a res- 
ult, harmonics of the fundamental frequency are detected in the power spectrum 
(e.g. Figure 3.2(a)| |. 



The accuracy of the location of a peak in the power spectrum is dependent on 
the length of the data set, and longer observations results in narrower peaks. While 
the obvious reason for constructing a power spectrum is to identify frequencies 
with power excess, the significance of a peak needs to be evaluated. Visually, this 
is often done by considering signals that have significantly higher power than the 
highest noise peaks in the power spectrum. However, a signal at low power can still 
be of interest, for instance, if it is detected in many data sets. A more quantitative 
approach to find the significance of a peak is to perform a randomisation test, which 
redistributes the data points in time while maintaining the original sampling pat- 
tern. In doing so, any periodic signal present in the data is deliberately destroyed, 
which allow us to find the probability of detecting excess power by random chance. 
Confidence limits can then be set on the chance of detecting excess power at a cer- 
tain frequency. However, this method is not valid if the data are dominated by red 
noise, i.e. if the data points are correlated. 

2.2.2 Trailed Spectrogram 

In Chapter m I present time-resolved optical spectroscopy of T Pyxidis. Such data 
can be effectively visualised by plotting the line profiles as a function of time. A 
trailed spectrogram is a 2-dimensional plot where the region around a chosen line, 
expressed in velocity units (V), is plotted as a function of orbital phase {(p), F{V, (p). 
The orbital motion of the gas in the binary will, due to Doppler shifts, generate 
spectral lines that follow the shape of a S-curve. The trailed spectrogram is then the 
superposition of all such S-curves. 

A trailed spectrogram can be very useful, especially in the early stages of the 
analysis, since it can give a rough idea of the orbital period and helps visualise 
the overall orbital behaviour of the lines. Also, if studying disc emission lines, 
the wings of the line profiles are expected to roughly track the motion of the WD, 
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Figure 2.6: Trailed spectrograms of the Hell line at 4686 A (left) and the Bozven blend at 
4640 A - 4650 A (right) in T Pyxidis. Data are phase binned and plotted over 2 orbital 
periods. 

and can therefore provide an estimate of the radial velocity of the WD (Ki). Also, 
the mean offset in velocity units between the measured line centre and the rest 
wavelength of the line immediately yields the systemic velocity (y). It is also useful 
to compare trailed spectrograms for different lines to each other. If there are phase 
shifts and /or amplitude differences between them, they most cert ainly originate 



frorn different regions within the CV system (for example. Figure 2 in lSteeghs & Casares 



2OO2I shows phase shifts between lines from the WD and donor). 

Figure |Z6l shows trailed spectrograms constructed from the double-peaked HII 
line at 4686 A (left) and the Bowen blend (consisting of NIII and CIII lines) at 4640 A ■ 
4650 A (right) in the spectrum of T Pyxidis. The data were binned in phase and ve- 
locity, and the systemic velocity was removed. Here, the trails are shown repeated 
over 2 orbital cycles for clarity. The colour scale represents normalised fluxes. The 
relative amplitudes of the two peaks in the Hell line vary over the orbital cycle, 
with first one peak, and then the other dominating. This could be due to a non- 
Keplerian velocity distribution of the emission (e.g. a slightly elliptical disc). The 
trailed spectrograms for both the Hell line and Bowen blend show spectral fea- 
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tures, moving together in phase, which indicate that they most likely are formed in 
the same region. 



2.2.3 DoppLER Tomography 



' Sopp ier tomography is an indirect imaging method developed by 



Marsh & Horne 



by 



( 19881) to isolate the different components of a CV iji velo city space (see reviews 



Marsh 



2001 



2005 : 



Morales-Rueda 



20041: 



Steeghs 



20041) . This method is most 



commonly used t o gain irisight into the structure of accretion discs and other ac- 
cretion flows (e.g. ISteeghsll2003l) and is analogous to the method used in medical 
X-ray tomography. Doppler tomography uses the fact that the spectral lines follow 
S-shaped radial velocity curves as a function of orbital phase as described in the 
previous section (e.g. Figure |231) . A Doppler tomogram therefore presents the flux 
distribution from a spectral line in velocity space. The velocity centre in the map 
corresponds to the centre of mass, which will be close to the WD. Doppler tomo- 
gram are often rotated so that (/) = is at the superior conjunction of the WD (where 
the donor is closest to us). In this orientation, the donor star appears at (p = 0.5, 
where orbital phases are defined to increase in the clockwise direction. However, 
before the data can be put on a 2-dimensional velocity map, the wavelengths and 
orbital phases needs to be converted into the velocity components Y^f A po- 
sition within the binary has an associated velocity vector in the co-rotating frame, 
reflecting the motion of the gas. The observed velocity V(0) is then the radial velo- 
city component of that gas vector at a given binary orientation, also correcting for 
any systemic component (7), 



V{(p) ^y-V^ cos{2n(p) + Vy sin{2n(p). 



(2.3) 



There are several methods for converting the spectral components into velocity 
space, i.e. a Doppler tomogram. The simplest method compares predicted S-curves 
for every possible velocity in the map with the data in the trailed spectrogram, by 
assigning the integrated flux D along each S-curve in the trailed spect rograms to 



the corresponding V^c, point in the Doppler map (M). 
this in the following way 



Horne 



( I992I ) presented 
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Figure 2.7: Doppler tomogram ofT Pyx constructed from the He II line at 4686 A. The 
different components in the system are marked. 



^D(V(V^,Vy,(P),4>)cv{cp)d,p 
J a){(p)d(p 



(2.4) 



where w is a we ighting function. This meth od is also referred to as back-projection 



(as described by 



Marsh 2001 



Steeghsl2004h . In this method each spectrum is projec- 



ted onto the map, rotated with an angle corresponding to the phase of the observa- 
tion. An extension to the method described above is to apply an iterative approach 
to const ruct ar\ optim al map, such as the maximum entropy regularisation, first 
used by lHornd ([19851) when performing eclipse mapping. 

Figure IZTl presents a Doppler map constructed from the trail of the Hell line seen 
to the left in Figure 12.61 The most prominent features, such as the positions of 
the WD, donor star, bright spot and accretion disc/ stream are marked. Since the 
outer part of the disc has lower velocity than the inner parts, the disc is turned 
inside out in velocity space. The velocity semi-amplitude of the WD (Ki) and the 
donor (^2) could, in principle, be found directly from the map (however, to find 
K2, one obviously would need to study lines which actually show emission from 
the donor). Also, the velocity at the outer disc radius, VRj._,^sini (marked as a black 
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circle) can be measured directly from the Doppler map (see Section [2.3.3|) . 

Doppler tomography is not restricted to the mapping of accretion discs, and 
has been used, for i nstance, for studying the accretion flow in magnetic CVs (e.g. 



Schwope et al.lll997n . It is possible to convert the velocity map into a spatial map. 



but this is a poorly constrained problem. In order to make this conversion, one must 
make assumptions about the distribution of the velocities of the material in the 
system, for example, it is often assumed that material in the accretion disc moves 
in a Keplerian fashion. Since important information might be lost in the conversion 
process, we might as well study the Doppler tomograms as they are, in velocity 
space. 



2.3 Estimating System Parameters 
2.3.1 Orbital Period 

In the context of CV evolution, the orbital period (P) is the most important system 
parameter for a CV since it indicates the evolutionary state of the system. There are 
several methods that can be used to obtain the orbital period, the most reliable one 
being eclipse timing, though this is obviously only available in systems viewed at 
sufficiently high inclination (> 75°). 

In non-eclipsing systems, orbital modulations from the bright spot are normally 
detected, and the orbital period can be found from extracting the frequency content 
of any photometric light curve obtained for a longer time than P (see Section l2.2.Hl . 
Any signal of orbital nature appearing in a power spectrum should be stable in 
frequency, and therefore appear as a narrow peak. For systems with low accretion 
rates, the first harmonic of the orbital frequency is commonly found to have higher 
power than the orbital frequency itself: these kind of double-humped waves are, 
for instance, seen in the power spectrum for both of CVs presented in Chapter |3l 

A spectroscopically determined orbital period is generally more reliable than 
a photometrically determined one. This is due to the fact that one can never be 
sure that the period detected in the light curve is of orbital nature, or that the cor- 
rect alias has been identified in the resulting power spectrum. In order to find P 
spectroscopically, radial velocity curves can be constructed from the time-resolved 
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spectroscopic data. Such a study often involves fitting Gaussian profiles to the 
Doppler-broadened disc emission lines. The line centre of each spectrum is then 
retrieved as a function of time, yielding a sinusoidal radial velocity curve with a 
period equal to the orbital period for the system. Radial-velocity curves can also 
be constructed by using more advanced techniques, such as the double-Gaussian 
method described below in Section l2.3.2l which specifically aims to track the inner- 
most disc regions, and hence the orbital motion of the WD. 

When a radial-velocity curve has been constructed, P is found either by fitting 
a sinusoid to the curve, which has a period correspond to the orbital period, or by 
constructing a power spectrum of the radial-velocity data. 

Spectroscopically determined orbital periods also suffer from the problem with 
aliasing in the resulting power spectrum. In Chapter HI I use bootstrap simulations 
to find the most likely orbital frequency in T Pyxidis (described in Sections 12.2.11 
andlHH). 

2.3.2 Systemic Velocity and Velocity Semi- Amplitude 

The systemic velocity (y) of a system is the velocity at the centre of mass of the 
binary. The velocity semi-amplitudes of the WD (Xj) and the donor (i^i), are pro- 
portional to the masses of the two components {Ki/Kj = M2/M1 - q), and can be 
measured as the amplitudes of their individual radial-velocity curves. 

In Chapter HI I estimate y and Ki from the radial-velocity curves obtained 



from disc emission 



ines in T Pyxidis, by using the double-Gaussian technique 



([Schneider & Younglll980h . where the velocities are extracted from the wings of the 
lines. These are thought to be formed in the innermost regions of the accretion disc, 
closest to the WD (see Figure IZ3l l. As a result, one might hope that they will track 
the orbital motion of the WD. In this method, a chosen double-peaked emission 
line profile is convolved with two Gaussian profiles, separated with a distance a. 
The weighted flux average in the red part of the line, with the weighting function 
being the Gaussian profile, is compared to that of the blue part of the line, with 
the aim to find for which value of the radial velocity, the two weighted averages 
are equal. Assuming that the double-peaked line profile is symmetric, this radial 
velocity corresponds to the Doppler shift. This is done for each of the spectra in 
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Figure 2.8: Illustration of the double-Gaussian technique, showing two different separa- 
tions (a) between the two Gaussians (here with a fixed FWHMj. The red line shows the He 
II line in T Pyxidis. 



the time series, for a wide range of central positions of the Gaussian profiles, and 
values of the FWHM of the two Gaussians. In Figure 12.81 I illustrate the double- 
Gaussian technique performed on the He II line in T Pyxidis, where the separation 
a is changed for a fixed FWHM. 

Resulting radial-velocity curves are created for each set of a and FWHM, and a 
sinusoid of the form V{(p,a) = y{a)-Ki{a)sin{(j)), where V is the velocity at phase (p, 
y is the mean velocity, and Ki is the amplitude, is fitted (for instance, by nonlinear 
least-squares fitting) to the data. In order to find what combination of a and FWHM 
that corresponds to the best output radial-velocity curve, a diagnostic diagram can 
be constructed. In this diagram, y, Ki and (p are plotted as a function of a for 
different values of the FWHM. Such a diagram, constructed from radial velocity 
curves of the Hell line at 4685 A in T Pyxidis, is seen in Figure l47l The best set of a 
and FWHM are found from inspecting for which a and FWHM each of the variables 
(p, y &c Ky are stable. The corresponding best values of a and FWHM are then used 
to construct a final radial velocity curve from which final values of y, K and (p are 
extracted. 

There are additional methods to find the K velocities of the components in a CV, 
for instance, Ki can be found by locating the centre of symmetry of the disc emis- 
sion i n a Doppler tomogram. The centre should be located at V^. = and Vy - -Ki 
(e.g. ISteeghs & Casaresll2002n . First, a Doppler tomogram is constructed from a 
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Doppler map 
(constructed from trail) 





Figure 2.9: Top left figure shows a trailed spectrogram for the Hell line at 4686 A in T 
Pyxidis. Its corresponding Doppler tomogram is seen in the top right figure. Bottom left 
figure shows the symmetry map. In the bottom right figure, the asymmetry map constructed 
from subtracting the symmetry map from the original map is plotted. 

trailed spectrogram. After masking out any regions of strong asymmetric non-disc 
emission, the centre of symmetry is found by taking a trial point, smoothing the 
map azimuthally about this point, and then subtracting this smoothed map from 
the real map. If the trial point is far from the centre of symmetry, the resulting differ- 
ence map will contain large asymmetric residuals. The optimal centre of symmetry 
is thus the trial point for which the difference map exhibits the smallest residuals 
(as measured by the standard deviation of the difference map). 

Figure |Z9l illustrates how Ki can be obtained using the method described above 
for the He II line at 4686 A in T Pyxidis. The upper left frame shows the trailed 
spectrogram used to construct the Doppler tomogram presented in the upper right 
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frame. The corresponding symmetry map centred on V^, Vy, is presented in the 
lower left figure. The lower right figure show the asymmetry map, constructed by 
subtracting the symmetry map from the original map. 

In the Doppler tomogram, the donor star is positioned and Vy - +K2. 

However, in order to locate the position of K2 in maps which are not constructed 
from donor lines, one needs to know JCi and q {K2 - Ki/q). 

2.3.3 Velocity at the Outer Disc Radius 

The projected velocity at the outer disc radius (vRj.^^sinz) can be found by measur- 
ing half of the peak-to-peak (AVpgaj,.to.pga]^) separation in the double-peaked emis- 
sion lines, and/ or by measuring the radius of the inner part of the disc structure in 



the Dopp 



ISmak 



er tom ogram (see Figure IZ7)l . 



([19811) investigated the exact relationship between the true VR^._,^sinz and 
the measured AVpgaj^.to-peak/ taking the effects of instrumental broadening and dif- 
ferent disc emissivity distributions into account. He found that, for a wide range of 
realistic situations, VRj._,^sinz - 0.5 AFpgak-to-peak/"/ where m = 1.05 + 0.05. 

2.3.4 Mass-Ratio 



Warned (|1973l) showed that the radial velocity of the WD {Ki ) divided by a char- 
acteristic disc radius (vR^sinf), is a function of only the mass-ratio (q). Such a 
relation implies that q can be estimated from time-series spectroscopy, where Ki 
and vr^ sin i can be me asured as des cribed in Sections 12.3.21 and 12.3.31 In order to 



derive a relation for q, IWarnen (|1973l ') assumed that a specific particle leaving the 
Lagrangian point L^, will conserve its angular momentum in the orbit around the 
primary WD. Thus, the circularisation radius is used as a characteristic disc radius. 
From Kepler's third law for circular orbits, a particle in an orbit around the WD has 
a velocity v 

2nai {l + q)^ 

^ = uf2i \ 3 — ' (^-^^ 

Pf^{q) 

where fli is the distance between the centre of the primary and the centre of mass 
for the binary, and f{q) is the distance between the centre of the primary and the 
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Li, measure d in units of a l. Combining Equations 12.51 with the expression Ki 



Inui s'mi/P, 



Warnerl (|l973[) found the following relation 



sm I 



{1 + q) 



2 ■ 



(2.6) 



In ChapterlH I estimate q for T Pyxidis in a similar manner. However, instead of us- 
ing a characteristic disc radius (vj^ sin i), I use the velocity at the outer disc radius 
(vR^j^_.sinf), which can easily be measured as described in Section l2.3.3l assuming a 
circular accretion disc. Such an assumption is plausible for any lat e-type binary sys 



tern since the tidal circularisation time scale must be very short (|Zahn & Bouchet , 



19891 ). Also, in the case of T Pyxidis, I assume that the disc radius extends all way 



out to the tidal radius. 



Harrop-Allin & Warneii (|l996l) performed empirical meas 



urements of the radii of accretion discs for eclipsing systems, and found that sys- 
tems with high M have larger discs than systems with lower M. Since T Pyxidis is a 
high M system, a tidal limitation can be expected to be a reasonable ap proximation. 



thus, its disc radius is given by the following equation (|Warner 



1995a|), 



0.60 
(1 + ^)' 



(2.7) 



These assumptions, together with Kepler's third law, yield the expected relation- 
ship (please see Appendix A for a full derivation of this expression). 



K 



- = 0.77 



1+q 



(2.8) 



In Figure 12.101 1 compare the ratio K^/vTj^ sini as a function of q, as derived from 
the expression in Equation 12.61 from IWarneii (|l973l ) (although this is not specifically 
the o uter disc radius), with the expression I present in Equation 12.81 jUthas et al 



2010h . The function / in Equation 12.61 is highly dependent of q, and was here de- 



Plavec & Kratochvil 



rived from tabular values of Roche-lobe models presented by 
( 19641) . High values of q result in a smaller value of the quote Kj/v^^. sint for 
Equation 12.61 compared to Equation 12.81 For example, the system T Pyxidis, has 
a measured value of Ki/v^^^^sini - 0.130, resulting in q = 0.20 if obtained from 
Equation |2.8[ while the same quote would result in a higher q (« 0.32), if calculated 
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Figure 2.10: The relation between q and the quote Xj/vr^^, sinz is presented as derived in 
Equation \2.6\ (red line) and in Equation \2.8\ (black line). 



from Equation 12.61 Since IWarneii (|1973l) used the fact that the disc radius cannot 
be smaller than the circularisation radius, this means that Equation 12.61 effectively 
yields an upper limit on q. 

2.3.5 Inclination 

The inclination (i) of a CV can be accurately determined for eclipsing systems, 
where the duration of the eclipse is a function of i, q and P. In order for an ec- 
lipse to occur, i has to be ^ 75°. For non-eclipsing systems, i £ 75°, and can be 
estimated if P, and q are known, from Ki = vi sini, where Vi is defined as 

27tfl q 

^'-—uiiry 

where the binary separation a can be obtained from Kepler's third law; see Equa- 
tion IA.3I in Appendix A. The system inclination can also be estimated using the 
mass functions, assuming that P, Ki and either the WD or donor mass is known 
(Equations |A.5| and |A.6| in Appendix A). 
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2.3.6 Temperature and Metallicity of the White Dwarf 

As noted above, in most CVs, the primary WD is difficult to observe since it is 
often obscured by light from the accretion disc and bright spot. Nevertheless, in 
some systems with orbital periods close to the minimum orbital period, the ac- 
cretion discs are faint and the WD contributes significantly to the overall light 
from the system, as evidenced by the presence of its absorption lines in the spec- 
trum. Magnetic CVs also provide a good opportunity to study the primary WD 
since they do not have pronounced accretion discs. Also, in some CVs, the ac- 
cretion temporarily switches off from time to time, providing a good opportunity 
for studying the WD without interference from the accretion disc. Such low states 
have been observed, for instance , in VY Scl stars - a subtype of nova-like variables 



(e.g. iTownsley & Gansickell2009l) . Since WDs are hot, they are bright in the UV. 
Therefore, UV observations are commonly used to study WDs in CVs. By com- 
paring the UV data spectrum with synthetic model-spectra, a number of WD para- 
meters can be established, such as the effective temperature (Tgff), surface gravity 
(logg), metallicity {[Fe/H]) and rotational velocity (vsinf). 

The Tgff and logg affect the gradient and overall shape of the spectrum and are 
highly correlated with each other. The two top panels in Figure l2Tll show synthetic 
model-spectra produced by Synspec at three different Teff (14000 K [green], 16000 
K [red] and 19000 K [blue]) and at three different values of logg (7.5 [green], 8.0 
[red] and 8.5 [blue]). All other parameters were kept at fixed values. As seen from 
the figures, higher effective temperatures and lower surface gravity produce bluer 
spectra. The two top panels in Figure 12.111 also demonstrate that, at these sorts of 
temperatures, it is much easier to determine Tgff and logg in the blue end of the 
spectrum compared to the red end (i.e. the differences between two Tgff values, or 
two log^ values, are much larger in the blue end of the spectrum). Also, there is 
a higher possibility of contamination from the accretion disc and bright spot at the 
red end of the spectrum. As a result, the far UV-region may be expected to provide 
more reliable estimates of Tgff and \ogg. 

The parameter [Fe/H] determines the depth of the spectral lines while vsini 
affects both the depth and shape of the lines. As a consequence, these two para- 
meters are also correlated. The bottom left panel in Figure 12.111 shows models at 
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Figure 2.11: In the top figure, synthetic model-spectra for T^j^, 14000 K (green), 16000 K 
(red) and 19000 K (blue), is shown (other parameters are fixed at log(g) = 8.0, vsini = 
200 kms~^ and [Fe/H] = -1). The middle panel show three different values oflog(g), 7.5 
(green), 8.0 (red) and 8.5 (blue). Here, T^ = 15000, vsinz = 200 kms'^ and [Fe/H] = 
-1. The bottom left panel show [Fe/H] at, (green), -1 (red) and -2 (blue) (Tg^ = 15000, 
log{g) = 8.0 and vsini = 200 kms~^). The bottom right panel present vsini at, kms~^ 
(green), 200 kms-^ (red) and 500 kms''^ (blue) (Tgff= 15000, log(g) = 8.0 and [Fe/H] = 
-1). Models are produced by Synspec. 

three different values of [Fe/H] ([green], -1 [red] and -2 [blue]). The bottom right 
panel shows models at three different values of vsinz (0 kms~^ [green], 200 kms~^ 
[red] and 500 kms'^ [blue]). As seen in the two figures, the lines become shallow 
for lower values of [Fe/H], and wide and shallow for higher values of v sin i. When 
making a model fit to a data spectrum, even the best fit might not perfectly agree 
with the data both in overall gradient, width and depth of the lines. Therefore, 
it might be better to focus on fitting to narrow line regions to be able to establish 
reasonable values for [Fe/H] and vsini. Also, to make such a fit reliable, both Tgff 
and logg should be kept kept at fixed values, determined from a fit to the overall 
spectrum (as described above). 

The models presented in Figure IZTT] were used to fit the UV spectrum of SDSS 
J1507+52, as described in Chapter |5l Also, the correlation between the parameters 
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is explored in Chapter |5]by constructing contour plots (see Figure l5!4l l. 



2.3.7 Distances 



From Trigonometric Parallax 



There are many ways of estimating the distance to a source, the best known being 
trigonometric parallax measurements. In this method, a star's annual movement 
over the sky as the earth moves around the sun, is compared to more distant stars 
in the background that appear not to move. The distance (d) to the star in parsecs is 
given hy d = l/p, where the parallax (p) is the angle by which the star has moved. 
However, this method is only possible to use for stars that are not too far away, and 
the limitation if using a telescope from earth is an angular resolution of ~ 0.01 arc 
seconds, corresponding to a distance of ~ 100 pc. In order to measure a smaller 
parallax, space telescopes are needed. For instance, the Hipparcos satellite (an ESA 
mission in operation during the years 1989 -1993) could measure a parallax as small 
as 0.0014 arc seconds, and find di stances for stars out t o ~ 700 pc (an overview of 



the Hipparcos Mission is given by 



Ferryman et al. 



1997h . 



Due to the fact that CVs are faint objects, historically it has been a rather diffi 
cult task to find distances towards CVs from trigonometric parallax m easurements. 



Nevertheless, such studies are nowadays quite common. For instance. 



Thorstensen 



|2003l) developed a method that uses measured parallaxes together wi th proper mo 



tions and absolute magnitude s to find distances towards CVs (see also 



Harrison et al. 



2004 



Thorstensen et al 



2008h . 



From Donors 

The distance towards CVs can also be estimated by assuming that the donor star 
has the same surface brightness that an isolated main-sequence star with the same 
spectral type would have (for which the distance is known, for instance, obtained 
from parallax measurements). The flux in a CV is often dominated by the accre- 
tion disc, and therefore the surface brightness of the donor star is often measured 
from obtained observations in the IR wavelength region, where the donor is at its 
brightest. 



52 



2.3. ESTIMATING SYSTEM PARAMETERS 



Bailey! ([19811) developed a method to find distances to CVs from IR observa- 
tions, by measuring the mass, spectral type or temperature of the donor, and com- 
paring this with main-sequence stars to find absolute magnitudes and hence also 



distanc e s. This method has been further dev eloped and presented by 



|2005aD; 



BCniggd ||200d); 



BCnigge et al 



Patterson et al 



( 20 111 ). For instance, observations in the K- 



band yield 



logd 



^ K-Mk,2{P) + 5 



(2.10) 



where K is the apparent magnitude and Mj^ 2 is the absolute magnitude, which can 
be estimated for norm al, unevolved CVs just from the orbital period (see Figure 17 



m 



Knigge et al 



201 ih . 



From Dwarf-Nova Outbursts 

During every dwarf nova eruption, the accretion disc is expected to become optic- 
ally thick with a surface temperature of about T = 20000 K. Therefore, t hese erup- 



tions can act as standard candles and be used for distance measurements. 



Patterson 



( 20111) presented an empirical study of the relation between the absolute mag- 
nitudes and orbital periods during dwarf-nova outbursts for 46 short-peri od CVs 



This study was an irn provement on the previous relation presented by 



Warner 



(h98Z). 



PattersonI (|201ll) finds that for erupting dwarf novae, (My)„^ax = 5.70 - 0.287 



P(hr), resulting in a final error on the correspon ding distance of about 10% - 30%. In 
particular, for superoutburst in SU UMa stars. iPattersonI (|201ll) finds that {Mv)max 
= 4.95 - 0.199P(hr). SU UMa stars enter long-lasting high-states between the su- 
peroutbursts, and therefore these stars are particularly good for distance estim- 
ates, where the distance in units of parsec is logrf = (V^piat - -Ay - 0.5)/5, where 
Ay = 3.1£(B- V) < 0.2 mag. Using this method, the errors on the distances es- 
timates are « 15%. 



From WDs 

If the light from an accreting WD can be isolated (from that of the accretion disc 
and bright spot), we can fit stellar-atmosphere models to the WD spectrum (as de- 
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scribed in Section [2.3.6l l to estimate distances by comparing data and model fluxes. 
The relation between observed fluxes (F) and model fluxes (here Eddington fluxes 
[H]) are found from 



where J?i is the radius of the WD. Therefore, the distance is a function only of Ri 
and the normalisation factor (N) between Eddington fluxes and observed fluxes 



The radius of the WD is not known from the stellar atmosphere models, but the fits 
do provide an estimate of the surface gravity, logg. Since WDs follow a mass-radius 
relation that depends only weakly on temperature, a fixed value of log g also fixes 
Ml and Ri, via 



where Mi is the mass of the WD. In Chapter |5l I use this method to estimate the 
distance towards SDSS J1507, as well as the mass and radius of the accreting WD 
in this system. However, these estimates of d, Ri and My, rely on the fact that the 
atmospheric parameters found from the model fitting are determined accurately, 
for instance, that the spectrum is purely WD with no contamination from disc or 
bright spot. For eclipsing systems, the mass and radius of the WD may be available 
directly from eclipse modelling, yielding a more accurate estimate of the distance. 



F = 



4=nRjH 



(2.11) 




(2.12) 
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2.3.8 Finding Errors on Parameters 

Errors on Individual Measurements and their Mean Value 

The data points for any light curve, spectrum or radial velocity curve have formal 
errors calculated by the programs used to produce them (I have used IraJ^ for pro- 
ducing light curves and MOLL^Ofor calculating radial velocity curves). Formal stat- 
istical errors on each data point can be estimated by calculating the statistical prop- 
erties of multiple measurements, where the standard deviation can be expressed 
as 



a - 



^ I -^2 

1 ! = 1 



where x is the mean value and N is the number of points. Generally, formal para- 
meter errors can be calculated on the assumption that measurement uncertainties 
follow the usual normal distribution (Gaussian distribution). Note that at least one 
of the methods used to find parameter errors in this thesis (bootstrap simulations), 
do not make this assumption since it does not take the errors on individual data 
points into account (see Section l2.3.8D . The normal distribution shows the probab- 
ility of finding a value within a certain o. The \-o error range from the mean value 
corresponds to about 68% of the area of the normal distribution. About 99.9% of all 
points are expected to be found within the 3-a error range. The normal distribution 
is defined by the probability density function /(x) 

f{x) = —=e^. (2.15) 

The error on the mean value derived from multiple number of measurements is 
a/yfN. 



IRAF is a astronomical data analysis software for photometric and spectroscopic data, developed 



by the National Optical Astronomy Observatories (NO AO) Ihttp:/ /iraf.noao.edu/ 1. 
^MOLLY is a program by Tom Marsh for analysis of time-series spectroscopy 

(http: / / deneb.astro.warv^ick.ac.iik/phsaap/software/ molly /html/INDEX.html). 
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X^-fitting to data 



When comparing a model to a dataset, the goodness-of-fit is used to describe how 
well the model matches the data, and is usually evaluated on the basis of the 
statistic. The X'^ is defined as 



The 1-cr error on any parameter can be estimated by identifying the parameter 
ranges that correspond to x'^ = X^in- 1- The mean of a x^ distribution is its number 
of degrees of freedom (v). The reduced X'^ ixl) is defined as the X'^ divided by v, 
and can be expressed as 



If the value of xl < h the model is likely over fitting the data, and the initial errors 
are overestimated. A value of > 1 indicates a poor fit to the data, and it is likely 
that the data errors are underestimated, or that the model is inadequate to describe 
the data. Throughout this thesis, in those cases where xl gave a poor fit, I modified 
the data errors so that xl = ^ (thus making the errors more realistic). Such an 
adjustment of errors is conservative, since it yields larger errors on the parameters. 
This was done in two different ways, either by adding an intrinsic dispersion term 



to the errors, so that the total error corresponds to Ofotal = y + (^Inf scaling 
the errors with a multiplicative factor. Flux /magnitude errors on light curves were 
rescaled by comparison with a multiple sine-wave model, while errors on radial 
velocities were modified by comparison against a single sinusoid fitted to the radial 
velocity curves. 

Correlation between Parameters 

There is at least one additional consideration that has to be taken into account when 
finding errors on parameters. If two parameters are likely to be correlated, the cov- 
ariance should be included in the total error. The correlation coefficient p between 
the two parameters A and B can be expressed as 




(2.16) 




(2.17) 
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covjAB) E[(A-xa)(Y-Xb)] 
Pa,b = = / U-io) 

where E is the mean value operator. The value of 04 r varies bet ween + 1, where a 



198% . Instead of calcu- 



value of indicates that A and B are uncorrelated (|Barlowl l 
lating the covariance matrix, 2-dimensional contour plots of the two parameters in 
space can visually reveal any correlations between parameters (e.g. Figure l5!4l) . 

Finding Errors from Data Simulations 

In this thesis, the final errors on almost all system parameters presented in Chapters 
millandlSl are found from Monte Carlo or bootstrap simulations. When performing 
Monte Carlo simulations (first proposed by Enrico Fermi in the 1930s), the meas- 
ured values are randomly redistributed within their 1-a error range. This process 
is iterated many thousands of times, and for each time, the analysis is repeated 
yielding a new result. The distribution of all results from all simulations is then 
used to find the mean value and its corresponding 1-a error, either from fitting 
a Gaussian to the output distribution or directly from statistical calculations (i.e. 
Equation l2.14|l . When performing bootstrap simulations, mock datasets are created 
from points that are drawn at random with replacement from the original dataset. 
In the case where the data are time series, the mock data sets do not share the ex- 
act same sampling pattern as the original data. In bootstrapping, typical about 
2/3 of the original data will appear in each mock data set, and hence, some values 
will appear twice. For the error analysis presented in Chapter|3j I have performed 
both Monte Carlo and bootstrap simulations, to ensure they produce similar res- 
ults. Since the bootstrap method is independent of the errors on the data points, 
this provides a cross-check that the errors are sensible. For example, if the errors 
are under estimated, the Monte Carlo method would yield a lower resultant error 
than the bootstrap method. 



Systematic Errors 

The error ranges for any of the system parameters presented in this thesis were 
carefully defined by using the methods described in this section. Sometimes, many 
different analysis methods were used to find one parameter. For instance, Ki can 
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be foimd both from the sinusoidal radial velocity curve and from the true centre of 
symmetry in the Doppler tomogram. Another example occurs when retrieving res- 
ults of the same parameter from different time spans or wavelength regions. In case 
these different methods and/ or choice of data do not yield the same result within 
their statistical formal errors, the range in which the different results are found 
could be a better measure of how well we are able to determine the parameter. This 
will give a lower bound of the error, resulting in a systematic error rather than a 
formal, purely statistical one. 

2.4 Obtaining Data 

Generally, optical photometric observations of CVs can be done from rather small 
telescopes (< 1 m), while spectroscopy normally requires slightly larger telescopes 
in order to achieve an acceptable signal-to-noise ratio (S/N). When observing CVs, 
the main contributing factor to an increased data quality is time. However, there 
is a trend at the world-leading observatories to close down smaller telescopes in 
favour of larger ones. A larger photon counting area is often an unnecessary ex- 
pense of both time and money. However, larger telescopes are sometimes needed, 
for instance when observing very faint systems or when a high time resolution or 
spectral resolution is required. 

When planning observations, one has to make sure that the data being obtained 
match the science goals. For instance, if the goal is to measure the orbital period 
from phase-resolved spectroscopy, the time resolution has to be small enough com- 
pared to the orbital period not to smear the orbital signal. As a consequence, the 
chosen optical setup, such as the grism and slit, must allow for enough throughput 
to obtain a reasonable level of S/N at this required time resolution. Normally, calib- 
ration frames such as bias frames, dark frames and flat fields, are obtained close in 
time to the observations. For flux calibration, a flux standard star observed close in 
time and at the same airmass at the data observations is required, but can only be 
obtained if the night has photometric conditions. Below follows a brief description 
of the calibration frames required for standard data reduction. 



58 



2.4. OBTAINING DATA 



2.4.1 Dark and Bias Frames 

On a CCD, after each exposure the analog signal containing the niimber of electrons 
in each pixel is converted into a digital signal, often expressed in ADU (Analog 
Digital Units). The conversion from photons to electrons is not perfect, and various 
sorts of noise are present on the resulting frame. Each pixel charge carries a fixed 
offset voltage, a bias. There is also noise caused by thermally generated electrons, 
the dark current. Thus, even if obtaining an exposure of seconds, there is still 
a signal varying from pixel to pixel. To compensate for these effects, an exposure 
of seconds (a bias frame) is taken to determine the bias pattern. A dark frame is 
obtained with the same exposure time as the science frames, but with the shutter 
closed. 

Proper cooling of the CCD is necessary to reduce the amount of dark current. At 
larger telescopes the CCD detectors are normally cooled with liquid nitrogen. As 
a result, the dark current only becomes significant for very long exposure times. If 
the CCD is not cooled at all, the dark current dominates the noise (together with the 
sky background), while the noise contribution from the bias might be negligible. 

2.4.2 Flat Fields 

Each pixel on the CCD has a slightly different sensitivity to the incoming photons. 
Flat fields are images obtained to correct for these pixel-to-pixel variations. A twi- 
light flat uses the twilight sky to imiformly illuminate the CCD, while a dome flat 
uses that from a bright (usually quartz/Halogen) lamp. The pixel sensitivity de- 
pends on wavelength, and therefore a set of flat fields must be obtained for each 
filter that is being used. The exposure time varies with optical setup, but the pixels 
should be filled to about half of their full well capacity, or about 60% of the satura- 
tion level, which normally ensures a good S/N in each flat without causing satur- 
ation in bright areas of the field. Small-scale variations in the frame, such as dust 
appearing in the optics, can sometimes also be removed by flat fielding, as well as 
larger-scale variations, such as vignetting. A vignetting effect is produced when 
the centre of the CCD is brighter than the corners, due to uneven distribution of 
light. Also the use of filters which have a higher transmission in the centre causes 
dark corners on the resulting frame. Other effects, for instance fringing, can also be 
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diminished by a flat-field correction. Fringing is seen as an interference pattern in 
the output image and is thought to be caused by light being reflected as it enters 
the CCD. Photons of higher energy are absorbed by the silicon in the CCD and 
have short absorption depth. However, photons with lower energy have a larger 
absorption depth, and when this depth becomes equal or larger than the thickness 
of the CCD, the photons go through the CCD and reflects back in the (backside illu- 
minated) CCD, creating interference pattern both in photometry and spectroscopy. 
Fringing becomes a problem for longer wavelengths and thin CCDs. 
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Chapter 3 

sdss j1457+51 and bw sculptoris 

- Two NEW ACCRETING AND PULSATING WDS - 



I would like to start this chapter by acknowledge the contribution from my collaborators. In 
■particular, I want to thank Joseph Patterson for providing me with data for BW Sculptoris, 
and for being responsible for the analysis and interpretation of this data. I have formalised 
and integrated his preliminary work with my own analysis of the target SDSS J1457+51. 
Due to the similarities between these two systems, comparison with the BW Sculptoris data 
were invaluable and provided many insights during my own analysis of SDSS J1457+51. 
For completeness, I here include the analysis of BW Sculptoris, even though this is not 
primarily my own work. 



Uthas et al. 



tOllW . with 



The main contents of this chapter will appear in MNRAS as 
co-authors; Joseph Patterson, Jonathan Kemp, Christian Knigge, Berto Monard, Robert 
Rea, Greg Bolt, Jennie McCormick, Grant Christie, Alon Retter and Alex Liu. 



N on-radial pulsations (NRPs) are commonly found in isolated white dwarfs 
of DA type, often called ZZ Ceti stars. These stars have hydrogen-rich atmospheres. 
Pulsations occur as the WD cools and p ass through a phase o f pulsational instabil- 
ity, and are detected mainly as g-modes (jGianninas et al.l2006h . Similar signals have 
also been detected in faint CVs, which are interpreted as NRPs from the primary 
WD. All CVs displaying this type of signal also show absorption from the imderly- 
ing WD in their spectrum (se e Sectionll. 2.311. Here, the se systems are called GW Lib 



stars, after the first discovery (jWarner & van Zyl 



19981). 



In this chapter, I present time-series photometry of two CVs that appear to be 
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certifiable members of the GW Lib class. Both systems have very low accretion lu- 
minosity and show signatures from the white dwarf in their optical spectra. Also, 
both systems show double-humped orbital signals and non-commensurate peri- 
odic signals, indicating non-radial pulsations. One is SDSS J1457+51, which has a 
photometric wave suggesting an orbital period of 77.885+0.007 minutes. The other 
is BW Sculptoris, with P^j-b - 78.22639 + 0.00003 minutes. Both stars show pulsa- 
tions mainly near 10 and 20 minutes. These rapid signals drift slightly in frequency, 
and may consist of several, finely spaced components. BW Sculptoris also shows 
a remarkable photometric variation at 87 minutes, which could be explained as a 
quiescent superhump, possibly arising from a 2:1 orbital resonance in the accretion 
disc. 



3.1 SDSS J1457-F51 



SDSS J1457+51 (hereafter J1457) was first identified in the Sloan Digital Sky Survey 
by Szkody et al.l ( 2005 ). They obtained spectroscopy that showed the broad absorp- 
tion characteristics of a white dwarf, indicating a system of low accretion rate. The 
source was found to be faint (g « 19.5). Due to the double-peaked nature of the 
emission lines, the y sugge s ted th e system to be of high inclination. 



The same year. 



Dillon 



(|2008[) obtained both spectroscopy and photometry of 
J1457, but no orbital period could be established. She found some evidence that the 
system might harbour non-radial pulsations and also identified a periodic signal at 



85 minutes, which she concluded must be close to the orbital period. 



3.1.1 Observations, Data Reduction and Analysis 

Time-resolved photometry of J1457 was obtained with the 1.3 m and 2.4 m MDM 
telescopes at the Kitt Peak observatory, Arizona, during April and May 2010. The 
star was observed during 14 nights in total, spread over 47 days. With all data 
coming from the same terrestrial longitude, we were not immune from aliasing 
problems and therefore strove to obtain the longest possible nightly time series 
(generally ~ 4 - 7 hours). Weather conditions such as clouds, snow and full moon 
prevented us from obtaining times series over more than 4 consecutive nights at 
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a time. The time resolution was in the range of 20 s - 30 s. Table ITT] presents a 
log of the observations, where the column for quality indicate the quality of the 
night (in terms of weather and contribution from moon light) A clear filter with 
a blue cutoff was used to minimise differential-extinction effects and allow for a 
good throughput. 

The data reduction was done in real time during the observations, using stand- 
ard IRAF routines. The data consisted of differential photometry with respect to 
the field star USNO A2.0:1350-08528847. The search for periodic signals was ini- 
tially done for single nights separately, and Lomb-Scargle periodograms (jScargle 



19821) were constructed. Formal flux errors were rescaled so that the xl ~ ^- This 



was done by fitting a fake light curve to the original data, composed of multiple 
sine waves with periods corresponding to the strongest signals found in the single- 
night power spectrum. Monte Carlo simulations were performed on every peak 
of interest in the power spectrum to find the period and its error. In this method, 
the peak errors are found by randomly re-distributing the points in the light curve 
within their errors, a repeated number of times, and constructing a Lomb-Scargle 
periodogram each time. The 1-cr error is then found by fitting a Gaussian to the 
output distribution of the peaks found in the periodograms. Data from several 
nights were then combined to allow the search for signals with lower amplitude, 
and also to improve the frequency resolution. Monte Carlo simulations as well as 
bootstrap analysis was performed to find errors on the peaks, and to distinguish 
between the most likely aliases. The two methods have different advantages and 
were used to provide a cross-check of the results. Monte Carlo simulations use 
the whole dataset and are only truly reliable when the errors are well estimated, 
while bootstrap analysis (formally bootstrap-with-replacement) does not consider 
the errors on the data points, and so can be more reliable in cases where the errors 
are not well estimated. Bootstrap analysis can also help to distinguish between the 
most likely aliases, since the sampling pattern of the data is changed by creating 
a mock dataset where a subset of the original dataset is randomly created and so 
with repeated use this will efficiently destroy or weaken the alias pattern. Results 
from the two methods agreed well with each other. 
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Date 


HJD 


Telescope 


Length 


Quality 




(-2455000) 




(hours) 




100414 


301 


1.3m 


2.30 


medium 


100416 


303 


1.3m 


7.89 


excellent 


100417 


304 


1.3m 


4.32 


bad 


100419 


306 


1.3m 


4.72 


medium 


100424 


311 


1.3m 


5.42 


medium 


100425 


312 


1.3m 


7.82 


bad 


100503 


320 


1.3m 


7.92 


good 


100504 


321 


1.3m 


8.26 


medium 


100506 


323 


1.3m 


8.40 


good 


100507 


324 


1.3m 


3.50 


medium 


100508 


325 


1.3m 


5.47 


medium 


100528 


345 


2.4m 


4.97 


bad 


100529 


346 


2.4m 


8.11 


good 



Table 3.1: Summary observing log for SDSS J1457+51. Data was obtained at the MDM 

observatory during April and May 2010. The column for quality is a measure of how good 
the night was (in terms of weather and contribution from moon light). 
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3.1.2 Light Curve and average Power Spectrum 

The mean power spectrum, averaged over the six best nights, is shown in Fig- 



ure 



3l(a)| The 18.48 cd'^ and 36.97 cd'^ signals are almost certainly the orbital 
frequency (cOg) and its first harmonic (IcOq). This kind of variations at 2cOg are 
commonly seen in the orbital light curves of CVs with low accretion rates (for in- 



stance in WZ Sge; see iPatterson et al.ll2002bi . The signals found in the range 142 



cd"^ - 148 cd^^ are non-commensurate with the orbital frequency. These peaks 
vary slightly in period and amplitude when present at all in the nightly power 
spectrum. The complex structure around them indicates either an unresolved fine 
structure or periods varying from night to night; it is discussed in detail below. 
Also, during a few nights, peaks were found at 135 c d~^ and 72 c d~\ which are also 
non-commensurate with the orbital frequency. The lower region of the power spec- 
trum shows strong peaks at 4 c d"^ - 6 c d"\ corresponding to the typical length of 
a nightly observing run. The unit cycles per day (cd^^), is used throughout this 
chapter since it is the natural unit for the sampling pattern of the multi-day light 
curves. Also, it clearly shows the natural daily alias pattern. 



The average power spectrum in Figure 3.1(a) has low frequency resolution, 
since each night is less than 8 hours long. The power spectrum of a spliced light 
curve spanning several nights (the coherent power spectrum) is in principle better, 
since the resolution is always near 0.1 T"^ c d~^, where T is the duration in days. It 
does, however, make the assumption that a candidate periodic signal is constant in 
period, phase, and amplitude over the duration of the observation. Power spectra 
will be difficult or impossible to interpret correctly when this assumption is grossly 
violated. 



Figure 3.1(b) shows the normalised and smoothed light curve for a sample 
night. A model light curve constructed from the four strongest peaks found in 
the power spectrum that night (including both co^ and 2cOg), is plotted on top of the 
smoothed light curve. Peaks found at frequencies higher than 40 c d"^ are not rep- 
resented in the model light curve. During one of the observing nights, we obtained 
multicolour data and found the brightness of of the star to be V = 19.2 + 0.2. The 
mean brightness on each night was constant within the measurement error of 0.03 
mag. There was essentially no evidence of flickering, a nearly universal and ob- 
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vious feature of CV light-curves. This, together with the fact that absorption lir\es 
are seen in the optical spectrum implies that the total light in the range 4000 A - 
7000 A is dominated by the white dwarf. 

3.1.3 The Orbital Signal 

A dominant, stable feature at « 37 cd^^ is always present in the nightly power 
spectra. In four of the nights, a weaker, but stable signal is also present at half that 
frequency, ~ 18.5 cd"^. These two signals are interpreted as the orbital frequency 
and its first harmonic, cOq and la)^. As mentioned above, double-humped orbital 
waves are quite common among CVs, especially those of very low luminosity (for 
instance in WZ Sge and AL Com). A power spectrum composed of 11 nights, span- 
ning 45 days, yields = 18.4888 + 0.0017 cd^^ and 2cOo = 36.9740 + 0.0005 cd"^. 
Errors are calculated from bootstrap simulations as described in Section 13.1.1 1 



Figure 3.2(a) shows the low-frequency portion of the full 11 night power spec- 
trum. A zoom of the region around the orbital frequency is plotted on top. A model 
power spectrum constructed from two artificial sinusoids at cOq and Icv^, using the 



exact same sampling as for the original dataset, is shown in Figure 3.2(b) When 
comparing model versus data, the surrounding picket-fence pattern is similar in 
structure and height. This indicates that the orbital signal indeed maintains an es- 



sentially constant amplitude and phase. In Figure 3.2(c) data from all 11 nights are 



folded onto the orbital frequency, showing the double-humped orbital wave at cOg 
and IcVq. 

The lower frequency range of the power spectrum was further investigated to 
rule out the possibility of signals hiding in the noise (see Section D.2.5| for the case of 
BW Sculptoris). The power spectrum was cleaned from the strongest signals at cOq, 
IoJq and also from the high-amplitude peaks between 4 cd"^ - 6 cd"^. However, 
no additional peak was found in this region or in the vicinity of the orbital period. 



3.1.4 High-Frequency Power Excess 

The complex signals spanning 142 c d~^ - 148 c d~^ move slightly in frequency and 
are non-commensurate with the orbital frequency. These signals could be inter- 
preted as non-radial pulsations arising from the primary white dwarf. In addition. 
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18-48 36,97 



20 



10 





100 

Frequency (c/d) 



200 



(a) The mean Lomb-Scargle periodogram for J1457, composed of data from the six nights of 
best quality. The orbital period at 18.48 cd~^ and its first harmonic are plotted as solid lines. 
0.2 




320.65 



320.7 



320.75 



320.8 320.85 
HJD-2455000 



320.95 



(b) The smoothed and normalised light curve of J1547 from one single night. A model 
light-curve constructed from the four strongest periods (including cOg and 2wg) found in 
Tisure \S^l(a)\ is plotted together with the data. 



Figure 3.1: Mean power spectrum and single-night light curve ofJ1457. 
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(a) The power spectrum of ]1A57 from 11 nights, showing the or- 
bital signal cOg and its first harmonic, Icog. A zoom of the region 
around cog is plotted on top. 
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(b) A model power spectrum ofJ1457 constructed from two sine 
waves at Wg and 2wg, using the same sampling pattern as in Fig- 
ure \3.2(a)\ 
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(c) Light curve of J1457 folded on the cOg frequency, showing a 
double-humped orbital wave. 



Figure 3.2: The orbital signal in J1457. 



star 


Frequency (cd ^) 


Period (min) 


Comments 


J1457 


18.4888 + 0.0017 


77.92 


Orbital period [cOq) 




36.9740 + 0.0005 


38.95 


2(Vo 




71.9 (nightly mean error: 0.5) 


20.0 


NRP, low amplitude, non-stable 




135.2/144.3/147.9 (nightly mean error: 0.8) 


10.7/9.9/9.7 


NRP, non-stable 


BWScl 


~ 16.50 + 0.01 


87.27 


Quiescent superhump, slightly non-stable 




18.40811 + 0.03 


78.23 






32.98 ± 0.01 


43.66 


Harmonic of quiescent superhump 




36.81622 + 0.03 


39.11 


Icoo 




50.5 + 0.5 


28.5 


period signal, weak 




69.55 + 0.03 


20.70 


NRP, Ci>i, non-stable 




103.55 + 0.03 


13.90 


coj - 2ct)o, non-stable 




140.37 ± 0.03 


10.26 


NRP, a>2, non-stable 




121.99 + 0.03 


11.80 


0)2 - <^or non-stable 




153.0 + 0.5 


9.4 


period signal, transient 




307.0 + 0.5 


4.7 


Harmonic of signal at 153 c d~^, transient 



Table 3.2: Summary of frequencies found in SDSS J1457+51 and BW Sculptoris. 
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during five of the observing nights, a broad, low-amplitude peak appeared at 135 
c d"\ along with a signal at 72 c d"\ neither of which are of orbital origin. 

In order to study these signals in more detail, a coherent power spectrum was 
constructed by concatenating the data from nights. However, this did not pro- 
duce clean signals, even though the nightly power-spectrum windows were al- 
ways clean. The broadness and complexity of the signals indicate a slight shift in 
amplitude and frequency on the time scale of a few nights and/ or an internal fine 
structure unresolved by our observations. Therefore, analysis was performed on 
power spectra from separate nights in comparison with the overall mean power 
spectrum. 

During the first observing night, a peak at 142.6 + 1.6 c d~^ was detected. Over 
the nights that followed, the signal was found at higher and higher frequencies 
and two weeks later, a peak was visible, at 147.8 + 0.4 c d~^. When combining data 
collected over 43 nights, there is a broad power excess around the frequencies, 135.2 
cd^^ (low amplitude), 144.3 cd"^ and 147.9 cd"^. This splitting is evident also in 



the mean spectrum shown in Figure 3.1(a) and is always seen when combining 
nights from the start and the end of the campaign. When studying the nightly 
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Frequency (c/d) 

Figure 3.4: The two top panels show the power spectrum of J1457 from the same set of 
nights shown in the mean spectrum (see Figure \3.1(a)i . The two bottom frames show 
models that was constructed from the same sampling pattern as the data. The following 
frequencies was included in the model; 18.48 cd~^(coo), 36.97 cd~^(2a)g), 71.9 cd'^, 135.2 
cd~^ , 144.3 cd'^ and 147.9 cd'^ . All main signals along with their line widths seen in 
the data can he reproduced fairly well by the model, indicating that there is power excess at 
these frequencies. However, if any of these frequencies are removed from the model, there is 
no longer a good match between model and data. 



power spectra in this range, there was (in general) only one peak present at the 
time. However, during one observing night, three peaks were seen simultaneously 
at 135.6 + 0.4cd-\ 144.6 + 0.3 c d^^ and 148.9 + 0.3 c d~i (Figure 1331). The nightly 
mean error for any peak appearing at 135 cd^^ - 148 cd"^, is about 0.8 cd"^. 

The upper two frames in Figure 13.41 show the combined power spectrum of 
the six nights of best quality. The spectrum spans 43 nights in total, and shows 
the regions around the signals at 72 cd~^ and 135 cd~^ - 148 cd~^. The bottom 
two panels shows model power-spectra constructed from sine waves at 18.48 c d"^ 
(wo), 36.97 cd-i (2wo), 71.9 cd^^, 135.2 cd'^, 144.3 cd^^ and 147.9 cd^^, and the 
same sampling pattern as the data. The model is able to re-construct the overall ap- 
pearance and widths of the signals seen in the data reasonably well, indicating that 
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Figure 3.5: Significance levels (red = 99%, green = 99.9% and blue = 99.997o) in the power 
spectrum constructed from the combined data of the six nights of best quality. The peak at 
72 cd"^ appears to be below the significance level, but since it is present in many individual 
data sets, it is still of interest. 



there is power excess at these frequencies. If any of these frequencies are omitted 
from the model, the data power spectrum cannot be reproduced. The signal at 71.9 
c d^^ is about half that of 144.3 c d~\ but not exactly, indicating that the signals are 
not constant in amplitude and phase. 

Figure 13.51 shows a power spectrum of the six nights of best quality, where the 
significance levels are marked as horizontal lines (red= 99%, green = 99.9% and 
blue = 99.99%). The peak at 72 c d"^ is below the significance level. However, since 
it is repeatedly seen in many data sets, and since a signal at this frequency has also 
been detected in similar targets (see Section l3.2.4l for BW Sculptoris), it should still 
be considered as significant. 

In Figure l3!6l all peak frequencies detected in the range spanning 130 c d^^ - 150 
c d"^, are plotted as a function of time. Horizontal lines mark the mean frequencies 
of the power excess found in the combined power spectrum. The figure shows 
that the signals do not appear to be intrinsically stable. However, as indicated 
in Figure 13.41 there could be an underlying stable structure at those frequencies 
indicated by the solid lines, varying in amplitude, that our sparse data are not able 
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Figure 3.6: The distribution of frequencies found in J1457, in the range 130 cd - 150 
cd~^, plotted as a function of time. 



to resolve in detail. 

I conclude that there are persistent periodic signals at the mean frequencies 71.9 
cd~\ 135.2 cd~\ 144.3 cd"^ and 147.9 cd""\ which are non-commensurate with 
the orbital period. The signals at 72 c d"^ and 144 cd~^ are not exactly in a 2:1 ratio 
(see Section [3.2.41 for the case of BW Sculptoris). For a complete summary of the 
frequency analysis, see Table 13.21 



3.2 BW Sculptoris 



BW Sculptoris (hereafter BW Scl) is a 16th magnitude blue star which was found 
to coincide with RXJ2353-0-3852 in th e ROSAT bright-source catalogue, and then 



identified as a cataclysmic variable by 



Abbott et al 



1993). 



Augusteijn & Wisotzki 



(|l997n independently discovered the star as a blue object in the Hamburg/ESO sur- 
vey for bright quasi-stellar objects (QSOs). These two studies established the very 
short orbital period of 78 minutes. In addition to broad and doubled-peaked H and 
He emission lines, BW Scl also shows very broad Balmer an d Lyman absorptions . 



signifying the presence of a WD of a temperature ^ 15000 K (jGansicke et al. 



2005) 



If half the visual light comes from such a white dwarf, then the WD has V - 17.3 
and My ~ 12, implying a distance of only 110 pc. This also agrees with the 
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large proper motion found in the USNO catalogue (105 msyr~^: lGirard et al.ll2004l) . 
These considerations (a nearby star of very short Porb)/ arid the possibility to study 
the underlying white dwarf, was the motivation to carry out observations of the 
star. 



3.2.1 Observations 

In total, BW Scl was observed for about 1000 hours spread over about 200 nights, 
mainly using the globally distributed telescopes of the Cen ter for Backyard Astro- 
physics fCBA: lskillman & Pattersonl l I993I: lpattersonlll998ai . A summary observing 
log is presented in Table 13.31 In order to maximise the signal and optimise the 
search for periodic features, usually no filter or a very broad filter spanning the 
range 4000 A - 7000 A was used. Occasional runs were obtained in V and I band- 
passes to provide a rough calibration, and BW Scl was found to remained within 
~ 0.3 mag of V - 16.6, throughout the campaign. In order to study the periodic 
behaviour, data was obtained densely distributed in time, preferably with contri- 
bution from telescopes widely spaced in longitude (in order to solve problems as- 
sociated with daily aliases). Most of the analysis was based on long time series 
from telescopes in New Zealand, South Africa, and Chile, and hence not afflicted 
by aliasing problems. 



3.2.2 The 1999 Campaign 

The upper frame of Figure 13^71 shows the light curve from one night during 1999. 
The general appear ance is typical of all nights, as well as seen in p reviously pub- 



lished light curves (|Augusteijn & Wisotzki 



1997 



Abbott et al. 



19971). As those pa- 



pers demonstrated, a 39-minute periodicity is always present, suggesting that the 
fundamental period is actually double this, i.e. 78 minutes (confirmed by spec- 
troscopy). The middle frame shows the 1999 double-humped orbital light curve, 
similar to that of J1457 (see Figure 3.2(c)| |. 

The average nightly power spectrum, the incoherent sum of the 11 nights of 
best quality, is shown in the bottom frame of Figure 13.71 Power excesses near 72 
cd-^ and 143 cd^^ are evident. In order to study these higher-frequency signals 
in more detail, adjacent nights were added together, and the coherent power spec- 



74 



3.2. BWSCULPTORIS 



Year 


Spanned 

(days) 


Observer 


Telescope 


Nights 

(hours) 


1999 


13 


Kemp 


CTIO 91cm 


12/61 


2000 


38 


McCormick 


Farm Cove 25cm 


7/35 






Rea 


TsTplson 35rm 


2/9 


2001 


77 


Rea 


n 


16/87 






Kemp 


CTIO 91cm 

X xv_y A^XX L 


14/84 






Woudt 


SAAO 76cm 


1/4 


2002 





Kemp 




1/4 


2004 


7 


Monard 

X T XVyX L^pCX Wf. 


Pretoria 25cm 


7/38 


2005 


55 


Rea 


Nelson 35cm 


16/66 






CViriofjp 

LX XL? LXV^ 


Auckland 35cm 

J. XU.\_XVXC1.X L^l. k-/\_/\>XX L 


12/48 






Retter /Liu 




6/32 






Monard 


Pretoria 25cm 

X X ^ CV^X X^pX *-b\-/ V^XX 1. 


4/22 










3/10 


2006 


90 


Rea 


X >l V^XCJWX I v^v^'V^XX I 


26/112 






Monard 

X T xvyx l^pCX 


Pretoria 35cm 

X X ^ cvyx xcx \^\-/^xx i- 


3/21 






McCormick 


Farm Cove 25cm 


5/15 


2007 


21 


Rea 


Nelson 35cm 


7/28 






Richards 


Melbourne 35cm 


1/3 






Allen 


Blenhem 41cm 


1/3 


2008 


71 


Monard 


Pretoria 35cm 


22/88 






Rea 


Nelson 35cm 


13/90 


2009 


50 


Rea 


It 


15/94 






Monard 


Pretoria 35cm 


11/62 



Table 3.3: Summary observing log for BW Sculptoris. 
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Figure 3.7: BW Scl in 1999. The upper frame shows an 8-hour light curve, dominated by 
the orbital wave. The middle frame presents the mean orbital light curve, showing a double- 
humped waveform. The lowest frame shows the mean nightly power-spectrum, averaged 
over the 11 best nights. Significant features are labelled with their frequency in cd~^ (with 
an average error of± 0.7). 



triun was constructed. However, even though the 72 c and the 143 c d~^ signals 
were always present, they were always broad, complex, and slightly variable in 
frequency (similar to J1457). This is usually a sign that the actual signals violate the 
assumptions of Fourier analysis; constancy in period, amplitude and phase. 

3.2.3 The 2001 Campaign 

During 2001, observatories at three longitudes in Chile, New Zealand and South 
Africa, contributed with good nightly coverage of BW Scl. A very low flickering 
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Figure 3.8: BW Scl in 2001. The upper two frames show the mean nightly power spectrum, 
averaged over the 15 best nights. Features of obvious or possible significance are labelled 
with their frequency in cd~^ (with an average error of+ 0.7). 
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background is shown in the light curves. Variations as small as ~ 0.002 magnitudes 
could be detected due to the very extensive campaign that was carried out. 

The upper frames of Figure l3!8l shows the nightly mean power spectrum, with 
significant signals marked to an accuracy of 0.5 cd"^. Both the orbital period and 
the signal at 72 c d"^ were present along with two other signals at low frequency: a 
powerful signal at 16.6 cd~^, and a weak signal at 50.5 cd~^. At higher frequency, 
signals are detected at 153 c d"^, 307 c d"^, and 724 c d"\ though the latter is likely 
to be caused by instrumental effects. Many telescopes have worm gears which 
turn with a period of exactly 120 sidereal seconds, and this period was reported in 
research on many types of stars during 1960 - 1990, i.e. during the photolectric- 
photometer era. CCDs are much less prone to this error. However, since 724 + 2 
cd~^ corresponds to 120 sidereal seconds (to within the measurement error), we 
interpreted the signal to be caused by this instrumental effect. 

The signal at 153 cd~^ can be interpret as a pulsation frequency, with a signi- 
ficant first harmonic. Examination of individual nights showed this signal to be 
slightly variable, at least in amplitude. This behaviour was seen on 7 of the 15 
good-quality nights. Since the orbital frequency is known precisely, and since its 
photometric signature is powerful and constant, the first harmonic (and also the 
second harmonic when detected) was subtracted from the central 12-night time- 
series, prior to analysis. The resultant power spectrum at low frequency is seen 
in the bottom panel of Figure 13.81 A weak signal appears at the orbital frequency, 
and stronger signals at 16.34/16.64 c d"^ and 32.98 c d"^ (with an accuracy of + 0.01 
c d~^). It seems likely that these are superhump signals in the quiescent light curve. 
Specifically, the 16.34/16.64 pair can be interpreted as signifying an underlying fre- 
quency of ~ 16.5 cd \ as such splitting can be produced by amplitude and/or 
phase changes. The precession frequency Q can be expressed as 16.50 = 0)^-0., 
where cOq is the orbital frequency (implying that 32.98 = 2{a)o-D.)). The precession 
frequency Q is then equal to 1.9 c d^^ . 

The primary s ignal near 16 . 5 cd~ ^ has been seen before. It was the dominant 



signa l reported by lAbbott et al.l (|l997l) , but was then discounted by lAugusteijn & Wisotzki 



( 1997n as probably the result of cycle-count error. The data presented here certainly 
have no ambiguity in cycle count and reveals this signal very clearly. 



78 



3.2. BWSCULPTORIS 



—1 — I — 1 
69 


1 r 1 1 1 1 r 

7 


1 1 1 1 1 1 ' 1' 1 I 1 1 1 1' 1 1 1 1 1 1 1 1 1 >■< ' 

BW SCL 2009 






NIGHTLY AVERAGE 
(19 nights) 


).,!,■ 


139 

104- 

IvW 

L,l ) I 1, \ 1 .,1 


9 

1 1 1 1 1 1 1 1 1 1 1 t 1 ) 1 1 1 1 ■ I 1 1 1 1 1 



40 80 120 160 200 240 280 320 360 400 
FREQUENCY (cycles/day) 




40 80 120 160 

FREQUENCY (cycles/day) 



Figure 3.9: BW Scl in 2009. The top frame shows the mean nightly power spectrum, 
averaged over the 19 best nights, features are labelled with their frequency in cd~^ (+ 0.6). 

Other frames show power spectra of ~ 5 day intervals with particularly dense coverage, 
with significant features labelled (+ 0.03). 
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3.2.4 The 2009 Campaign 

The upper frame of Figure 13.91 shows the nightly power spectrum, averaged over 
the 19 nights of coverage obtained in 2009. Signals near 70 cd~^ and 140 cd~^ are 
present, as well as weaker signals near 104 cd"^ and 122 cd"^. In the three lower 
frames, intervals of particularly dense coverage are shown, each with a frequency 
resolution of + 0.04 cd"^. In these frames, the 70 cd^^ and 140 cd^^ signals (both 
normal features of BW Scl) show their variability in frequency and amplitude (sim- 
ilar to J1457). If we define cu2 = 140 c d"\ the frequencies at 104 c d^^ and 122 c d~^ 
are displaced by exactly C02 - Iwq and 0^2 - of the orbital frequency. 



3.2.5 SUPERHUMPS IN BW SCULPTORIS? 



The signal at 87 minutes in BW Scl, displaced by ~ 1 1 % from Poj.^ seen in the 2001 
power spectrum (Figure |3^ , is a transient but repeating feature in the light curve. 
This feature was always strong when detected, but moved slightly in frequency, 
even on timescales of a few days. This qualitatively describes a common super- 
hump, which is a well-known feature in CV light curves (see Section [L22]|. How- 
ever, there are a few substantial differences: superhumps are found in outburst 
states, typically at My ~ 5, and not at quiescence when My ~ 12. Also, super- 
humps occur with fractional period excesses, (Pgh ~ ^orb)/^orb/ near 3%, not 11%. 
Systems with quiescent superhumps are rare, but has been reported be fore. In 



AL Co m, a nearly identical signal was reported and extensively discussed by 



Patterson et al 



( 1996h. A simi : 



Patterson et al. 



ar signal has also 



in SDSS J1238-03 (lAviles et al. 



jeen reported in CP Eri (s ee Appendix Al of 



20q2bh, V455 And (lArauio-Betancor et al. 



20051) , and possibly also 



2OIOI) . Except for CP Eri (which is excluded from 



consideration because of its heliu m content), all t hese stars are all period-bouncer 



candidates (see Tables 3 and 5 of 



Patterson 



20111) . These candidates were chosen 



due to their low donor-star masses (or low q = Mj/Mi). A possible account of 
how stars of very low a might be a ble to manufacture quiescent superhumps has 
been given in [Patterson et al.l (|l996l) , and includes the idea that a low donor mass 
implies a larger Roche lobe surrounding the disc. Weak tidal torques could then 
allow the quiescent disc to extend to the 2:1 orbital resonance, where an eccentric 



instability could drive a fast prograde precession (viz., at Q = 1.9 cd ), resulting 
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in a superhump with 0) = cOg-Q. 



3.3 Interpretations as Non-Radial Pulsations 

Both BW Scl and J1457 are quiescent CVs of very low accretion luminosity. Their 
spectrum show evidence of the primary WD, as do their light curves which con- 
tain rapid, non-commensurate signals. This is the general signature of the GW Lib 
stars, where the periodic signals are believed to represent the non-radial pulsations 
of the underlying white dwarf. However, no proof of this has ever been found, not 
for GW Lib or any other of the 10 - 15 members of the class. In order to explain 
these signals, two other hypotheses deserve consideration: first, quasi-periodic os- 
cillations (QPOs; see Section [LZ2)| arising from the accretion disc, and second, the 
spin frequency of a magnetic WD (DQ Her modulations; see Section [1.2.31 1. QPOs 
are typically seen as very b road peaks in the pow er spectrum, dv/v -0.5 (see for 



instance Figure 11 and 12 in [Patterson et al.ll2002a ), while the signals at 70 c d ^ and 



140cd-i (cui and a;i),inbothBW Scl and J1457 have 0.01. Therefore, QPOs 

are most likely not the origin of these pulsations. A WD spinning with a period 
of 20 minutes could possibly explain the signals seen at ~ 70 c d^^ and 140 c d~^ in 
both objects (and could possible also explain the orbital sidebands seen in BW Scl). 
However, such modulations usually result in a much more coherent signal, while 
those seen in BW Scl and J1457 are broad and shift in frequency and amplitude on 
a time scale of a few days. Also, in both BW Scl and J1457, 0)2 is not exactly Ixco^. 

With these considerations taken into account, I conclude that J1457 and BW Scl 
most likely are new members of the GW Lib class. Common to all known GW Lib 
stars is that they are low-M CVs, showing low-amplitude and non-commensurate 
periodic signals in their light curves. The signal is roughly constant over a few 
days, but is slightly shifting in frequency and amplitude on longer timescales. Also, 
these signals often show strong first harmonics. 
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3.4 Summary 

1. Two more CVs with very low accretion rates have shown rapid non- commen- 
surate signals in quiescence, which makes them new members of the GW Lib 
class. Both J1457 and BW Scl show a complex power spectrum with main 
signals near 10 and 20 minutes. 

2. The pulsation frequencies in both stars wander by a few percent on a times- 
cale of days. In addition, the power-spectrum constructed from multiple- 
night light curves, show broad peaks, which might be due to the frequency 
wandering, or a fine structure not resolved by our data. 

3. BW Scl shows several peaks displaced from the main pulsation frequency by 
exactly C02 + (Oq (02 + 2(0^, where cOq is the orbital frequency « 18.4 c d"\ 
2a>o w 36.8 c d~^ and CO2 ~ 140 c d~^ Similar behaviour is evident also in other 
GW Lib stars (SDSS J1507+52, V386 Ser and SDSS J1339+48). The origin of 
this phenomenon is still unknown. 

4. The orbital light curves of both stars show double-humped waves. From 
these waves, precise periods are found at Porb = 0.05432392(2) days for BW 
Scl, and 0.054087(5) days for J1457. Similar systems showing signals non- 
commensurate with the orbital period, such as SDSS J1339+4847, SDSS J0131- 
0901 and SDSS J0919+0857, all have orbital periods at « 80 minutes. BW Scl 
sometimes shows a transient wave with a period, Pgh = 0.060 days, which 
is interpreted as a quiescent superhump. It thereby joins a small group of 
stars who manage a superhump at quiescence, all of which are likely to have 
very low mass ratios. This might arise from an eccentric instability at the 2:1 
resonance in the disc. 
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The Recurrent Nova T Pyxidis 

- Orbital Period and System Parameters - 



All work presented in this chapter w as carried out tos,e ther with Christian Knigge and 



Danny Steeghs and was published as \Uthas et al\ uQlW- In particular, the photometric 



ephemeris was calculated by Christian Knigge. I also want to thank Joe Patterson for valu- 
able comments and for kindly providing us with the a compilation of recent photometric 
timings for T Pyx. Also Gijs Roelofs is thanked for obtaining part of the data. 

T Pyx is a luminous cataclysmic variable in which the donor star is transferring 
mass at a very high rate onto its white-dwarf companion, resulting in unusually 
frequent thermo-nuclear runaways (TNRs) on the surface of the primary. Between 
the years 1890 - 1967, T Pyx has undergone five such nova eruptions, with a re- 
currence time of about 20 years between the eruptions; it was therefore classified 
as a member of the recurrent nova subclass. However, the last eruption was in 
1966, which means that T Pyx has now passed its mean recurrence time by more 
than 20 years. The eruptive behaviour of RNe in comparison with classical novae 
is thought to be due to a high mass-transfer rate in combination with a massive 
prir nary white dwarf (see Section [L2?T)- 



Patterson et al. 



(|1998n carried out an extensive photometric study of T Pyx and 
found a stable, periodic signal at P = 1.83 hours that was interpreted as the likely 
orbital period. This would place T Pyx below the CV period gap and suggests 
a donor mass around M2 ~ 0.1 Mq. This is surprising. According to standard 
evolutionary models, a CV below the period gap should be faint and have a low 
accretion rate driven primarily by gravitational radiation (GR). Yet T Pyx's quies- 
cent luminosity and status as a RN both imply that it has a high accretion rate of 
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Figure 4.1: Accretion rates of cataclysmic variables plotted as a function of their orbital 
periods. The approximate position of T Pyxidis is marked with a black star, showing that 
the system clearly has a higher accretion rate comp ared to other systems at the same orbital 



period (adapted from figure by 



Howell et ah 



2003 



> IQ-^ Mq yr-^ ^Patterson et al. 



1998 



Selvelli et al. 



20081) . This is about two orders 



of magnitude higher tha n expected for ord inary CVs at this period (illustrated in 
Figure 14. 1[ adapted from 



Howell et al 



200l|). 



Assuming that the determination of the photometric orbital period is correct, 
the existence of T Pyx is interesting for at least two reasons. First, unless the system 
is seen in a transient evolutionary state, its lifetime would have to be very short, 
T ~ Mj/Mj £ 10 Myrs. This would imply the existence of an evolutionary channel 
leading to the fast destruction of at least some short-period CVs. Second, TNRs 
frequent enough to qualify as RNe are thought to be possible only on high-mass 
accreting WDs [Mi ^ 1 M©). Moreover, RNe are the only class of novae in which 
the WD is expected to gain more mass between eruptions than it loses during them. 
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This would make T Pyx a strong candidate Ty pe la supernova pro genitor. 



However, the recent study of the system by 



Schaefer et al 



|2010D(seealso 



Selvelli et al 



2008h suggests, first, that T Pyx is, in fact, in a transient evolutionary state, and. 



second, that, integrated over ma ny nova erupt i ons, i ts WD does lose more mass 



than it gains. More specifically. 



Schaefer et al. 



Il2010j) suggest that T Pyx was an 



ordinary cataclysmic variable until it erupted as a nova in 1866. This eruption 



trigg ered a wind-driven supersoft X-ray phase (as first suggested by iKnigge et al. 



20001) , resulting in an unusually high lum inosity and a c cretio n rate. However, un- 



like in the original scenario proposed by iKnigge et al.l (|2000h . the supersoft phase 
is not self-sustaining, so that the accretion rate has been declining ever since the 
1866 nova eruption from M ~ 1 0^^ M© yr^^ to 1 0^^ Mq yr~^ As a result, T Pyx has 
faded by almost 2 magnitudes since the nova eruption (jSchaefer et al.ll2010l ). Based 
on this, and th e fact that T Pyx has a lready passed its mean recurrence time by more 



than 20 years, ISchaefer et al.l (|2010[) argue that T Pyx might no longer even be a re- 
current nova. If these ideas are correct, T Pyx is not a viable SN la progenitor, and 
its remaining lifetime can be substantially longer than a few million years. How- 
ever, if all its ordinary nova eruptions are followed by relatively long-lived (> 100 
yrs) intervals of wind-driven evolution at high M, its secular evolution may nev- 
ertheless be strongly affected, with significant implications for CV evolution more 
generally (see also Section |4!4]| . 

A key assumpti on in virtually all of t hese arguments is that the photometric 



period measured by 



Patterson et al 



19981) is, in fact, the orbital period of the sys- 



tem. So far, there has orv ly been one attempt to obtain a spectroscopic period for T 



Pyx, by IVogt et al. 



19901) ■ who reported a spectroscopic modulation with P - 3.44 



hours. Such a long orbital period above the CV period gap would be much more 
consistent with the high accretion rate found in T Pyx. In this study, we present the 
first definitive spectroscopic determination of the orbital period of T Pyx, showing 
that it is, in fact, consistent with Patterson et al.'s photometric period. We also use 
our time-resolved spectroscopy to estimate the main system parameters, such as 
the velocity semi-amplitude of the white dwarf (iCi), the mass-ratio (q), the masses 
(Ml and Mj) and the orbital inclination (i). Finally, we discuss the implications of 
our results for the evolution of T Pyx and related systems. 
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4.1 T Pyx - Observations and Their Reduction 
4.1.1 VLT Multi-fibre Spectroscopy 

Multi-fibre Spectroscopy of T Pyx was obtained during five nights in 2004 and 2005 
with the GIRAFFE /FLAMES instrument mounted on the Unit Telescope 2 of the 
VLT at ESQ Paranal, Chile. The data were taken in the integral field unit mode. The 
total field of view in this mode is about 11.5" x 7.3" and thus covers most of T Pyx's 



10" diameter nova shell ( Williams 



119821 ). We used the fibre system ARGUS, which 
consists of 317 fibres distributed across the field, of which 5 point to a calibration 
unit and 15 points to sky. The grating order was 4, which gives a resolution of 
R=12000. The wavelength range was chosen between 4501 A to 5078 A, so that the 
emission-line spectrum would include the Bowen blend at 4645 A - 4650 A and 
the Hell at 4686 A. With this setup, the dispersion is 0.2 A/ pix, corresponding to 
X 12.5 kms"Vpix. The full widths at half-maximum (FWHMs) of a few spectral 
lines obtained simultaneously with the science spectra from fibres pointing to the 
calibration unit indicate that the spectral resolution is about 0.4 A. A log of the 
observations can be found in Table l4T] 

The initial steps in the data reduction were performed using the ESQ pipeline 
for GIRAFFE. The pipeline is based on the reduction software BLDRS from the 
Observatory of Geneva. The basic functions of the pipeline are to provide master 
calibration data and dispersion solutions. The pipeline also provides an image of 
the reconstructed field of view, which can be used to associate a specific fibre to a 
given object. In order to extract the spectrum from the desired fibre, and to correct 
for the contribution from the sky background and cosmic rays, the output from the 
pipeline was processed further in IRAF. The PSF of our target, T Pyx, is spread out 
over several fibres and 6 single fibre spectra containing significant target flux were 
extracted and median combined after weighting each spectrum by the mean flux in 
the region 4660 A - 4710 A (covering Hell at 4686 A). Cosmic rays were removed, 
by first binning the spectra over 7 pixels (the cosmic rays have a typical width of 
2-6 pixels). The smoothed spectra were then subtracted from the corresponding 
fibre spectra to only leave the residuals and the cosmic rays. Cosmic rays were 
then removed and the residuals added back to the target spectra. Fibres containing 
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Date 


Tel/Inst. 


No. of Exp. 


Exp. (s) 


041210 


VLT/GIRAFFE 


28 


180 


041224 


VLT/GIRAFFE 


28 


180 


050127 


VLT/GIRAFFE 


29 


180 


050128 


VLT/GIRAFFE 


28 


180 


050131 


VLT/GIRAFFE 


29 


180 


080316 


MAGELLAN/IMACS 


37 


240 


080317 


MAGELLAN/IMACS 


62 


240 


080318 


MAGELLAN/IMACS 


62 


240 


080319 


MAGELLAN/IMACS 


39 


240 



Table 4.1: Log of the observations (the date is according to UT time at the start of the night). 

the sky were extracted and combined to create a master sky spectrum. Finally, the 
master-sky was subtracted from the combined science spectrum. 

4.1.2 Magellan Long-slit Spectroscopy 

T Pyx was observed again during four nights in March 2008, this time on the 6.5 
meter Baade telescope Magellan I at Las Campanas, Chile (see Table 14.11 for a log 
of the observations). Long-slit spectroscopy was obtained with the instrument 
IMACS using the Gra-1200-17.45 grating with a 0.9" slit, resulting in a dispersion 
of 0.386 A/pix, corresponding to « 25 kms~Vpix. We estimate that the spectral 
resolution of these observations is about 1.5 A, as measured from the FWHMs of a 
few spectral lines in the arc-lamp spectra. The overall spectra span over four CCDs 
and cover a total wavelength range of 4000 A - 4800 A. All frames from the four 
CCDs were treated separately during both the 2D and ID reduction steps. 

The spectra were reduced in IRAF using standard packages. A master bias 
produced from combining all bias frames obtained during the four nights was sub- 
tracted from the science frames, and the overscan region was used to remove the 
residual bias. A master flat-field, corrected for illumination effects, was produced. 
The spectra were then flat-field corrected and extracted. Line identifications of the 
ThArNe lamp spectra obtained during the nights were carried out with the help of 
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the NO AO Spectral Atlas Central and the NIST Atomic Spectra Database. The spec- 
tra were cosmic-ray rejected using the method described previously for the VLT 
data set. A flux calibration was done using data of the white-dwarf flux-standard 
star EG274, obtained during the same nights as the targe t spectra. Also, exti nction 



Hamuy et al. 



( 1992 ) were 



data from the site and flux reference data of EG274 from 
used. 

Analysis of the reduced VLT and Magellan data were carried out using the 
packages MOLLY and DOPPLER, provided by Tom Marsh. 



4.2 Analysis of T Pyx 



4.2.1 The Overall Spectrum 



T Pyx has a high-excitation emission-line spectrum that is unusual when compared 
to systems with a low M, but the overall spectrum is similar to that for nova-like 
variables, supporting the idea that T Pyx has a high accretion rate. Its brightness 
has been found to be fading for the last cent ury, indicating a de crease in M, and 
its 2009 magnitude in blue was B = 15.7 mag jSchaefer et al.ll2010l) . Double-peaked 
lines originating from the disc are seen in the spectrum, and the bright disc out- 
shines any spectral signature of the primary WD and the donor star. The strongest 
feature is the doubled-peaked Hell line at 4686 A. Double-peaked lines identified 
as Hel at 4713 A, 4921 A and 5015 A, as well as the Balmer lines, are also present 
(Figure l42l and Figure 1431) . 

The Bowen blend at 4640 A - 4650 A is clearly visible and consists of several 
lines of NIII and CIII. The Bowen blend is also seen in other CVs and in low-mass X- 
ray binaries with high- accretion rates. In some c ases, it has been related to emission 
from the donor star (eg. lSteeghs & Casaresl2002l '), where it is thought to be produced 
by a fluorescence process as the front side of the donor is strongly irradiated by the 
ho t accretor. The Bowen blen d was clearly visible in the spectrum of T Pyx obtained 
by lMargon & DeutschI (|l998l) . and, based on this, we were hoping to detect narrow 
components associated with the donor. However, no narrow donor star features 
were found in the blend. 

Data analysis was carried out for all the strongest lines, but our final results are 
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based on analysis of the Hell line at 4686 A and the Hel line at 4921 A. 



4.2.2 The Orbital Period 



The Photometric Ephemeris 



The key goal of our study is to obtain a definitive determination of T Pyx's orbital 



period based on time-resolved spectrosc opy. More specifical 



the photometric modulation reported by 



The photometric ephemeris is given by 



Patterson et 



Patterson et al 



(|l998h . 



y, we wish to test if 



(19981) is orbital in nature. 



Minimum light = HJD 2446439.428 + 0.0762233£ + 3.5 x lO^^E^ 



The quadratic term is highly significant and indicates a very high period deriv- 
ative of P ^ 9 X 10"^° (this corrects a typo in the original paper). The period is 
found to be increasing, which is not expected for a system that has not yet passed 
the minimum period. If the photometric signal is orbital, the current evolutionary 
time-scale of the system is only t ~ P/P -'3x10^ yrs. Such a short period change 
time-scale is highly unusual for a CV. 

The Center for Backyard Astronomy (CBA) has continued to gather additional 
timings in the period 1996 - 2009. This new photometry, along with an updated 
ephemeris, will be published in due course (Patterson et al., in preparation). How- 
ever, since this new data set is much more suited to a comparison against our spec- 
troscopic observations (since it overlaps in time), Joe Patterson has kindly made 
this latest set of CBA timings available to us. Our own preliminary analysis of all 
the existing timings yields a photometric ephemeris of 



Minimum light = HJD (UTC) 2451651. 65255(35) + 0.076227249£(16) + 2.546(54)x 
10-ii£2. 



The numbers in parenthesis are the uncertainties on the two least significant digits 
and were estimated via bootstrap simulations. Taking the period change into ac- 
count, the average period at the time of the optical spectroscopy was P = 0.076228860+ 
3.0x 10"^ d(/ = 13.1183912 + 5.2 X IQ-^cd-i). 
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Figure 4.2: The normalized and averaged spectra constructed from 140 individual expos- 
ures obtained at the VLT telescope. 
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Figure 4.3: The flux-calibrated average of 200 spectrum observed with the Magellan tele- 
scope. 
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100 - 




13,11 13.115 13.13 13.135 13.13 



Frequency c/d 

Figure 4.4: Power spectrum of the Hell radial-velocity data, obtained from the combined 
VLT and Magellan data sets. The two bottom panels are zooms of the panel above. The 
vertical grey line in the bottom panel represents the best period found by Patterson et al. 
1998, along with its error. The red line represents the recalculated period using the latest 
photometric ephemeris (see Section \4.2.2i where the period derivative has been taken into 
account (Patterson et al. in preparation). 

This new ephemeris confirms that the period continues to increase at a fast rate. 
More specifically, the current estimate of the period derivative is P = (6.68 + 0.14) x 
10"^", corresponding to a timescale ofP/P=:3xlO^ yrs. It is important to take 
this period derivative into account when comparing spectroscopic and photometric 
periods. 
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The Spectroscopic Period 

In order to establish the orbital period spectrosco pically, we carried out a r adial- 
velocity study using the double-Gaussian method Schneider & Young ( 1980 ). This 
method is most effective in the line wings, which are formed in the inner regions 
of the accretion disc. As discussed further in Section I4.2.4[ we tried a variety of 
Gaussian FWHM's, as well as a range of separations between the two Gaussians. 

The results shown in this section are for a FWHM = 450 kms~^ and a separa- 
tion of 300 kms~^ and are based on the full combined VLT + Magellan data set, 
spanning the time period December 2004 to March 2008. The radial-velocity study 
was carried out for several spectral lines but for the purpose of this section, we 
will focus exclusively on the radial-velocities derived from the strong Hell line at 
4686 A, since this line provides the most accurate results. The power spectrum 
( Scarglelll982h of this spectroscopic data set is presented in Figure l44l The top 
panel, covering a broad frequency range, confirms that there is a strong signal near 
the photometrically estimated 13 c d~^. No signal i s found at P = 3.44 hours, cor- 
responding to the period found by lVogt et al.l (|l990l) . The middle and bottom pan- 
els show close-ups of narrower frequency ranges around this power spectral peak. 
Note, in particular, that the red line in the bottom panel marks the photometrically 
predicted frequency for the epoch corresponding to the mid-point of our spectro- 
scopic data set. For comparison, we also show the frequency of the linear eph- 



emeris determined by iPatterson et al.l ((19981) for just the 1996 - 1997 timings, along 
with its error. The spectroscopic and photometric periods appear to be mutually 
consistent, but only if we account for the photometry period derivative. More specific- 
ally, the P-corrected photometric prediction lies extremely close to the peak of the 
most probable spectroscopic alias. 

This apparent consistency can be checked more quantitatively. The orbital fre- 
quency predicted by the photometry for the average epoch of our spectroscopy is 
/phot = 13.1183912+ 5.2 x 10"^ cd"\ where the error is once again based on boot- 
strap simulations. Similarly, we have carried out a bootstrap error analysis for the 
frequency determined from our spectroscopy. Since we are only interested in the 
agreement between the photometry and the spectroscopy, we only considered the 
spectroscopic alias closest to the photometrically determined frequency (this is also 
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the highest alias). This yielded /gpec = 13.118368+ 1.1 x 10"^ cd^^ The difference 
between these determinations formally amounts to just under 2a. We consider this 
to be acceptable agreement. 

The level of agreement between spectroscopic and photometric periods is im- 
portant. The spectroscopy on its own is sufficient to establish beyond doubt that T 
Pyx is a short-period system with Ppj-b -1-8 hrs. However, the fact that we can es- 
tablish reasonable consistency with the photometric data if (and only if) we account 
for the photometric period derivative suggests that (i) we can trust the much higher 
precision photometric data to provide us with the most accurate estimate of the or- 
bital period, and (ii) that the large period derivative suggested by the photometry 
is, in fact, correct. 



4.2.3 Trailed Spectrograms and Doppler Tomography 

In order to visualise the orbital behaviour of the spectral line profiles, trailed spec- 
trograms were constructed for Hell, Hel, Hj3 and the Bowen blend. These spec- 
trograms, shown in Figure 14.51 are phase binned to match the orbital resolution 
of our data (37 bins for the VLT data and 27 bins for the Magellan data) and are 
here plotted over 2 periods. All lines show emission features moving from the blue 
to the red wing, and localised emission is seen in the red wing at phase, (p « 0.6. 
Similar structures are seen in the spectrograms for both the Hell and the Bowen 
blend indicating that they originate in the same line-forming region, presumably 
the accretion disc. 



Doppler tomography ({Marsh & Hornd 19881) is an indirect imaging method where 



the emission-line profiles, depending on phase, are plotted onto a velocity scale 
(see Section l2.2.3|) . The method was used to visualise, in velocity space, the origin 
of the emission features found in T Pyx. In Figure 14.61 Doppler tomograms for the 
most prominent lines are presented. For the reconstruction of the tomograms, the 
systemic velocity y was set to zero (see Section l4.2.4|l . All lines show asymmetric- 
ally distributed emission from the disc. These asymmetric features are also seen in 
the trailed spectrograms (Figure |45l) . No emission can be connected to the bright 
spot, the accretion stream or the donor star. The map of the Bowen blend was con- 
structed from a composite of three lines, NIII at 4640.64 A and CIII at 4650.1 A and 
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(c) VLT, Hell (4686 A) (d) VLT, Hj3 (4861 A) 

Figure 4.5: Trailed spectrograms, phase binned to match the orbital resolution (37 bins for 
the VLT data and 27 bins for the Magellan data), are here plotted over 2 periods. Phases are 
plotted against velocities (in km s~^). 
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Figure 4.6: Doppler tomograms for the most prominent lines. The center of symmetry in 
each map is marked as a dot and the outer disc radius is plotted as a circle. Both axis are in 
kms~^. 
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4.2.4 A Radial-Velocity Study of T Pyx 
The Systemic Velocity 

We noticed early on in our analysis for both the VLT and Magellan data sets that the 
systemic velocity y, appeared to be varying with time. In order to confirm the real- 
ity of these variations, we extracted spectra from a few VLT /GIRAFFE /FLAMES 
fibres that were centred on knots in the nova shell surrounding T Pyx. The intrinsic 
radial velocities of these large scale knots are not expected to change appreciably 
over a time-scale of months, so they provide a useful check on the stability of our 
wavelength calibration. We found that the shifts seen in the radial velocity data 
for the central object are much larger than those seen for the knots, suggesting that 
these shifts cannot be explained by instrumental effects alone. However, since we 
have only five nights worth of VLT data spread over two months, our data set is 
too sparse to allow a careful study of trends in the apparent y velocities. For the 
purpose of the present study, we thus subtracted the mean nightly y from all radial 
velocities before carrying out any analysis. This should minimise the risk of biasing 
our results. 



The Velocity Semi-Amplitude 



In order to establish the velocity semi-amplitude, Ki , of the WD, further analysis 
was pursued on the radial-velocity data obtained from the double-Gaussian method 
A least-squares fit was made to the radial velocity curves by keeping the orbital 
period fixed, but allowing the three parameters, the systemic velocity, y, the velo- 
city semi-amplitude, Ki, and the phase, (pQ, to vary. Phase zero corresponds to the 
photometric phase of minimum light. The errors on the input radial velocity data, 
provided in MOLLY, were used to calculate the and then rescaled so that the 
reduced xl = ^- The technique was applied to several strong lines in both the VLT 
and Magellan data sets, but the radial velocity data measured from the Hell line at 
4686 A, and in particular from the Magellan datase t, gave the most reliable fits. The 



res ults were plotted i n a diagnostic diagram (see 



Shaffer 



19831: 



Thorstensen 



1986 



and 



Horne et al. 



19861). However, as can be seen in Figure l47l the key parameters 



{y, Ki and (po) did not converge convincingly for any combination of the FWHM 
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Figure 4.7: Diagnostic diagram for the Hell line, Magellan. None of the parameters (y, K 
and (Pq) are converging for any combination of the FWHM and separation. 

and separation. Normally, it is thought that the bright spot or other asymmetries in 
the disc tend to be responsible for distorting the diagnostic diagram. However, in 
T Pyx, the disc appears to be rather symmetric, and no contribution from the bright 
spot is seen (see Figure l431 l. We conclude that the radial velocity curves obtained 
from the double-Gaussian method are, most likely, not tracking the velocity of the 
primary WD.Thus, Ki, and (pQ carmot be estimated reliably using this method. 
In order to improve on the results obtained from the double-Gaussian method, 
we used two other techniques to obtain an estimate of Ki. First, we measured 
the velocity centre of the lines by fitting Gaussians to the individual spectral lines. 
Two Gaussians were fitted to the double-peaked lines, keeping the FWHM and 
the rest wavelength of the line fixed, but allowing the peak offsets and the peak 
strengths to vary. Radial velocity curves were then plotted for both blue and red 
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wings. As can be seen from the trailed spectrograms. Figure 14.51 the red wing (in 
particular for the Hell lines), does not show the same smooth orbital signal as the 
blue wing. Consequently, only the data from the blue wing was taken into account 
when measuring . Several lines from both the VLT and Magellan data sets were 
analysed. The data were phase binned to the orbital resolution, velocity-binned in 
a ~ 100 pixel wide region around the line, and the nightly mean y was subtracted. 
Hell at A4686 is the strongest line and was ultimately used to estimate our best-bet 
value of the velocity semi-amplitude of the WD: Ki = 17.9 + 1.6 kms"^. Figure |48l 
shows the radial velocities obtained for the Hell line in the Magellan data set. Data 
presented in black are phase-binned to the orbital resolution afforded by our time- 
resolved spectroscopy, while all individual data points are shown in grey. Over- 
plotted onto the radial velocity curve is the best sinusoidal fit, indicating a phasing 
close to zero, (pQ - -0.03 + 0.03. Figure |4!9] shows the final trailed spectrograms 
from the Hell at A4686 in the Magellan data set that was used to determine Ki . 

Second, Ki can also be measured from the Doppler tomograms, by locating 
the centre of symmetry of the disc emission that dominates the maps (this method 
is described in Section [2.3.21 1. This procedure yielded estimates for Ky that were 
consistent with those obtained from the Gaussian fits {lA < Ki < 19 kms"^). We 
thus adopt Ki = 17.9 + 1.6 kms~^ as our best-bet estimate. 

In closing this section, it should be acknowledged that our final estimate of Ki 
is subject to considerable systematic uncertainties. In particular, it is difficult to rule 
out that whatever is preventing the double-Gaussian method from converging may 
also bias the Gaussian fits to the line peaks and the centre of symmetry of the tomo- 
grams. Strictly speaking, our estimate of Ki should thus perhaps be viewed as a 
lower limit. However, we are reasonably confident that our measurements do trace 
the orbital motion of the WD in T Pyx. This is mainly because the presence of per- 
sistent double-peaked lines and ring-like structure in the tomograms implies that 
there is a strong accretion disc component to the line emission. This disc compon- 
ent must ultimately trace the motion of the WD, and our two preferred methods 
are designed to exploit this in as direct a fashion as possible, while simultaneously 
minimising contamination from asymmetrically placed emission regions. 
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Figure 4.8: Radial velocity curve of the Hell line, Magellan. Data in black is phase binned 
while all data are shown in grey. The best sin fit is over-plotted onto the phase-binned data. 



The Velocity at the Outer Disc Radius 

The peak-to-peak separations of the double peaked emission lines can be used to es- 
timate the projected velocity at the outer disc radius, sin i. We thus fitted Gaus- 
sians to the phase-binned, double-peaked H^, Hel and Hell lines in the Magellan 
data set. The peak-to-peak separations varies significantly with line, but the mean 
value is close to that for Hel. We therefore use the Hel peak-to-peak separation as 
our best-bet estimate of VRj.^^sinz and adopt half the full range as a rough estim- 
ate of the as sociated error, AVpeaij.to.peak = 290 + 26 kms"^ Taking the relation by 



SmakI (|l98ll ') into account (see Section [2.3.3L we find V]^.^^sin f = 1 38 + 1 5 km s ^ . 



4.3 Spectroscopic System Parameters 

If we are willing to assume that our estimates of Ki and vi^._,^sinz are valid, we can 
use our data to estimate several other system parameters. 
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Figure 4.9: Trailed spectrograms for the Hell line, Magellan. 



4.3.1 Mass Ratio 



Warnei 



The m ethod we use to estimate the mass ratio, q, is similar to that described in 
jl973 ). The basic idea is that, subject to fairly benign assumption, the ratio VRj._,^sin i/K^ 
should be a function of only q. More specifically, we assume that the disc is circu- 
lar, Keplerian, and, since T Pyx is a high-M system, that it extends all the way to 
the tidal radius (Equation 14.111 . It is then possible to show (see Appendix A for a 
detailed derivation of our modified relation) that 



= 1.2909 

Xl CI 



(4.1) 



Based on this, we estimate a mass ratio of q = 0.20 + 0.03. Given a mass ratio and 
disc radius, it is possible to estimate the expected phasing of the bright spot in 
the system, if we make the usual assumption that the bright spot lies at the point 
where th e ballistic accretion stream from the Li point meets the disc edge (see, for 



example. 



Patterson et al. 



2008 . 



2^. crying eu. .his calcuIaHon re, T pyx sug- 



gests that the bright spot should be located at orbital phase (p =^ -0.03. This is 
interesting, because it is consistent with the negligible phase offset we have found 
between photometric minimum light and the red-to-blue crossing of the WD radial 
velocity curve. Thus, based on phasing alone, the bright spot could be the source of 
the photometric modulation. However, there are also other possible interpretations 
for the photometric signature. For example, the orbital signal could be dominated 
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Figure 4.10: Doppler tomogram for the Hell line, Magellan. The location of the WD, Ki, 
the outer disc radius as well as the location of the donor star, K2, and the accretion stream 
is plotted onto the map. 



by a reflection effect, i.e. the changing projected area of the irradiated face of the 
donor as it moves around the orbit. This would explain why inferior conjunction of 
the secondary corresponds to minimum light. On the other hand, it raises the ques- 
tion why we do not see the signature of the irradiated donor in our spectroscopy. 
As noted above, there is no sign of the narrow Bowen fluorescence lines that one 
might expect in the case of strong donor irradiation. 



4.3.2 Component Masses 

Recurrent novae are expected to harbour a more massive WD than most ordinary 
CVs and even than most classical novae. In a classical nova, the mass of the WD 



vationally. Based on this, both 



is predicted theoretically to be ~ 1 (lYarort et al.ll2005), which is also seeri obser- 



Selvelli et al 



|2008Dand 



Schaefer et al. 



(|2ninh agree 



on a plausible mass r ange for the prim ary WD in T Pyx of 1.25 Mq - 1.4 Mq (using 



theoretical models by 



Yaron et al. 



20050 . From the theoretical point of view, a high- 



mass WD is needed to achieve a sufficiently short outburst recurrence timescale. 
In T Pyx, the bright accretion disc outshines any spectral signature from the 
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donor star, ar\d we can therefore not constrain the donor mass spectroscopically. 
How ever, we can us e the period-density relation for Roche-lobe filling secondary 



stars 



EggletonI (Il983h to set some constraints. 



As shown by 



Patterson et al 



(|2005bl ) andlKniggd (|2006h . the donor stars in or- 
dinary CVs below the period gap are inflated by approximately 10% due to mass 
loss, relative to ordinary main sequence stars of the same mass. Given T Pyx's pe- 
culiar evolutionary state, it is not clear if this level of inflation is appropriate for T 
Pyx, and we therefore adopt a conservative range of 0% - 20% inflation. In order 
to estimate the mass of the donor star, we thus tak e the theoreti c al ma in-sequence 



mass-radius relation from the 5 Gyr isochrone of iBaraffe et al.l (|l998l) . adjust the 
stellar radius to account for inflation, and then find the secondary mass that yields 
the correct density for T Pyx's orbital period. This yields M2 = 0.14+0.03 M©- Note 
that the inferred donor mass decreases with increasing levels of radius inflation. 
The corresponding mass of the WD is then Mi = 0.7 + 0.2Mo. 

Since this estimate of the WD mass is lower than expected for a recurrent nova, 
we can also turn the problem around. If the mass of the WD is > 1 Mq, as theoretical 
nova models imply lYaron et al.l (|2005t) , the mass of the donor becomes > 0.2Mo, if 
our estimate of the mass ratio is correct. 

All of these estimates should, of course, be taken with a grain of salt. In par- 
ticular, they rely on (i) the correctness of our measured and (ii) the assumption 
that the disc extends all the way to the tidal radius. Despite these uncertainties, our 
calculations highlight an important point: it is highly unlikely that the donor star in 
T Pyx is already a brown dwarf. For example, if we retain the assumption of max- 
imal disc size, Ki <7 kms"^ would be required in order for M2 < 0.07 M© for any 
Ml > 1 Mq. We thus rule out the possibility that T Pyx is a period bouncer in the 
usual sense, i.e. that it had already reached the minimum period for ordinary CVs 
70 min theoretically, or ^ 80 min observationally), in which case its secondary 



would now be well below the Hydrogen-burning limit 







Pyx, is, of course, a "period bouncer" in the basic sense that its orbital period is currently 
increasing. 
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Ki - 17.9 + 1.6 kms~^ velocity semi-amplitude 

V sin i = 1 38 + 1 5 kms"^ velocity at the outer disc radius 

(| = 0.20 + 0.03 mass-ratio (M2/M1) 

M2 = 0.14 + 0.03 Mq donor mass 

Ml = 0.7 + 0.2Mq WD mass 

f = 1 + 2 degrees inclination 



Table 4.2: The spectroscopic system parameters for T Pyx are presented here, where is 
found from the amplitude of the radial velocity curve of the HII line at 4686 A, and vs'mi 
is estimated from the peak-to-peak separation of the Hel line at 4921 A. The mass-ratio 
q is calculated from Ki and vsini, and i is calculated from q. Component masses were 
estimated from q and from the mass-radius relation accounted for an inflated donor radius. 



4.3.3 The Orbital Inclination 



The orbital inclination for T Pyx is thoug 
profiles and low radial velocity amplitude. 



it to be low due to the sharp spectral 



Shahbaz et al 



( 19971) suggested a lower 



'yx of i ~ 6° based on the peak-to-peak separ- 



Patterson et al 



1 19981 ) estimated an inclination 



limit of the orbital inclination for T ' 
ation of the Ha line in their spectra. 
of i ~ 10° - 20° due to the low amplitude of the orbital signal, and 
( 2008 ) estimated i ~ 20° - 30°. 



Selvelli et al 



The inclination, i, can be estimated from our data via Ki - sini, where Vi = 
(2na)/P {q/(q + 1 )) and a is the distance between the two components obtained from 
Kepler's third law. This provides a constraint on the system inclination of i = 1 + 2 
degrees, for any reasonable combination of the component masses. We stress that 
the real uncertainty on the inclination is bound to be larger, because the formal 
error does not include systematic uncertainties associated with, for example, pos- 
sible bias in our Ki measurement and the assumption of tidal truncation in our 
derivation of q. The complete set of spectroscopic system parameters for T Pyx is 
summarised in Table 14.21 
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4.4 Summary and Discussion 



The main result of our study is the spectroscopic determination of T Pyx's orbital 
period, Porb - 1-83 hours. This confirms that the system is a CV below the period 
gap and implies that its current accretion rate is at least 2 orders of magnitude 
higher than that of an ordinary CV at this period. We also find that our spectro- 
scopic orbital period is consistent with the photometric ephemeris found for T Pyx 



I Patterson et al 



19981, an updated version is given in Section l4.2.2|l . This means not 



only that photometric timings can be used as a more precise and convenient tracer 
of the orbital motion, but also that the large period derivative required by the pho- 
tometric ephemeris marks a genuine change in the orbital period of the system. In 
fact, the spectroscopic data are consistent with the photometric period only if the 
period derivative is accounted for. The period derivative obtained here from the 
latest co mbined photometric d ata (P = 6.7 x 10"^°) is slightly lower than that ob- 



tained by 



Patterson et al. 



(1998D from data up tol997_(P^^ lO"^" ). A decrease in 



Schaefer et al. 



( 2010 ) scenario that T 



the rate of period change would be in line with|l 
Pyx's days as a high-M recurrent nova are numbered (at least until its next ordinary 
nova eruption). In any case, the average time-scale for period change found across 
all the photometric data is about 3x10^ yrs. 

We have also used our spectroscopic data to obtain estimates of other key sys- 
tem parameters, most notably the radial velocity semi-amplitude of the WD, - 
17.9 +1.6 kms"\ and the mass ratio, q - 0.20 + 0.03. The latter estimate rests on 
three key assumptions: first, that our determination of is correct, second, that 
our estimate v(Rdisc)sinf is correct, and, third, that the accretion disc around the 
WD in T Pyx is tidally limited. Taken at face value, this relatively high value of the 
mass ratio implies that the donor star in the system is not a brown dwarf. Thus 
T Pyx is not a period bouncer. If we assume that the radius of the secondary is 
1 0% + 1 0% inflated relative to an ordinary main sequence star of the same mass, we 
find that its most likely mass isM2 = 0.14 + 0.03 Mq. 

Overall, the physi cal picture that emer ges from our study is consistent with the 



scenario proposed by lSchaefer et al.l (|2010t) . In particular, they suggest that, prior to 
the 1866 eruption, T Pyx was an ordina ry CV. That eruption then triggered a high- 
M wind-driven phase, as suggested by .Knigge et al.l (|2000|) to account for T Pyx's 
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exceptional luminosity. However, this phase is not quite self-sustaining, so that T 
Pyx is now fading and perhaps not even a RN anymore. In line with this picture, 
we find that the mass ratio and donor mass we derive are not abnormally low for 
a CV at its orbital period. This shows that the present phase of high-M accretion 
cannot have gone on for too long already. The hint of a declining period derivative 
may point in the same direction, but this needs to be confirmed. 

Does all this mean that the phase of extraordinarily high accretion rates T Pyx 
is currently experiencing will have no lasting impact on its evolution? Not neces- 
sarily. In a stationary Roche-lobe-filling system, the orbital period derivative and 
total mass-loss rate from the donor are related via Equation ll.lOl In T Pyx, which 
has an increasing orbital period, we use a mass-radius index a ^ -1/3 (see Sec- 
tion| 1.3.4l) , which is app ropriate for adiabatic mass-loss from a fully convective star 
(e.g. Knigg e et al.ll2000h . We thus expect that M2/M2 - ^orb/^orb/ which suggests a 
typical mass-loss rate from the donor of M2 ~ 5 x 10""^ Mq/jt in its current state. 
(The accretion rate onto the WD can be lower than this, since, in the wind-driven 
scenario, much of this mass escapes in the form of an irradiation-driven outflow 
from the donor.) If every ordinary nova eruption in T Pyx is followed by ~ 100 
yrs of such high-M2 evolution, the total mass loss from the donor in the luminous 
phase is ~ 5 X 1 0"^ Mq between any two such eruptions. This needs to be com- 
pared to the mass lost from the donor during the remaining part of the cycle. If 
this is driven by gravit ational radiat i on, th e mass loss rate from the donor will be 



M2 - 5 X 10"^^ Me/yrlKniggeetal 



The rec urrence time of or dinary nova 



Yaron et al 



(|2005h . so the total 



eruptions for such a system is on the order of 1 0^ yrs[ 
mass lost from the donor during its normal evolution (outside the wind-driven 
phase) is 5 X 1 0"^ Mq. This shows that the long-term secular evolution may be dominated 
by its high-M wind-driven phases, even if the duty cycle of these phases is very low (e.g. 
0.1% for the numbers adopted above). 

It is finally tempting to speculate on the relevance of T Pyx-like evolution for 
ordinary CVs. At first sight, the numbers above suggest that the evolution of a CV 
caught in such a state may be accelerated by about an order of magnitude. This is 
interesting, since it could bear on the long-standing problem that there are fewer 
short-period CVs and period bouncers in current samples than theoretically expec- 
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ted (e.g. 



Patterson et al. 



1998 



Pretorius et al. 



2007a 



Pretorius & Knigge 



2008a 



b) 



Moreover, this channel need not be limited to systems containing high-mass WDs 
that are capable of becoming RNe. After all, it is not the recurrent nova outbursts 
that are of interest from an evolutionary point of view, but simply the existence of 
a prolonged high-M phase in the aftermath of nova eruptions. In CVs containing 
lower-mass WDs, such a phase may still occur, although its evolutionary signific- 
ance could still depend on the WD mass. For example, the inter-outburst time- 
scale is longer for low-mass WDs, and the duration of the high-M phase could also 
scale with WD mass. One obvious objection to this idea is that, observationally, 
most nova eruptions are not followed by centuries- (or even decades-) long high-M 
phases. However, this need not be a serious issue. Most observed novae are long- 
period systems, so if the triggering of a wind-driven phase requires a fully convect- 
ive donor, most novae would not be expected to enter such a phase. It may be rel- 



evant in this context that at least one other short-period nova - GQ 
an ex c eptionally long post -outburst supersoft phase of ~ 1 years 



Mus - exhibited 



( 1993h . 



Greiner et al. 



Oegelman et al 



( 20031) have also suggested that the duration of the supersoft 
phase in novae may scale inversely with orbital period. 

However, there is another important consequence to the idea that the evolu- 
tion of many short-period CVs is accelerated by T Pyx-like high-luminosity phases. 
The orbital period-derivative is positive in the high-M phase, but negative during 
the remaining times of GR-driven evolution. But since Porb/^orb - Mj/Mj in both 
phases, the sign of the secular (long-term-average) period derivative will gener- 
ally correspond to the phase that dominates the secular evolution. So whenever 
wind-driving dramatically accelerates the binary evolution, the direction of period 
evolution will also be reversed. Is this a problem? Perhaps. The rec ent detection of 



the lo ng-sought period spike in the distribution of CV orbital periods 



Gansicke et al 



( 20091) suggests that there is, in fact, a reasonably well-defined minimum peri od for 



CVs. as has long been predicted by the standard model of CV evolution (e.g. 



Kolb 



19931 ). On the other hand, the location of the observed spike (^ 83 mins) is signi- 



ficantly different from the expected one 65 min - 70 min). Is it possible that the 
observed period minimum corresponds to the onset of wind-driving in most CVs? 
T Pyx, with P^j-b - 1 1 min would then have to be an outlier, however, perhaps 
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because of an unusually high WD mass. 

The idea that T Pyx-like phases may significantly affect the evolution of many 
CVs is, of course, highly speculative, and we do not mean to endorse it too strongly. 
However, it highlights the importance of understanding T Pyx: until we know what 
triggered the current high-luminosity state, it will remain difficult to assess the 
broader evolutionary significance of this phase. Note that the apparent uniqueness 
of T Pyx is not a strong argument against such significance. For example, if the duty 
cycle of high-luminosity phases is ~ 0.1%, as suggested by the numbers above, we 
should not expect to catch many CVs in this state. Thus T Pyx could be the tip of 
the proverbial iceberg. 
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Chapter 5 

SDSS J1507+52 

- An eclipsing period bouncer in the Galactic halo - 



The main parts of this chapter will appear in MNRAS as \Uthas et al\ ilOllca) . I would 
like to thank my collaborators, Christian Knigge, Knox S. Long, Joseph Patterson and John 
Thorstensen for guidance and valuable comments. 

SdSS J1507+52 (hereafter J1507) is an eclipsing cataclysmic variable that consists 
of a cool, non-radially pulsating white dwarf and an unusually small sub -stellar 



secon dary. The system was first identified from the Sloan Digital Sky Survey (|Szkody et al 
20051 ) and was quickly recognised to be odd due to its short orbital period of about 
67 minutes, which is well below the minimum period characteristic of normal CVs. 
A few other systems are occasionally found to have periods below the usual min- 
imum period, but, in general, these systems show evidence for abnormally hot 
and bright donor stars in their optical spectra, suggesting that their secondaries are 
nuclear-evolved objects. This is not the case in T1507, whose secondary is not vis- 



ible at all in optical spectroscopy (jSzkody et al 



20051: 



Littlefair et al 



20071). Together 



with its short orbital period, this indicates a system with a faint disc, a low accretion 
rate and a relatively unevolved, low-mass donor. In Figure 15.11 1 present a scaled 



mo del painting of 



Littlefair et al. 



1507. 



( 20071) performed eclipse analysis of J1507 and found a donor 
mass of 0.056 + 0.001 Mq, clearly below the hydrogen-burning limit. However, be- 
cause of its anomalously short orbital period, the sy stem is not consi stent with the 



stanc 



ard mass-period relation for CVs (lKniggel2006l ). In line with this. 



Littlefair et al. 



( 20071) found the donor radius to be smaller than predicted by standard CV evolu- 



tion sequences. They therefore suggested that the donor might be unusually young. 
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Figure 5.1: My impression ofSDSS J1507+52. The model is scaled from u sins: the estim- 
ated b inary separation and the radii of the primary and donor as presented by lLittlefair et al. 



i.e. that J1507 might be a CV that formed recently from a previously detached WD - 
brown dwarf close binary system. A young brown dwarf has a higher density, and 
therefore a smaller radius, since the donor has not yet had a chance to expand in 
response to mass loss. This would explain the short orbital period found in J1507, 
since according to the period-density relation, a higher density and subsequently 
a smaller radius, implies a shorter orbi tal period. The list o f photometrically in- 



ferred system parameters provided by 



Littlefair et al. 



( 2007 ) also includes the ef- 



fective temperature of the white dwarf, which they estimate to be Tgff = 11000 + 500 
K. 

Concurrently, 



Patterson et al 



1 20081 ) showed that the system has an unusually 
high space velocity, similar to the velocities of stars in the Galactic halo. A typical 
star of 1 Mq in the Galactic disc has a space velocity below 50 kms"^, while J1507 
has a velocity of about 167 kms~^. If J1507 is a member of the Galactic halo, its 
donor will be a Population II object with low metallicity. Due to their lower atmo- 
spheric opacity, such objects have significantly smaller radii than their solar metal- 
licity Population I counterparts. Membership of the Galactic halo would therefore 
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also account for J1507's small donor radius and short orbital period. However, Pop- 
ulation II stars represent only about 0.5 % of the stars in the S olar neighbourhood. 



which would make J1507 a rare system. iPatterson et al.l (|2008n estimated a slightly 



higher effective temperature of the WD in J1507, Tgff = 11500 + 700 K. They also 
found the system to exhibit multi-periodic variability on time-scales of minutes, 
which they interpreted as non-radial pulsations originating from the primary white 
dwarf. 

Both hypotheses imply that J1507 is an interesting and important system that 
can shed light on poorly understood aspects of CV evolution. However, since the 
models are very different, it is clearly important to determine which - if either - 
of them is correct. As noted above, eclipse-based estimates for the effective tem- 
perature of the WD suggested Teff ^ 10000 K - 12000 K. At these temperatures, the 
ultraviolet (UV) spectrum of a WD accreting metal-rich material would be expected 



19941) . 



to show strong absorption features due to Fe II and Fe III (c.f. Horne et al. 
In this chapter, I present far- and near-UV spectroscopy of J1507 that was obtained 
with the COS and STIS instruments onboard the Hubble Space Telescope, with the 
aim of measuring the metallicity of the system and establishing its status as either 
a Pop I (disc) or a Pop II (halo) object. 



5.1 Observations AND Reduction 

Far- and near-UV observations of J1507 were performed in February 2010, using 
the HST. The observations were obtained in the so-called time-tag mode. In this 
mode, each photon event is tagged with a corresponding event time, along with 
their wavelength and spatial position on the detector. This allows us to rebin the 
data over any time interval and wavelength region. 

Far-UV observations were carried out with the Cosmic Origins Spectrograph 
(COS) for a total exposure time of about 3 hours, covering 3 eclipses in 3 successive 
HST orbits. The maximum time resolution of the COS TIME-TAG mode is 32 ms. 
The G140L grating was used, which gave a spectral range of 1230 A - 2378 A and a 
spectral resolution of 0.5 A. 

Near-UV observations were obtained by the Space Telescope Imaging Spectro- 
graph (STIS) for a total exposure time of about 1.4 hours, at a maximum time resol- 
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Orbital Phase 



Figure 5.2: The Gaussian smoothed and flux-calibrated light curve obtained by COS/HST, 
constructed from the monochromatic flux averaged over the region 1425 A - 1900 A. Bot- 
tom panel shows the phase-folded light curve in the vicinity of the eclipse. 

ution 125 ^s. The G230L grating was used, resulting in a spectral range of 1650 A - 
3150 A and a spectral resolution of 3.16 A. The data set consists of 3 sub-exposures 
(of length 1600 s, 1918 s and 1708 s). No near-UV data was obtained during the ec- 
lipses. Also, there are gaps in the data caused by interruptions coming from Earth 
occupation. 

All data were reduced and calibrated using the Pyraf package Stsdas, provided 
by STScI (the Space Telescope Science Institute]^ 

5.2 Analysis and Results 

Figure 15.21 shows the Gaussian smoothed and flux-calibrated light curve from the 

far-UV data. This light curve was constructed by dividing the larger exposures into 

one-second bins, and by averaging the fluxes over the wavelength range 1425 A - 

^Stsdas and Pyraf are products of the Space Telescope Science Institute, which is operated by 
AURA for NASA. 
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Figure 5.3: The two top panels show the mean out-of-eclipse spectrum for both the far (left) 
and near (right) UV regions. A model constructed from parameters obtained from the far 
UV-region is plotted in red (see Table \5.1h A model corresponding to values from Littlefair 
et al. (2007), Tf,ff= 11000 K and logg = 8.5, is plotted in blue. Underneath each wavelength 
region is a plot showing the residuals. The dips at (1460, 1530, 1600 and 1670) A coincide 
with regions of bad pixels. The 1-a flux errors are plotted in green. 

1900 A. The light curve shows sharp, square-shaped eclipses, close to zero flux at 
mid-eclipse, and a very flat out-of-eclipse continuum, all strongly suggesting that 
the UV flux is dominated by the accreting WD. The folded light curve construc- 
ted from the combined data for all three HST orbits is presented in the bottom 
panel. All further analysis was carried out on the unsmoothed data, excluding the 
eclipses. Ground-based optical observations were performed at the same time as 
the HST observations, yielding a magnitude for J1507 of g = 18.44 + 0.02. This is 
consistent with SDSS, showing that the system was in quiescence during the HST 
observations. 

The two top panels of Figure 15.31 show the mean out-of eclipse spectrum (in 
black), for both the near- and far-UV. The reduction was performed separately for 
the COS and STIS data, and the mean spectrum for both regions were joined at 
1900 A without overlap for a total range of 1268 A - 3150 A. The COS data were 
rebinned slightly, to 0.2 A/ pix, to increase the S/N and provide sufficient sampling 
of the G140L grating's spectral resolution. Different dispersions were retained for 
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the the far- and near-UV regions (0.2 A/pix and 1.58 A/pix, respectively). The 
overall UV spectrum for J1507 is relatively sparse in spectral lines. We identify CII 
at 1335 A, Sill at 1526.7 A, Hell at 1640.5 A, AlII at 1671 A and Mgll at 2800 
A. The strong emission feature near 1300 A is almost certainly due to geocoronal 
emission in the O I 1304 A line that could not be cleanly removed by the pipeline 
background subtraction. 



5.2.1 Spectral Modelling 

A model grid spanning the four key atmospheric parameters, effective temperat- 
ure (Tgff), surface gravity (logg), metallicity {[Fe/H]) and rotational velocity {v sinz), 
was constructed in order to fit the far- and near-UV spectrum. Atmospheric struc- 
tures were calculated with Tlusty, with the sp ectral synthesis being done with 



SynSPEC§ jHubeny 



1988 



Hubeny & Lanz 



19950 . The grid spanned the region; 



12500 K < Tgff < 20000 K in steps of 500 K, 7.5 < log^ < 9.5 in steps of 0.25, -2.0 
< [Fe/H] < 0.0 in steps of 0.25, and kms"^ < vsinf < 1000 kms"^ in steps of 50 
kms"^. Models at intermediate parameter values were constructed by linear in- 
terpolation on this grid. Pixels with non-zero data quality flags set were excluded 
from our fits, as well as the regions around the spectral lines, OI at 1304 A CII at 
1335 A and the quasi-molecular H feature at ~ 1400 A. 

Before the models could be compared to the observed spectrum, they were con- 
volved with a Gaussian filter to the spectral resolution appropriate to the COS (< 
1900 A) and STIS (> 1900 A) data, and linearly interpolated onto the observational 
wavelength scale. In general, our model fits did represent the main features of the 
data quite well, but the formal was somewhat high. In order to obtain more 
realistic formal parameter errors, we therefore added an intrinsic dispersion term 
to the flux errors in such a way that xl = ^- This dispersion corresponds to about 
1% of the flux at 1900 A. Formal 1-a parameters errors were then estimated in the 
usual way, by identifying the parameter ranges that correspond to x^ - Xmin- ^■ 

In order to explore the systematic uncertainties affecting our fits, we tried fit- 
ting many different wavelength regions, such as the far-UV (COS), near-UV (STIS), 



^Tlusty and Synspec are programs developed for modelling of stellar atmospheres: 
|http:/ /nova.astro.umd.edxi/ 
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log(g) [cin/s^] V siii(i) [km/s] 



Figure 5.4: Contour plots in x'^space showing the correlation between the fit parameters. 
The left panel shows a strong, almost linear correlation between Tg^and logg. The panel to 
the right shows that the parameters [Fe/H] and v sinf are roughly uncorrected. 

regions with a high line concentration and also the whole wavelength range span- 
ning both the far- and near-UV. 

The correlations between fit parameters were investigated by constructing X'^ 
contour plots for each possible parameter pair. The left panel in Figure 15.41 show 
that there is a clear correlation between Tgff and \ogg. The correlation between 
[Fe/H] and vsini (right panel) is not pronounced. 

5.2.2 Temperature and Metallicity 

A fit to the far UV-region 1268 A - 1900 A yields T^k = 14200 + 50 K, \ogg = 8.2 + 
0.04,[fe/H] = -1.2+0.05 and vsini = 180 + 20 kms"^. These fit parameters provide 
a good match also for the near-UV spectrum out to 2800 A. We adopt these para- 
meters as our best-fit parameters. The quoted errors are formal errors from the X'^ 
fitting, shown in Figure 1531 Table ISJl presents the best parameters along with their 
uncertainties. These uncertainties are larger than the formal errors, because they 
include our best estimates of the systematic uncertainties associated with different 
choices of fitting windows. However, there are several reasons for adopting the 
parameters inferred from the fit to the COS far-UV data set as our best-bet estim- 
ates. First, this region includes almost all line features, and is also most sensitive 
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Figure 5.5: x^-fiiii-^g of the parameters to the COS far-UV spectral region. 

to changes in Tgff and log g. Second, the far-UV region is more likely to represent 
pure light from the WD, while towards redder wavelengths, the spectrum could be 
affected by the disc and the bright spot. Third, fitting only the COS data removes 
any possibility that a mis-match in the flux calibration of COS and STIS could affect 
our results. 

The two top panels in Figure 15.31 show the best-fit model in red, along with 
the mean out-of-eclipse spectrum for J1507. The corresponding residuals between 
model and data are plotted underneath each of the COS and STIS spectral regions, 
together with the 1-a error range marked in green. The dips in the residuals at 
1460 A , 1530 A , 1600 A and 1670 A are caused by bad pixels in the data, and these 



regions were exc 
meters found by 



uded from our fits. A s a comparison, a model based on the para- 



Littlefairetal 



(2007D, Teff = 11000 K and logg = 8.5, is plotted in 
blue. Clearly, this temperature is far too low to match the data. 

The top panel in Figure \57\ shows the mean spectrum together with models at 
three different values of [Fe/H]: (green), -1.2 (red) and -2 (blue). It is clear that a 
model of solar metallicity (green line) does not match the data and overestimates 
both the number and strength of visible absorption lines in the spectrum. The left 
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2; 




Figure 5.6: Figure from \Carney et al\ 09961) . showing the metallicity distribution for 
Galactic halo stars. 

panel shows a zoom of the region around the Si II line (the coloured lines corres- 
pond to the same values of the metallicity as given for the top panel). To the right is 
a plot showing the dependence of the model on v sin i, where the coloured lines in 
green, red and blue represent rotational velocities of kms~\ 180 kms"^ and 500 
kms"\ respectively. In all panels, the red lines show the best-fit parameters. The 
metallicity and rotational velocity affect both the line depths and widths. There- 
fore, as a check, we fixed Tgff and logg to our best values and made model fits to 
smaller wavelength regions specifically chosen because they contain strong and/ or 
many spectral lines, with vsinz and [Fe/H] as free parameters. We find consistent 
values of vsinf and [Fe/H], independent of the wavelength region we performed 
the fits on. 

The metallicity found for J1507 can be compared with the metalli c ity di stribu- 



Patterson et al 



120081) calcu- 



tion for Galactic halo stars at the same space velocity, 
lated a space velocity for J1507 of 167 km s"^ , which can be broken into th e Galactic 
veloci ty components U, V and W, where ^/U^ + -139 km s"^ for J1507. 



Carney et al 



( 1996h presented a study of the metallicity distribution for the Galactic halo stars. 



They found that for stars with velocities in the range 100 < Vt/^ + < 140, the 
metallicity distribution peaks around -0.7, but has a lon g tail towards metallicit- 



ies as low as -1.5 (see Figure |5^ from ICarney et al 



19961 ). This is consistent with 



our metallicity estimate of -1.2. We therefore conclude that J1507 is most likely a 
Galactic halo CV. 
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Figure 5.7: The top panel shows the mean out-of-eclipse spectrum for the far UV region, 
where models at three different [Fe/H], (green), -1.2 (red) and -2 (blue), are plotted on 
top. The bottom left panel shows a zoom of the region around the Si U, where the coloured 
lines correspond to the same values of the metallicity as in the top panel. In the bottom right 
panel, the colours correspond to vsini at kms'^ (green), 180 kms~^ (red), 500 kms~^ 
(blue). 
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5.2.3 Possible Extinction Effects 

Up to now, we have assumed that the observed far-UV spectrum of SDSS J1507 
is unaffected by interstellar extinction. Based on the absence of the well-known 
2175 A absorption feature in the data, we can set E(B-V) < 0.05 as a fairly conser- 
vative limit on the amount of reddening that may be present. In order to test the 
effect of extinction at this level on our conclusions, we carried out additional model 
fits after dereddening the data by E(B-V) = 0.05. We find that none of our inferred 
parameters would change significantly if extinction at this level were present. 



5.2.4 Non-Radial Pulsations 



The high effective temperature found for J1507 (14200 + 500 K) means that the sys- 
tem is pushed out of the instability strip for non-accreting WDs (12000 + 1200 K; 
see Figure 3 of Gianninas et al.ll2006 ). However, observations show that the temper- 
atures found for pulsating WDs in CVs, cannot be defined by this narrow instabil- 
ity stri p, and their temper atures are found in a wider range, often towards higher 



values (ISzkody et al. 



2010). 



Nilsson et al 



( 20061) in 2005 , 



Patterson et al 



120081) 



Time-resolved photometry of J1507 was obtained by 
but no periodic signals were found. However, also in 2005, 
found signals in the light curve for J1507 at 76 c d^^, 130 c d~^ and 175 c d~^ (+ 0.5 
cd~^), which they identify as non-radial pulsations from the primary WD. These 
signals were found to have moved to longer periods when measured a year later, 
and was then found at 75.3 cd"\ 127.8 cd"i and 171.3 cd"^ (+0.6 cd'^). 

I obtained observations of J1507 at the MDM telescopes during one night in 
May 2010 and recovered two of these signals, now at 90 cd"^ and 130 cd"^ (+1 
cd~^). The amplitudes of these pulsations are expected to vary with wavelength, 
and in the COS far-UV data, a strong peak is, in fact, found 85.3 + 0.5 d"^ (See 
Figure |5^ . All eclipses were removed from both the optical and UV light curves 
before constructing the Lomb-Scargle periodograms. A complete search for the 
UV counterpart of the optical WD pulsations will be presented in Uthas et al. (in 
preparation). However, it is already clear that both the new optical and the UV 
data support the classification of J1507 as a GW Lib pulsator - the first such object 
in the halo and thus with a low-metallicity WD envelope. 
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Figure 5.8: Lomb-scargle periodogram constructed from the monocromatic integrated 



5.2.5 Distance Estimates 

The theoretical relationship between observed flux (f ) and the Eddington flux (H) 
provided by the Synspec models can be used to estimate the distance {d) towards 
J1507. More specifically, F - AnR^^H/ d^, where jRwd is the WD radius. This means 
that the normalisation factor needed to optimally match a model spectrum to the 
data is a direct measure of R^^/d'^. WDs also follow the well-known mass-radius 
relation (with only a weak temperature sensitivity), so the surface gravity of the 
model, g - GM^^/R^^, uniquely fixes Kwd- Thus for a WD model with given Tgff 
and logg, the normalisation factor required to fit the data is a unique measure of 
distance. 

In practice, the distance was estimated by first fitting a linear function to the 

relationship between M^^ and logg in Pierre Bergeron's WD cooling model at 

Tgff = 14500 K. Next, this function was used to estimate M^^/ srid hence Kwd/ 

for a given logg. The resulting estimates of the WD mass and radius are found 

■^Cooling models by Pierre Bergeror\ are found at: 
|http:/ /www.astro.umontreal.ca/ ~bergeron/CooiingModels/| 



fluxes from thefar-UV region from 1268 A - 1900 A (eclipses were removed). 



120 



5.3. DISCUSSION AND SUMMARY 



Teffi 14200 + 500 (50) K 

logg: 8.2 + 0.3 (0.04) 

[Fe/H]: -1.2 ± 0.2 (0.05) 

vsini: 180 + 20 (20) km s'^ 



Table 5.1: Best-fit parameters for J1507 obtained from model fitting to thefar-UV region. 
The errors are defined as the whole range for which good fits are found, irrespective of 
wavelength region. Formal 1-a errors are given in parenthesis. 



to be; = 0-75 + 0.15 Mq and R^^ = 0.011 + 0.002 Rg. These values are not 
very precise because the surface gravity is only poorly co nstrained by our spe ctro- 



scopy. In any case, ou r numbers are consistent with both 



Patterson et al. 



Littlefair et al 



( 2007 ) and 



(|2008l) . The normalisation constant of the model combined with 
-Rwd then finally yield the corresponding distance estimate. The distance towards 
J1507 is thus found to be d = 250 + 50 pc (the effect of reddeni ng is allowed for in 



the quoted error). This is consistent wi th 



Patterson et al 



distance given by 



Littlefair et al. 



120081) but larger than the 



(|2007l) . since the logg^ presented here (logg = 8.2) 



is lower than the logg implied by their estimates of Mi and Tgff (logg = 8.5). 



5.3 Discussion and Summary 

In this chapter, I have presented HST observations in the UV spectral-range of the 
cataclysmic variable J1507 that was obtained with the aim of measuring its metalli- 
city to determine whether or not the system is a member of the Galactic halo. 

By comparing the observed spectrum to synthetic spectra described by the four 
parameters, Tgff, logg, [Fe/H] and vsinz, a best fit is found at an effective temper- 
ature of 14200 + 500 K. This value of Tgff is based on the assumption that the WD 
is the only component contributing to the UV flux in J1507. If a hot, optically thick 
boundary layer would be present, it would bias the estimate of Tgff. However, the 
single-temperature WD models presented here seem to provide a good fit to the 
data, and the boundary layer in a low M-system (such as J1507) is expected to be 
optically thin and thus not contribute significantly to the UV flux. 
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The best model-fit give a Tgff that is much higher than previous estimates, which 
implies that the accretion rate is also higher than previously thought. More spe- 
cifically, since M oc T"*, the 30% incre ase in the estimated Tpff c orresponds to an 
increase in M by almost a factor of 3 (jTownsley & BildstenI 120031) . This may have 
significant implications for the evolution of this system, including the question of 



whether gravitationa 



Townsley & Gansicke 



radiation alone is suffic ient to drive this accretion rate (see 



2009 



Knigge et al.ll201lh . 



At the high Tgff presented above, Fe II and III are no longer dominant contribut- 
ors to the atmospheric UV opacity, making it more difficult than expected to meas- 
ure the metallicity of the system. Nevertheless, model fits to the data clearly favour 
a significantly sub-solar metallicity, [Fe/H] = -1.2 + 0.2, comparable to the typical 
metallicity found for halo stars at the same high space velocity as J1507. This im- 
plies that J1507 is most likely is a halo CV. 

This result raises the question if there are any other confirmed Galactic halo 
CVs. Observationally, they are expected to be diffi cult to ident i fy, sin ce they are 



faint and typically will be located at large distances. IStehle et al.l (|l997h performed 
theoretical population synthesis on Population II CVs, and found that halo CVs 
would have higher accretion rates and therefore shorter evolutionary timescales 
compared to CVs with a solar chemical composition. As a result, they would pass 
through the period gap more quickly and also have a shorter minimum period. 
This results in an even larger percentage of short-period systems for halo CVs than 
is expected for the Galactic disc population. Also, due to the lowe r metallicity. 



dwar f novae would be expected to exhibit higher outburst amplitudes dCannizzo et al 
1988|). So far, only a few surveys have directly ai r ned to identify halo CVs (IHowell & Szkody 



1990 : 



Mukai et al 



19901: 



Szkodv & Howell 



19921 ). 



Howell & Szkodyl (|1990D collected 



a sample of candidate halo CVs that consisted of 84 high-Galactic latitude objects. 
They were chosen on the b asis that they had latitudes > 30° and were very faint. 



indicating large distances. 



Howell &: Szkody 



( 19901) found a slightly larger frac- 



tion of short-period systems among their candidates. They also find higher out 
burst amplitudes. All this might indicate that some of these candida tes are in fact 



laloCVs, although later surveys using infrared photometry dat a, e.g. 



( 1993) , show that the method used by 



Howell et al. 



Howell &: Szkodvl (1990) overestimates the 
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distances, which would imply that their CV sample is probably situated closer to 
the Galactic plane than they estimated. In order to confirm a membership in the 
Galactic halo, one needs to measure the space velocities and/ or metallicities of any 
candidates. In the absence of such measurements for any other candidate halo sys- 
ter ns, T1507 is the most li kely halo CV known to date. 



Patterson et al 



(|2008|) found periodic signals in J1507, which they identified as 



non-radial pulsations originating from the primary WD. N on-accreting CVs show 



pulsations within the instability strip, 10900 K - 12200 K jGianninas et al 



20061). 



The effective temperature found for J1507 is well above the instability strip for non- 
accreting pulsating WDs. However, it is already becoming clear that non-radially 
pulsating WDs in CV s are found across a wider temperature range, often towards 
higher temperatures (jSzkody et al.ll2010h . The fact that non-radial pulsations are 
present in J1507 is certainly interesting and potentially important, since a higher 
metallicity in the outer envelope of th e accretors rnight be able to push the instabil- 



ity strip toward higher temperatures (jArras et al 



20061). In this context, it is inter- 



esting to note that J1507 has a very low metallicity, but is nevertheless managing to 
pulsate at an effective temperature of above 14000 K. 
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Chapter 6 

Summary & Discussion 



I 



have presented analysis of four cataclysmic variables with orbita l periods close 



to the minimum period for CVs (« 83 minutes: iGansicke et alJl2009l ). All four sys- 
tems are interesting not just in their own right, but also from the point of view of 
CV evolution. In order to understand the importance of these systems, and to high- 
light what impact they may have on our current understanding of CV evolution, 
I have also given an overview of the theoretical and observational state-of-the-art 
in CV research in general. In this chapter, I will provide a brief summary and dis- 
cussion of the main results. The discussion will, in particular, focus on how the 
study of individual systems - especially ones that are seemingly rare or unusual - 
can provide unique and vital insights that cannot easily be gained by other means. 
In addition, I will also briefly point out some promising directions for future work 
suggested by the research described in this thesis. 



6.1 Summary 

6.1.1 sdss j1457+51 and bw sculptoris 

In Chapter |3l I introduced two systems that appear to follow the standard evolu- 
tionary track for CVs; SDSS J1457+51 and BW Sculptoris. Bo th svstems are fain t 



2003), 



and show prominent white-dwarf components in their spectra (|Szkody et al 
indicating low accretion rates. I showed that both stars exhibit non-coherent pulsa- 
tions in their light curves and are likely new members of the GW Lib class of 
variable stars, in which the periodic and non-commensurate signals are widely 
believed to arise from non-radial pulsations in the underlying white dwarf. The 
power spectra of both stars show complex signals with primary periods near 10 
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and 20 minutes. These signals change in frequency by a few percent on a timescale 
of weeks or less and probably contain an internal fine structure unresolved by the 
observations. Double-humped waves were also found, marking the underlying or- 
bital periods, near 78 minutes for both stars. This is right at the observed minimum 
orbital period for CVs. Their short orbital periods together with their low accre- 
tion rates indicate that these stars are likely to be very old CVs. In addition, BW 
Sculptoris showed a transient, but powerful signal with a period near 87 minutes 
- a quiescent superhump. The ~ 11% excess over the orbital period is difficult to 
understand and may arise from an eccentric instability near the 2:1 resonance in 
the accretion disk. 



6.1.2 TPYXIDIS 

T Pyxidis is a luminous recurrent nova that accretes at a much higher rate than 
is expected for its photometrically determined orbital period of about 1.8 hours. 
In Chapter m I presented the first spectroscopic orbital period, P = 1.8295 hours, 
based on time-resolved optical spectroscopy obtained at the VLT and the Magellan 
telescopes, confirming T Pyxidis as a short orbital system. An upper limit of the 
velocity semi-amplitude of the white dwarf, Kl = 17.9 + 1.6kms"^ was derived, 
as well as an estimate of the mass ratio, = 0.20 + 0.03. If the mass of the donor 
star is estimated using the period-density relation and theoretical main-sequence 
mass-radius relation for a slightly inflated donor star, one obtains M2 = 0.14 + 0.03 
Mq. This implies a mass of the primary white dwarf of My - .7 + 0.2 Mr^,. If th e 



white-dwarf mass is > 1 Mq, as classical nova models imply (jYaron et al.l 120051) , 
the donor mass must be even higher. Therefore, T Pyx has, most likely, not yet 
evolved beyond the period minimum for cataclysmic variables. I constrained the 
system inclination to be z « 1 degrees, confirming the expectation that T Pyx is a 
low-inclination system. 



6.1.3 SDSSJ1507+52 

SDSS J1507+52 is an eclipsing cataclysmic variable consisting of a cool, non-radially 
pulsating white dwarf and an unusually small sub-stellar secondary. The system 
has a high space velocity and a very short orbital period of about 67 minutes 
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( Szkody et al.l. l2005l) , well below the usual minimum period for CVs. In order to 
explain the existence of this peculiar system, two theories have been proposed. One 
suggests that SPSS T1507+52 was fo rmed from a young, detached white-dwarf /brown- 
dwarf binary (|Littlefair et al.ll2007n . The other theory propose s that the system is a 



member of the Galactic halo population (jPatterson et al. 



2008) 



In Chapter |5l I presented ultraviolet spectroscopy of SDSS J1507+52 obtained 
with the Hubble Space Telescope with the aim of distinguishing between these two 
theories. I found that the UV flux of the system is dominated by its accreting white- 
dwarf. Fits to model stellar atmospheres yielded physical parameter estimates of 
Tgff^ 14200 + 500 K, logg = 8.2+ 0.3 cms"^ v sini = 180 + 20kms"i and [Fe/H] = 
-1.2 + 0.2. These fits suggest a distance towards SDSS J1507+52 of = 250 + 50 
pc. The quoted uncertainties include systematic errors associated with the adopted 
fitting windows and interstellar reddening. 

The strongly sub-solar metallicity found for SDSS J1507+52 is consistent with 
that of halo stars at the same space velocity. I therefore concluded that it is likely 
that SDSS J1507+52 is a member of the Galactic halo. In addition, I found peri- 
odic signals in the U V spectrum of SPSS T1 507+52 with frequencies close to those 
already identified by IPatterson et al.l (|2008|) in the optical. These signals could be 



interpreted as non-radial pulsations originating from the primary white dwarf. 



6.2 Discussion 



6.2.1 Importance of Individual CVs 

With a growing population of newly discovered CVs, it has become the norm to fo- 
cus on statistical analys es of large CV samp 



tion as binary stars (e .g. 



2009 



Patterson 



2011 



Littlefair et al 



2008 



es, in efforts to understand their evolu- 



Townslev &: Gansickel200' 



Gansicke et al. 



Knigge et al.ll201lh . Such studies are important and yield a 
fundamental understanding of the common evolutionary behaviour for CVs. How- 
ever, to quantitatively test the proposed theoretical models, careful studies of well- 
chosen individual systems, in addition to statistical studies, can also provide us 
with important and unique insights. Thus statistical studies and studies of indi- 
vidual systems are complementary. 
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In particular, the study of systems passing through key evolutionary phases, 
such as the minimum period, is of great importance for the general understanding 
of CV evolution. All four systems presented in this thesis have orbital periods close 
to the minimum period, indicating that they are in exactly this transitional phase. 
Also, studies of CVs that are able to provide us with precise system parameters, 
such as the well-defined properties found from eclipsing systems, are highly valu- 
able and can provide more stringent tests of theory (for instance SDSS J1507+52). 

Another example of individual CVs that deserve detailed study are systems 
displaying non-radial pulsations (NRPs) from the underlying WD (such as SDSS 
J1457+51 and BW Sculptoris). The study of NRPs in GW Lib stars may provide im- 
portant information on how the process of accretion is affecting the evolution of the 
WDs, compared to isolated ones. Although our current theoretical understanding 
of the mechanism for exciting non-radial pulsations in CVs is limited, modelling 
of these signals might ultimately offer the most precise ways of constraining key 
parameters (see Section l6.2.2D . 

More generally, the study of apparently unusual individual systems will always 



rem ain important, since al 



(e-g- 



Pretorius et al 



2007b; 



observational samples ar e afflicted by selection biases 



Pretorius & Knigge 



2008a 



b|). Given that even the best 



available samples are usually imperfect tracers of the underlying population, un- 
usual systems that do not appear to follow the common behaviour (outliers) could 
be more important than their apparent rarity would suggest. 

One example is T Pyxidis. With its unusually high luminosity and high ac- 
cretion rate compared to other CVs at its period, the system appears to be a clear 
outlier. However, epochs of enhanced mass transfer (like its present state) may 
accelerate or even dominate the overall evolution of the system, even if they are 
relatively short-lived. If most or all systems were to go through similar phases, 
this could be considered as an evolutionary short-cut that may be relevant to the 
evolution of cataclysmic variables more generally. Note that even if all CVs did go 
through such phases, and even if these phases were important to their overall evol- 
ution, due to the short evolutionary time scales associated with these phases, we 
might only ever see a handful of systems like T Pyx observationally. In this context, 
a considerable fraction of CVs might at some point in their evolution appear as out- 
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Hers for some short-lived part of their lives. As a cor\sequence, the systems we see 
as outliers today may very well nevertheless represent a common and important 
phase of CV evolution. 

SDSS J1507+52 is another example of a system that clearly appears to be an 
outlier. Only if one assumes that this system belongs to the Galactic halo popu- 
lation can its main physical properties be explained. Compared to CVs found in 
the Galactic disc, halo CVs are expected to have higher accretion rates, resulting 
in a shorter evolu tionary time scale, and therefore also a shorter minimum period 



(Stehle et al 



1997n . SDSS J1507+52 is actually interesting in two respects: first, it 
belongs to the rare sample of CVs that show evidence of harbouring sub-dwarf 
donors, and second, it is probably the first confirmed CV in the Galactic halo. This 
provides a unique opportunity to study the process of almost pure hydrogen accre- 
tion in a CV, and also give us the possibility to study NRPs from the primary WD 
without the complication of metallicity. 

6.2.2 Future Work 

Non-coherent signals in the light curves of the GW Lib stars, are primarily thought 
to arise from NRPs in the underlying WD, and have been detected for only a decen- 
nium. However, to exclude other possible origins of these signals, these group of 
CVs need to be studied in more detail. Future theoretical modelling of NRPs, pre- 
dicts more precise constrains of system parameters in CVs. For instance, the WD 
mass and WD spin periods could be measured in a more precise manner, as well 
as parameters that we are unable to measure today, such as the accreted envelope 
mass, wh ich would provid e us with a direct estimate of the time since the last nova 



eruption (jArras et al.ll2006h . 

In particular, the study of NRPs in SDSS J1507+52 is interesting since it provides 
an unique opportunity to study the pulsation pattern of an accreting WD in a low- 
metallicity system. Pulsating WDs in CVs are found in a wider temperature range 
than isolated WDs displaying NRPs (see Figure [L8] and [L9l l , which previously has 
been suggested could be a result of the hig her metallicity contents in the outer 



envelope of accreting WDs (jArras et al.ll2006l ). However, SDSS J1507+52 is a low 



metallicity system, but its WD is still found to be pulsating at an effective temper- 
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ature of > 14000 K. This implies that we might have to consider alternative explan- 
ations of the existents of pulsating and accreting WDs at temperatures above the 
instability strip. A detailed study of how the pulsations in SDSS J1507+52 vary 
with wavelength (in the UV region) will be presented by Uthas et al. (in prepara- 
tion). 



As mentioned above in Section 16.2.11 (see also Section I4.4|l , phases of enhanced 
mass-transfer (such as seen in T Pyxidis) may represent an evolutionary stage that 
is not yet included in any evolutionary model. Therefore, theoretical modelling 
of how phases of enhanced mass-transfer would affect the overall evolution of 
CVs, would be highly interesting and potentially important. For instance, popu- 
lation synthesis could indicate what fraction of systems we expect to undergo such 
phases. Also, the effect of irradiation of the donor star by the WD and /or accre- 
tion disc, is not well understood and needs to be studied in more detail. Such 
irradiation effects might be the cause of these prolonged phases of enhanced mass- 
transfer in the aftermath of eruptions, and may not only be important for systems 
with high-mass WDs, such as RNe, but also in systems with lower mass WDs (see 
Section l44ll. 

In addition to the spectroscopic data of T Pyxidis presented in ChapterlH phase- 
resolved spectroscopy of the knots in the shell surroun ding the system, was also 



obtained with the VLT in 2004/2005. 



Schaefer et al. 



(|2010h presented a detailed pho- 



tometric study of the nova shell in T Pyxidis, and concluded that the system most 
likely, was an ordinary CV before the event of the classical nova eruption in 1866, 
and that this event triggered the recurrent nova phase. I intend to examine the 
spectral composition of the nova shell in T Pyxidis, with the aim to study how 
these eruptions have affected the sy stem. A spectro scopic study of the nova shell 
in T Pyxidis was also performed by IWilliamsl (|l982n . However my initial analysis 



of the VLT spectrosc opic data c 



shell as presented by 



WiUiamsl ( 1982 ) 



oes n ot reveal the same spectral composition in the 
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Appendix A 



Formula Derivations 



A.l Mass function 

In this section, I will present brief derivations of the mass functions for both the 
donor star and the primary WD. We start by expressing the distance (fli) between 
the WD and the centre of mass of the binary system 

M1+M2 l+q 

where q = M2/M1 and a is the binary separation. Assuming circular orbits, the true 
orbital velocity of the primary Ui is 

27rfli 2na q . * 

vi = = —, (A.2) 

^ P P l+q ^ ' 

where P is the orbital period. Newton's generalisation of Kepler's third law yields 

(A.3) 



G(Mi+M2) 

If we combine Equations lA.2l and |A.3[ and assume that the observed radial velocity 
of the WD is Ki = sin i, we can eliminate a to obtain 

K^P , 

^— = sin^- 2 (A.4) 

InG (Mi+M2)2 ^ ' 

We can now express the mass function of the donor as 

which provides a lower limit on the donor mass, given only the orbital period of 
the system and an estimate of the radial velocity of the WD [Ki). If we are able to 
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measure the radial velocity of the donor star {K2), a lower limit for the WD mass 
can be estimated. The mass function of the WD is 



^ sin^'zMi where (7 = — ^ = — ^. (A.6) 

2nG ^\{l+q)j ^ Ml K2 



A. 2 Mass Ratio 

In this section, I present a derivation of a relationship that permits a spectroscopic 
estimate of the mass ratio q. In high M systems (such as T Pyxidis presented in 
Chapter m, we can estimate the di sc radius i^Hkr by assuming that the disc extends 



all the way out to the tidal radius (jWarner 



1995a ) 



_ 0.60 a 

-'"^disc ~ ^"^max — TV" T- (A. 7) 

[l+q) 

SmakI 1(19811) investigated the relationship between the true projected velocity at 
the outer disc radius (vr^^ sin i) and the measured spectral peak-to-peak separation 
(A^peak-to-peak)/ and found that VR^i^^sini = AVpeak-to-peak/j^/ where m = 1.05 + 0.05. 
The velocity at the outer disc radius, can then be estimated by 

AVpeak-to-peak = 1.05sini J + 0.05. (A.8) 

If we rewrite Equation I A. 71 by inserting a from Equation IA.3I we find, after some 
algebra. 



i?disc = 0.6(GMi)i/3|^__^| . (A.9) 



2n(l + q) 

Now, combining Equation lA.SI and |A.9l yields 



2/3 



1 / GMtc \^^^ 

AVpeak-to-peak = (1.05-^==)sinfM|/3(l+q)l/3|^_J +0.05. (A.IO) 

Finally, dividing Equation lA.lOl by the expression for Ki from Equation IA.4[ the 
ratio AVpeak-to-peak/^i can be expressed as a function of only q. This expression can 
be used to estimate q, given only observational estimates of Ki and of the peak-to- 
peak separation of disc-formed emission line, and we have 

AVpeak-to-peak ^ ^ 355 (iH \ ^ q 05 (^.11) 
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Appendix B 

SPERIOD 

- A Program for Period Analysis - 



For quick and easy analysis of the period contents in a light curve, I have de- 
veloped routines in the Python programming languagJll In these routines, the light 
curve is normalised and smoothed (I will explain the advantages of this below), 
frequencies with power excess are automatically identified (from the unsmoothed 
light curve), and the errors on these frequencies are calculated. Also, the window 
function is plotted for any desired frequency, and the light curve is folded on any 
period. There are two main programs, Speriod and Speriod-All, which are 
used for analysis of light curves from single nights and light curves from multiple 
nights, respectively. For each of the two programs, there is a suite of routines, 
which the user can use interactively (by answering simple questions that appear 
on the screen). All output data files, log files and figures are saved in folders cre- 
ated by the program. I used these programs for the analysis of SDSS J1457+51 and 
SDSS J1507+52 (presented in Section Hand©. 

Below follows a description of the subroutines in Speriod and Speriod-All. 
Both programs require that PYTHON is installed together with the python packages 
SCIPY and NUMPY. 



B.l Analysis of Single-Night Light Curves 
B.l.l Part 1 

The input file is a light curve with 3 columns; time, magnitudes or fluxes and their 
errors. In the first part of the program, all columns in the input file are read and 
^http:/ /www.python.org/ 
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the mean magnitude or flux is calculated and subtracted from the light curve to 
construct a normalised light curve (Figure B.l(a) and B.l(b)| . All examples in this 
section present data from one observing night of the target SDSS J1457+51. 



B.1.2 Part 2 

In part 2, the user is asked to specify the upper and lower limit for removing out- 



liers. The normalised light curve is plotted to screen (Figure B.l(c) I. Next, a Gaus 



sian smoothing is performed (Figure B.l(d) I, by passing a window along the light 
curve and assigning the weighted mean (according to a Gaussian distribution) to 
the bin in the centre of the window. Equal numbers of points are lost at the start 
and the end of the data set. As demonstrated in Figure B.l(d)[ a smoothing of the 
light curve can help to reveal periodic signals with periods long compared to the 
smoothing window but short compared to the total length of the data set. How- 
ever, this smoothed light curve is only for visually inspecting the data, and the 
unsmoothed, normalised light curve is used in the following steps. In addition to 



smoothing the data, a rebinned light curve (Figure B. 1(e) I is calculated by taking 
the mean of every consecutive, independent set of 3 data points. This effectively 
changes the sampling rate of the data. 



B.1.3 PART 3 



In part 3, a Lomb-Scargle periodogram (|Lomb 



Lomb 


1976; 


Scargle 


1982) 



from the normalised, unsmoothed light curve, with outliers removed. The pro- 
gram will ask the user to specify the frequency range and the oversampling to be 
used. A larger oversampling causes the program to search for power excess in finer 
frequency steps. The output power spectrum is printed to screen (Figure B.l(f)| |. 



B.1.4 PART 4 

Next, frequencies which have excess power (hereafter peaks) are automatically 
found by looping over all values in the Lomb-Scargle data file, to find the highest 
amplitudes. The user is asked to specify a lower search limit in units of power: 
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m 65 303 7 303 75 303.8 303.85 303.9 303.95 304 304.0. 
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(a) Raw light curve. 
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(b) Normalised light curve (mean is removed). 
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(c) Outliers are removed. 



-0.05 



+4.+ij. "*"-»- 



sd and rsbinned Ighfcurvf 



+ + 



2 
I 

fi -0.15 
-0.2 
-0.25 ■ 

"^3^3 65 303 7 303 75 303.8 303.85 303.9 303.95 304 304.0. 
HID 

(e) Rebinned light curve. 



(d) Gaussian smoothed light curve. 
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(f) Lomb-Scargle periodogram. 



Figure B.l: The plots above (showing data from a single night of observations of SDSS 
J1457+51) are produced in the initial steps of the analysis program, Speriod. The normal- 
ised light curve where outliers have been removed, Fi%ure \Kl(c)\ is used for producing the 



Lomb-Scargle periodogram in Figure B.l(f) 
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Figure B.2: "Exam-pie of an output from Speriod, showing the signals of highest amplitude 



found in the power spectrum in Figure B.l(f) 



Please specify lower search limit for finding peaks in the power spec- 
trum : 
»15 



All peaks above this value will be registered. A Gaussian is then fitted to each 
peak. However, this is only to construct a list containing all rough peak frequen- 
cies of interest. This list is later used to limit the frequency range when performing 
Monte Carlo simulations to find the exact frequencies of the peaks and their corres- 
ponding errors. Both the peaks and their Gaussian fits are plotted to screen together 
with the preliminary frequencies (see Figure IB^ . 



B.1.5 Part 5 

In part 5, the errors on the original magnitude and flux values are rescaled. This is 
done by fitting multiple sine waves to the normalised and unsmoothed light curve 
with outliers removed, and scaling the errors until the reduced ~ ^ (the fit is a 
composite of the three highest frequencies found in Part 4). This process is illus- 
trated in Figure IB. 31 This plot is also saved as an output file. 

Please wait - reseating the errors... 

The reduced chi2 is 0.989. Errors are rescaled to 87.0% of original er- 
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Figure B.3: The data light-curve with a fake light curve constructed from the highest fre- 
quencies found in the Lomb-Scargle periodogram, plotted on top. The least-square method 
is used to rescale the errors on the data light-curve. 

rors . 

The now normalised, unsmoothed light curve, with outliers removed and with 
rescaled errors, is saved in an output file to be used in further analysis. 

B.1.6 Part 6 

Next, Monte Carlo simulations are performed with the aim of finding the exact 
peak frequencies and their errors. In this method, the peak errors are found by 
randomly re-distributing the points in the light curve within their errors a repeated 
number of times, and constructing a Lomb-Scargle periodogram each time. The 1-a 
error is then obtained from the output distribution of the peaks found in the peri- 
odograms, both statistically, and by fitting a Gaussian. To speed up the program, 
simulations are only performed in a narrow range around each peak. The frequen- 
cies, amplitudes and errors for every peak are printed to screen and the Gaussian 
fits are plotted (Figure EH. 
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Figure B.4: Monte Carlo simulations of the peaks found in Figure lB^ 

B.1.7 Part 7 

Next, the user will be asked whether or not a spectral cleaning of the Lomb-Scargle 
periodogram should be performed: 

Do you want to do a spectral cleaning? Answer yes or no. (This will 
take less than a minute) 
» yes 

How far down in the power noise do you want to continue cleaning? (Give 

a lower power limit) 

»3 

Cleaning in progress... 

Cleaning 10% from peak with frequency 7.209 c/d 

If the user answers yes, the Lomb-Scargle diagram is imported and the maximum 
frequency is found. A fake light curve, a sine wave at this frequency, is produced 
and its corresponding window function is plotted (as a Lomb-Scargle periodo- 
gram). The peak amplitude of the window function is scaled so that it matches 
the original peak in the Lomb-Scargle periodogram. Next, the window function 
multiplied with a gain factor (10% of the peak amplitude) is subtracted from the 
original Lomb-Scargle periodogram. The same cleaning process is then repeated 



138 




Frequency (c/d) Frequency (cJd) 



(a) (b) 

Figure B.5: Figure [R5(fl)| shows the Lomb-Scargle periodogram of SDSS1457+51, con- 
structed from multiple nights. Figure \B.5(h)\ show the cleaned periodogram of the same 
data set. 

again, first finding the highest peak, now from the partially cleaned Lomb-Scargle 
periodogram, then a fake light curve is constructed together with its corresponding 
power spectrum, and so it continues. This will continue until the peaks are below 
the user-specified lower level. Left are now only the residuals with all the peaks 
removed. Finally, the peaks are added back to the residual periodogram, construct- 
ing the cleaned periodogram. The shape of the peaks in the cleaned periodogram 
corresponds to a virtual peak of the same width as the most narrow peak found 
in the original periodogram, and with an amplitude matching the total calculated 
area of what has been removed. It is important to realise that the cleaned spec- 
trum should not be used to finding frequencies, errors or amplitudes, and is only 
meant as a tool for better visual inspection of the frequency contents of the data. A 
cleaned power-spectrum can still be of much use, especially when inspecting data 
from multiple nights, where the frequency content is more complex. Figure IB.5I 
shows the multiple-night power spectrum of SDSS1457+51, before and after clean- 
ing. 

B.1.8 Part 8 

In part 8, the user is asked whether or not the window function should be cal- 
culated. If the user answers yes, a list of all frequencies found in previous steps 
appears on screen. After choosing a peak from the list, a fake light curve at the 
corresponding frequency and with the same sampling is produced, and a Lomb- 
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Figure B.6: The folded light curve ofSDSS 1457+51 from one single night. It is folded on 
the frequency 143.546 (cd^^). 



Scargle periodogram is constructed. 



B.1.9 Part 9 

Finally, the user is asked if folded light-curves should be produced. A list of all 
frequencies found in previous steps is plotted to screen, and the user is asked to 
choose which one to fold and plot. 



Do you want to fold the light curve? (answer yes or no) 
» yes 

Choose peak from the list below: 

Peak has a frequency of 3.596 (c/d) and a period of 0.278 (d) 
Peak 1 has a frequency of 143.546 (c/d) and a period of 0.007 (d) 
Peak 2 has a frequency of 36.586 (c/d) and a period of 0.027 (d) 
Which peak number do you want to fold? Type q when you want to quit. 



»1 

Folded light curves can be produced at several frequencies and the loop stops when 
the user types q. Figure IB.6I shows the folded light curve with a fitted sine wave 
plotted on top, according to the example above. 
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B.2 Analysis of Multiple Light Curves 



When running the program Speriod for combined light curves from several nights, 
I quickly realised that the program needed to be adapted to handle the longer data 
sets and account for the problems associated with aliases. I therefore developed 
Speriod-All for the analysis of more complex light curves. The input data file 
should have 3 columns; time, magnitudes or fluxes and errors. In addition to using 
Monte Carlo simulations to find frequencies with power excess and their associ- 
ated errors, bootstrap simulations are also performed as a comparison. Bootstrap 
simulations select a random subset of the light curve containing approximately 2/3 
of the the data values, and are useful for dealing with multiple-day datasets as the 
random re-sampling helps with reduction of aliases caused by the day-to-day data 
gaps. 

B.2.1 Part 1 

In part 1, a Lomb-Scargle periodogram is constructed from the input data file and 
the user is asked to choose which of the frequency regions should be investigated 
further. The user can zoom in and out in the Lomb-Scargle periodogram, and 
roughly define the left and right edges of every peak (Figure IB.2.HI . 

I got 4 options for you: 

1) If you want to zoom in on a peak, type 1 

2) If you want to zoom out to original size, type 2 

3) If you want to enter a range for the bootstrap analysis, type 3 

4) If you are done, type 4 
»1 

» 130,180 



B.2.2 Part 2 

Next, the user will be asked if Monte Carlo and bootstrap simulations should be 
performed. The oversampling is set to 64. Three plots are produced: 1) The ini- 
tial histogram of all peaks plotted over each user-defined frequency range. 2) The 
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(a) 



(b) 



Figure B.7: Figure^J(a)\shows the power spectrum for multiple nights, and Fi<iure \K7(b) 



shows a zoom of the same power spectrum in the range between 130 cd ^ - 180 cd ^ 
(according to the example in the text). 



highest alias found by the Monte Carlo simulations. 3)The highest alias found by 
the Bootstrap simulations (see Figure IB^ . Errors on the peaks are found from the 
output histograms of the peak frequencies found by the simulations. 

B.2.3 Part 3 

In part 3, a spectral cleaning is performed in the same way as described in the 
previous section (Figure |B31 |. 

B.2.4 Part 4 

Next, the spectral window functions are calculated as described in the previous 
section (Figure IB^ . 

B.2.5 Part 5 

The folding of the light curve is also performed in the same manner as described in 
the previous section. Finally, all results are printed to screen and to an output log 
file: 
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Figure B.8: The upper row shows a histogram of two user defined ranges where simulations 
will be performed. In the middle row, the highest alias found in the Monte Carlo simulations 
is plotted for each of the two user defined ranges. The bottom row shows the highest alias 
found in the bootstrap simulations. 
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OUTPUT - SPERIOD-ALL 
************************************************** 

Input file: mags.txt has 1876 data points 
Peaks are chosen at: 

Peak no. Is In range 0.0 - 10.0 c/d. 
Peak no. 1 Is In range 25.0 - 35.0 c/d. 

Bootstrap analysis gives... 
Statistics : 

Period: 5.06267578 +/- 0.00951901 c/d (1 slgma) 
Period: 28.79924526 +/- 0.01006367 c/d (1 slgma) 
Gausslanflt : 

Period: 5.06195725 +/- 0.01053709 c/d (1 slgma) 
Period: 28.79913552 +/- 0.01005182 c/d (1 slgma) 

Monte Carlo Simulations gives... 
Statistics : 

Period: 5.06213765 +/- 0.03155680 c/d (1 slgma) 
Period: 28.80118285 +/- 0.03214551 c/d (1 slgma) 
Gausslanflt : 

Period: 5.05995243 +/- 0.02726467 c/d (1 slgma) 
Period: 28.80196939 +/- 0.02545825 c/d (1 slgma) 

Good Bye! 
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Frequency (c/d) 



Figure B.9: Window function at frequency « 145 cd ^. 
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